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CONTEÚDO 


Capítulo 1 


Conceitos Fundamentais — Vetores 


Em qualquer teoria científica, e em mecânica em particular, é necessário começar com certos 
conceitos primitivos. Necessitaremos também fazer um certo número de suposições razoáveis. 
Dois dos conceitos mais básicos são espaço e tempo. No nosso estudo inicial da ciência 
do movimento, mecânica, vamos supor que o espaço físico das experiências comuns é des- 
crito adequadamente pelo espaço matemático tridimensional da geometria euclidiana. E em 
relação ao conceito de tempo, vamos supor que uma sequência ordenada de acontecimentos 
pode ser medida numa escala de tempo absoluta e uniforme. Vamos supor ainda que o espaço 
e tempo são entidades distintas e independentes. Mais tarde, quando estudarmos a teoria 
da relatividade, reexaminaremos os conceitos de espaço e de tempo e veremos que eles não 
são independentes e nem absolutos. Porém, isto é um assunto ao qual retornaremos depois 
de estudarmos os fundamentos clássicos da mecânica. 

Para definir a posição de um corpo no espaço, é necessário ter um sistema de referência. 
Em mecânica, usamos um sistema de coordenadas. O tipo básico de sistema de coordenadas 
que satisfaz aos nossos propósitos é o Sistema de Coordenadas Cartesianas ou retangulares, 
um conjunto de três linhas retas ou eixos mutuamente perpendiculares. Especificaremos a 
posição de um ponto em tal sistema através de três números ou coordenadas, x, y, e z. As 
coordenadas de um ponto móvel mudam com o tempo, isto é, elas são funções da quantidade 
t medida na nossa escala de tempo. 

Um conceito muito útil em mecânica é o de partícula ou ponto de massa, uma entidade 
que tem massa (o conceito de massa será discutido no capítulo 2), mas não tem extensão 
espacial. Rigorosamente falando a partícula é uma idealização que não existe — mesmo um 
elétron tem um tamanho finito — mas a ideia é útil como uma aproximação de um corpo 
pequeno, ou seja, um corpo cujo tamanho seja relativamente sem importância numa discussão 
particular. A Terra, por exemplo, pode ser tratada como uma partícula em mecânica celeste. 


1.1 Grandezas Físicas e Unidades 


Os resultados experimentais da física são expressos em termos de certas entidades fundamen- 
tais chamadas grandezas físicas — por exemplo, comprimento, tempo, força e outras. Uma 
grandeza física é algo que pode ser medido quantitativamente em relação a alguma unidade 
escolhida. Quando dizemos que o comprimento de um certo objeto é, digamos 7 polegadas, 
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queremos dizer que a medida quantitativa 7 é a relação (razão) do comprimento daquele 
objeto para o comprimento da unidade (1 polegada). Observamos que é possível definir 
todas as unidades das grandezas físicas da mecânica em termos de apenas três unidades 
básicas, ou seja comprimento, massa e tempo. 


A UNIDADE DE COMPRIMENTO 


A unidade padrão de comprimento é o metro. O metro era anteriormente definido como 
a distância entre duas marcas numa barra de platina guardada no Museu Internacional de 
Padrões Métricos, Sevres, França. Atualmente definimos o metro como a distância ocu- 
pada por exatamente 1.650.763,73 comprimentos de onda de luz da linha laranja do isótopo 
Kriptonio 86. 


A UNIDADE DE MASSA 


A unidade padrão de massa é o quilograma. O quilograma é a massa de um bloco de 
platina-irídio também guardada no Museu Internacional. 


A UNIDADE DE TEMPO 


A unidade básica para a medida de tempo, o segundo, era anteriormente definida em 
termos da rotação da Terra. Mas, como o metro, o segundo é atualmente definido em termos 
de um padrão atômico específico. O segundo é, por definição, o intervalo de tempo necessário 
para exatamente 9.192.631,770 oscilações de uma transição atômica particular do isótopo do 
Césio de número de massa 133. 


O sistema de unidades acima é chamado de sistema MKS. (Neste sistema existe uma 
quarta unidade, o Coulomb, que é usado para definir unidades elétricas). Os padrões 
atômicos modernos de comprimento e tempo neste sistema não são apenas mais precisos 
do que os padrões anteriores, mas são também universalmente reproduzíveis e indestrutíveis. 
Infelizmente, ainda não é no presente tecnicamente factível empregar um padrão atômico de 
massa. 

Na verdade, não existe nada particularmente especial a respeito das quantidades físicas 
comprimento, massa e tempo como um conjunto básico para definir unidades. Outros con- 
juntos de quantidades físicas podem ser usados. Os chamados sistemas gravitacionais usam 
comprimento, força e tempo. 

Além do sistema MKS, existem outros sistemas em uso comum, ou seja, o CGS, ou 
sistema centímetro-grama-segundo, e o PLS, ou pé-libra-segundo. Estes dois últimos sistemas 
podem ser vistos como secundários com relação ao sistema MKS porque suas unidades são 
especificamente frações definidas das unidades MKS: 
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1 cm = 107? m 

1 g = 10™° kg 

1 pé = 0,3048 m 

1 libra = 0,4536 kg 


1.2 Grandezas Escalares e Vetoriais 


Uma grandeza física que fica completamente especificada por um único número é chamada 
um escalar. Exemplos familiares de escalares são densidade, volume e temperatura. Mate- 
maticamente, escalares são tratados como números reais comuns. Eles obedecem a todas as 
regras algébricas regulares de adição, subtração, multiplicação, divisão e assim por diante. 

Existem certas grandezas físicas que possuem uma característica direcional, como um 
deslocamento de um ponto para outro do espaço. Tais grandezas necessitam uma direção, 
sentido e um módulo para as suas completas especificações. Estas grandezas são chamadas 
vetores se elas obedecem a regra da adição do paralelogramo como será discutido na seção 
1.7. (Um exemplo de uma grandeza direcional que não obedece a regra para adição é uma 
rotação finita de um objeto em torno de um dado eixo. O leitor pode prontamente verificar 
que duas rotações sucessivas em torno de eixos diferentes não produzem o mesmo efeito 
que uma única rotação determinada pela regra do paralelogramo. Por ora, todavia não 
examinaremos tais grandezas direcionais não vetoriais). Além de deslocamento no espaço, 
outros exemplos familiares de vetores são velocidade, aceleração e força. O conceito de vetor 
e o desenvolvimento de toda uma matemática de quantidades vetoriais tem sido indispensável 
ao desenvolvimento da ciência da Mecânica. O resto deste capítulo será em sua maior parte 
voltada a um estudo da matemática de vetores. 


1.3 Notação 


Quantidades vetoriais são representadas em imprensa por tipo em negrito, por exemplo A, 
enquanto tipo itálico representa quantidades escalares. Em trabalho manuscrito é costume 
usar uma flecha, A, para representar um vetor. 

Especificamos um vetor A por seu módulo, sua direção e sentido em relação a algum sis- 
tema de referência escolhido. Diagramaticamente, representamos um vetor por um segmento 
de linha direcionada, como mostrado na Figura 1.1. 

Um vetor pode também ser especificado relacionando-se suas componentes ou proje- 
ções sobre os eixos coordenados. O símbolo de componentes [A,, Ay, Az] será usado como 
uma representação alternativa de um vetor. O lado direito da equação A = [ As, Ay, Az] 
exprime o vetor A em termos de suas componentes num sistema de coordenadas particular. 
(Consideremos subentendido o uso do sistema de coordenadas cartesianas, a menos que seja 
dito o contrário). Por exemplo, se o vetor A representa um deslocamento de um ponto 
Pi(x1,y,21) até o ponto P(x2,y2,22), então A, = £2 — £1, Åy = Y2 — Y1, Áz = 22 — 21. Se A 
representa uma força então A, é a componente x da força, e assim por diante. Evidentemente, 
os valores numéricos das componentes escalares de um dado vetor dependem da escolha dos 
eixos coordenados. 
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xX 


Figura 1.1: Componentes de um vetor em coordenadas cartesianas. 


Restringindo a discussão a vetores contidos num plano, somente duas componentes são 
necessárias. Por outro lado, podemos definir um espaço matemático de qualquer número de 
dimensões. Então o símbolo [A;, As, As, ..., An] representa um vetor n-dimensional. Neste 
sentido abstrato um vetor é um conjunto ordenado de números. 


1.4 Definições Formais e Regras 
Começamos o estudo da álgebra vetorial com algumas definições formais sobre vetores. 


1. Igualdade de Vetores 


A equação À = B ou 
[As Ay, A.] = [Bx, B}, B.] 


é equivalente às três equações 


Isto é, dois vetores são iguais se, e somente se, suas respectivas componentes forem iguais. 


2. Adição de Vetores 


A adição de dois vetores é definida pela equação 
A +B = [As, Ay, 4.] + [Bo, By, B.] = [Ar + Br, Ay + By, A. + B.] 


A soma de dois vetores é um vetor cujas componentes são as somas das componentes dos 
vetores dados. 


3. Multiplicação por um Escalar 
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Se c é um escalar e À é um vetor, 
cÀ = ClAs, Ay, A2] = [CAs CA, CA,| = Ac 


O produto cA é um vetor cujas componentes são c vezes aquelas de A. 


4. Subtração de Vetores 


Subtração é definida como se segue: 


AB Ape nB=|A=-BaA;-BpAz=B/ 


5. O Vetor Nulo 
O vetor O = [0, 0, 0] é chamado de vetor nulo. A direção e sentido do vetor nulo não são 
definidas. De (4) segue-se que A — A = 0. Desde que não haja confusão quando o vetor nulo 
for representado por um “zero”, usaremos daqui para a frente a notação O = 0. 
6. A Lei Comutativa da Adição 
Esta lei é válida para vetores; isto é, 
A+B=B+4A 


desde que A, + B, = B, + Az, e da mesma forma para as componentes y e 2. 


7. A Lei Associativa 


A lei associativa é também verdadeira, porque 


A+(B+O = [A +(B. +05), 4,+(B,+C,), 4, +(B. +04) 
= (4, + Bz) + Cr, (Ay + By) + Cy, (A: + B2) + C] 
(A+B+C 


8. A Lei Distributiva 


A lei distributiva é válida na multiplicação de um vetor por um escalar porque, de (2) e 


(3) 


d(A+B) = dA, +B,A + Bn AEB] 

= [(A + Bi), (A, + By), (A. + B.)] 
CA + cBr,cAy + cBy, cA; 06] 
cĂ + cB 


Então os vetores obedecem as regras da álgebra ordinária enquanto estivermos conside- 
rando as operações acima. 
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1.5 Módulo de um Vetor 


O módulo de um vetor À, representado por IA] ou por 4, é definido como a raiz quadrada 
da soma dos quadrados das componentes, ou seja, 


A=|A|=,/A2+A2+ A? (1.1) 


Geometricamente o módulo de um vetor é seu comprimento, isto é, o comprimento da dia- 
gonal do paralelepípedo retangular cujos lados são A,, A, e Az. 


1.6 Vetores Unitários 


Um vetor unitário é um vetor cujo módulo é unitário. Vetores unitários são frequentemente 
representados pelo símbolo é da palavra alemã einheit. Os três vetores 


& = [1,0,0]  & = [0,1,0] &=[0,0,1] a2 


são chamados vetores unitários ou vetores base. Qualquer vetor pode ser expresso em termos 
de vetores base como uma soma vetorial de componentes como se segue: 


> 


AA AA) = ([A5,0,0]4-/0,4,,0] 4+[0,0,4,)] 
= A,[1,0,0]+ 4,[0, 1,0] + 4./0,0,1] 
= rAr +EyAy HEM, (1.3) 


Figura 1.2: Os vetores coordenados unitários TĪ k. 


Uma notação amplamente usada para vetores unitários Cartesianos são as letras T, J e k, ou 
seja 


ree J-e kea 


Daqui para a frente usaremos sempre essa notação. 


1.7. SIGNIFICADO GEOMÉTRICO DAS OPERAÇÕES VETORIAIS T 


As direções dos vetores unitários são definidas pelos eixos coordenados (Figura 1.2). Eles 
formam uma tríade orientada pela regra da mão direita ou pela mão esquerda, dependendo 
de qual sistema de coordenadas é usado. É costume usar sistemas coordenados orientado de 
acordo com a mão direita. O sistema mostrado na Figura 1.2 é orientado de acordo com a 
mão direita. 


1.7 Significado Geométrico das Operações Vetoriais 


Se representarmos um vetor por um segmento de linha orientado, interpretaremos facilmente 
as definições enunciadas acima como se segue: 


1. Igualdade de Vetores 


Se dois vetores são iguais então os vetores são paralelos e têm o mesmo comprimento, 
mas não têm necessariamente a mesma posição. A Figura 1.3, bidimensional para maior 
clareza, mostra dois vetores iguais. 


k A=B 
B À B 
A 
Ay 
= > 
O Ax Bx é 


Figura 1.3: Ilustrando igualdade de vetores. 


Observe que os vetores formam lados opostos de um paralelogramo. (Vetores iguais não 
são necessariamente equivalentes em todos os casos. Então duas forças vetorialmente iguais 
atuando em diferentes pontos de um objeto podem produzir diferentes efeitos mecânicos). 


2. Adição de Vetores 


A soma vetorial de dois vetores é igual ao terceiro lado de um triângulo, dois lados do 
qual são os vetores dados. A soma vetorial está ilustrada na Figura 1.4. A soma é dada 
também pela regra do paralelogramo, como mostrado na Figura. (A soma vetorial é definida 
de acordo com a definição dada na seção 1.4 (2), mesmo se os vetores não têm um ponto 
comum). 


3. Multiplicação de um Vetor por um Escalar 
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C=A+B=B+A 


Figura 1.4: Adição de dois vetores. 


O vetor cÀ é paralelo a À e tem c vezes o comprimento de A. Quando c = —1, o vetor 
—A tem a mesma direção e sentido oposto a A, como mostrado na Figura 1.5. 


Figura 1.5: O negativo de um vetor. 


1.8 O Produto Escalar 


Dado dois vetores A e Bo produto escalar ou produto interno, A. B, é o escalar definido 
pela equação 


A.B-AB,+A,B,+4,B, (1.4) 


Segue-se dessa definição que 


A.B=B.A (1.5) 


Desde que 4,B, = By Ay e assim por diante. Segue-se também que 


> 


A (B+O=ÀÃ.B+Ã. E (1.6) 
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porque se aplicarmos a definição (1.4) em detalhes 


A (BFO = AMB. CAMBRA B AC) 
= AB +4,B,+4,B, + AC + AO, + AL, 
As BÃO 


Da geometria analítica lembramos a fórmula para o cosseno do ângulo entre dois segmentos 


de linha 
ABA Boss 


VA + AZ + 42, B? + B2 + B2 


cos ĝ = 


Usando as equações (1.1) e (1.4), a fôrmula acima pode ser escrita 


MESE: 
cos =» 
ou e 
A . B = AB cos0 (1.7) 


A equação acima | pode ser vista como uma definição alternativa do produto escalar. Geome- 
tricamente, A. B é igual ao comprimento da projeção de À sobre B, vezes o comprimento 
de B. 

Se o produto escalar A. Bé nulo, então A é perpendicular a B, desde que nem A nem 
B sejam nulos. 

O quadrado do módulo de um vetor A é dado pelo produto escalar de A com ele mesmo, 


> 


= |P=. A 
Das definições dos vetores unitários T, J e k, fica clara a validade das seguintes relações 


= 1 
=0 (1.8) 


T 
T 


1.9 Alguns Exemplos do Produto Escalar 


1. Componente de um Vetor. Trabalho 


Como exemplo de um produto escalar, suponha que um objeto sob a ação de uma força 
constante (o conceito de força será discutido mais tarde no Capítulo 2), sofre um desloca- 
mento linear AS, como é mostrado na Figura 1.6. Por definição, o trabalho AW realizado 
pela força é dado pelo produto da componente da força F na direção AS, multiplicado pelo 
módulo AS do deslocamento, isto é, 


AW = (Fcos0)AS 
onde 6 é o ângulo entre Fe AS, isto é, 


AW =F. AS 
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9) 
AS 


Figura 1.6: Uma força submetida a um deslocamento. 


2. Lei dos Cossenos 


Consideremos um triângulo cujos lados são A, B e C, como mostrado na Figura 1.7. 


Então Č = A + B. Vamos fazer o produto escalar de Č com ele mesmo 
C.6 = (A+B).(A+8B) 
= ALA+2A.B+B.B 


A segunda passagem segue-se da aplicação das regras nas equações (1.5) e (1.6). Vamos 
substituir A . B por AB cos 0 para obter 


C? = A? + B? + 2AB cos 0 


que é a conhecida lei dos cossenos. Este é apenas um dos exemplos do uso da álgebra vetorial 
para provar teoremas em geometria. 


Figura 1.7: A lei dos cossenos. 


1.10 O Produto Vetorial 


Dado dois vetores À e B, o produto vetorial A x B, é definido pelo vetor cujas componentes 
são dadas pela equação 


AxB=[A,B, — A,B,, AB — AB, AB, — A Bil (1.9) 
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A interpretação geométrica do produto vetorial é dada na seção 1.11. Podemos mostrar que 
as seguintes regras se aplicam para a multiplicação vetorial. 


AxB=-BxÃ (1.10) 
Ax(B+0O)-AxB+AÃx (1.11) 
n(AxB)=n(A)xB=Axn(B) (1.12) 


As provas destas relações decorrem diretamente das definições e vamos deixá-las como exer- 
cício. 

De acordo com as definições algébricas dos vetores unitários, Equação (1.2), verificamos 
imediatamente que as seguintes relações para os produtos vetoriais são verdadeiras: 


Ixk=T=-kx] (1.13) 


Por exemplo 


> 


Tx = [0 -0,0—0,1 -— 0] = [0,0,1] = k 


De maneira idêntica verificamos facilmente as outras relações. 


1.11 Interpretação Geométrica do Produto Vetorial 


> 


O produto vetorial expresso em termos dos unitários TJ k é 


> 


AxB=TA,B,— AB +J(A.B, — AB.) + RAB; — AB) 


cada termo entre parênteses é igual a um determinante 


OE TE E a ro cal id A! caio id 
SB Bo) PIE, nb ag B, 
e finalmente A 
a e US E E 
ASBES ÃO o A (1.14) 
B, B, B, 


que é prontamente verificada pelo desenvolvimento. 

A forma determinante é uma ajuda conveniente para se lembrar a definição < de produto 
vetorial. Das propriedades de determinante pode-se ver que se A é paralelo a B, isto é, se 
A = cB então as duas últimas linhas do determinante são proporcionais e logo o determinante 
é nulo. Então o produto vetorial de dois vetores paralelos é nulo. 

Vamos calcular o módulo do produto vetorial. Temos 


À x BP = (A,B, — A;B)? + (AB, — AB) + (AB, — A,B)? 
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Com um pouco de paciência podemos reduzir esta expressão a 
JA x BP = (4? + 4? + A2)(B2+ B? + B?) — (ABe + AyBy + A-B}? 
ou usando a definição de produto escalar, a equação acima pode ser colocada na forma 
|A x B? = AB? — (A . BP 


Tomando a raiz quadrada de ambos os lados da equação acima e usando a Equação (1.7), 
podemos expressar o módulo do produto vetorial como 


|Ā x B| = ABv1 = cos? 0 = AB sen 0 (1.15) 


onde 6 é o ângulo entre AeB. 


Figura 1.8: O produto vetorial. 


Para interpretar geometricamente o produto vetorial, observamos que o vetor C=AxB 
é perpendicular a ambos Ae B, pois 


Da AC A CAE: 
= AAB,- A,B) FA (ABe -ABA FAA By -ABa 
= 0 


Da mesma forma B . Č = 0. Então o vetor Č é perpendicular ao plano que contém os vetores 
A eB. 

O sentido do vetor Č = À x B é determinado a partir da condição de que os três vetores 
A, Be Č formem uma tríade “orientada de acordo com a mão direita”, como mostrado na 
Figura 1.8. (Isto está consistente com o resultado previamente estabelecido de que na tríade 
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>, 


orientada de acordo com a mão direita TTK temos Tx J = k). Consequentemente, da equação 
(1.15) vemos que podemos escrever 


A x B = (AB sen 6) (1.16) 


Onde r é um vetor unitário normal ao plano dos dois vetores A e B. O sentido de à é dado 
pela regra da mão direita, isto é, a direção de avanço de um parafuso de rosca direita girando 
da posição positiva de A para a de B através do menor ângulo entre eles, como está ilustrado 
na Figura 1.8. A Equação (1.16) pode ser entendida como uma definição alternativa do 
produto vetorial. 


Exemplos 


> 


— k, B =1—)J+ 2k, encontre À. B e A x B. 


= 
a 
+ 
= 


1. Dados os dois vetores 


AB dr) es 
zilestr aeir 


2. Encontre o ângulo entre A eB. Temos, da definição de produto escalar, 


> 


Abe -10 1 
AB /2+2+( E 2+22 “OVO 6 


Logo 9 = arccos(— 1) = 99,6º. 


cos 0 = 


1.12 Um Exemplo do Produto Vetorial. Momento de 
uma Força 


A representação dos momentos constitui uma aplicação útil do produto vetorial. Considere- 
mos uma força F atuando em um ponto P(x, y,z), como mostrado na Figura 1.9, e vamos 
representar o vetor OP por T, isto é 


o momento N, ou o Torque, em relação a um dado ponto O é definido como o produto vetorial 
N=rxF (1.17) 


Então o momento de uma força em relação a um ponto é uma quantidade vetorial tendo 
módulo, direção e sentido. Se uma única força for aplicada em um ponto P de um corpo 
livre para girar em torno de um ponto fixo O como um pivot, então o corpo tende a girar. O 
eixo desta rotação é perpendicular à força F, e é também perpendicular à linha OP. Dessa 
maneira a direção e sentido do torque será ao longo do eixo de rotação. 
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N=rxFr 


3 
P` 
g 
s 
ad) 
' 
; 
< 
3 
; 


Figura 1.9: O momento de uma Força. 


O módulo do torque é dado por 
IÑ| = |? x F| = rF sen 8 (1.18) 


onde 6 é o ângulo entre T e F. Logo INI pode ser visto como o produto do módulo da força 
pela quantidade r sen O que é exatamente a distância perpendicular da linha de ação da força 
ao ponto O. 

Quando várias forças são aplicadas a um único corpo em pontos diferentes os momentos 
se somam vetorialmente. Isto decorre da lei distributiva de multiplicação vetorial, Equação 
(1.11). A condição para equilíbrio de rotação é que a soma vetorial de todos os momentos 


seja zero: q 3 


Uma discussão mais completa disto será dado mais tarde no Capítulo 5. 
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1.13 Representação de um dado Vetor em Termos do 
Produto de um Escalar e um Unico Vetor Unitário 


Consideremos a equação , , 
A=TA+JA,+HKA, 

vamos multiplicar e dividir o lado direito pelo módulo de A 

A A, >A, ) 


AeA jk 
foda 


A 3 . > 
Agora + = cosa, | = cos f e = = cosy são os cossenos diretores do vetor A ea, 8 ey 


são os ângulos de direção. Então podemos escrever 
A= A(Tcos a + Jcos 8 + k cos Y) = Af[cos a, cos 8, cos y] 
ou E 
A= AÑ (1.19) 
onde t é um vetor unitário cujos componentes são cos a, cos f e cosy. 
Considere qualquer outro vetor B. Sem dúvida, a projeção de B em À é exatamente 
B.A > 
Bcosð = — =B. ñ (1.20) 


onde 0 é o ângulo entre AeB. 


Exemplo 


Encontre um vetor unitário perpendicular ao plano que contém os dois vetores 
A=W+]-keB=r1-J+2k 
Do exemplo 1, temos 


logo 
n es AxB _— i-5J-3k 
[AxB] — vI245243? 
ad T Jo sk 


1.14 Produtos Triplos 


A expressão OO C NS 

A.(BxC) 
é chamada triplo produto escalar de À, Be É. É um escalar desde que é o produto escalar 
de dois vetores. Lembrando-se da expressão em forma de determinante do produto vetorial, 
Equação (1.14), vemos que o triplo produto escalar pode ser escrito 


Dee a MÃE E do 
A. (Bx0)=|B. B, B (1.21) 
Gs Cos 
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Usando a conhecida propriedade de determinantes de que a troca de termos de duas linhas 
ou de duas colunas muda o sinal mas não muda o valor absoluto do determinante, podemos 
obter a seguinte expressão bastante útil: 


> > 


A. (BxO=(AxB). € (1.22) 
A expressão 
A x (Bx ©) 
é chamada Triplo Produto Vetorial. Deixamos o leitor provar que a seguinte equação é válida 
para o triplo produto vetorial 


A x (B x © = (A . ©B — (å . B)C (1.23) 


1.15 Mudança de Sistema de Coordenadas. A Matriz 
Transformação 


Consideremos o vetor À representado pela tríade TT k, 


> 


A= TAr At kA 


i ; FFI : Ros sar 
relativamente a uma nova tríade 1 7 k tendo orientação diferente daquela de Y7 k, o mesmo 
vetor A é representado por 

A = VA, SEJA + k' A, 
Mas o produto escalar A . 1 é exatamente Ax, isto é, a projeção de A no vetor unitário v. 
Então podemos escrever 


Ap=Ãi=(tNA+GDA+(k DA, 
Ay=A.} = (@. A+. FA +k. DA: (1.24) 
Ave hor e EPA Gelo (Œ. 1)A 


Os produtos escalares (T. 1), (T . J), e os demais, são chamados de coeficientes de transfor- 
mação. Eles são iguais aos cossenos diretores dos eixos do sistema de coordenadas com linha 
relativamente ao sistema sem linha. Podemos expressar as componentes sem linha da mesma 
forma como, 


> 


AE Te DAE AA A) 
AÃ.T= (ONA + . Ay +E . DA, (1.25) 
A=A.k = (f. K)Ay +0 DA + (KA, 


Todos os coeficientes de transformação na Equação (1.25) também aparecem na Equação 
(1.24), porque T. 1 = Y . T, etc., mas aqueles nas linhas (equações) das Equações (1.25) apa- 
recem nas colunas de termos nas Equações (1.24) e ao contrário. As regras de transformação 


1.15. MUDANÇA DE SISTEMA DE COORDENADAS 17 


expressas nestes dois conjuntos de equações são propriedades gerais de vetores. Na verdade, 
elas constituem uma maneira alternativa de se definir vetores!. 

As equações de transformação são convenientemente expressas em notação matricial. 
Então a Equação (1.24) é escrita 


Au e Jor E A 
Ap h= dogs kr Ay (1.26) 
Az DE J.K k. K Az 


A matriz 3 x 3 na equação acima é chamada a matriz transformação. Uma vantagem da 
notação matricial é que transformações sucessivas são prontamente executadas por meio de 
multiplicação de matrizes. 

O leitor observará que a aplicação de uma dada matriz transformação em um certo vetor 
A é também equivalente formalmente a girar aquele vetor dentro do sistema de coordenadas 
(fixo) sem linha, as componentes do vetor girado sendo aquelas dadas pela Equação (1.26). 
Então rotações finitas podem ser representadas por matrizes. (Observe que o sentido de 
rotação do vetor neste contexto é o oposto do sentido de rotação do sistema de coordenadas 
no contexto anterior). 


Exemplos 


1. Escreva o vetor A = 37 + 27 + k em termos da tríade 1'7k' onde os eixos x'y' estão 
girados de 45º em torno do eixo z. z e z” são eixos coincidentes, como mostra a Figura 1.10. 
Baseando-nos na figura, temos para os coeficientes de transformação 


T.l=1/V2 T.í=1/V2 k.í=0 
rJe qa kg 
r.k=0 T.k=0 k.k=1 
Destas relações tiramos 
3 2 5 3 2 1 
Ax = + = An == + = mr As=1 
v2 v2 42 Co ya v20 42 
Então, no sistema linha, o vetor A é dado por 
z 5o. a 
A = —' — —J' + k’ 
VET 


2. Encontre a matriz transformação para uma rotação do sistema de coordenadas linha de 
um ângulo q em torno do eixo z. 
(O exemplo anterior é um caso especial deste). Temos 


T. =J. J = coso 
J.Y = -T.J = seno 
k.k=1 


1Veja, por exemplo, L.P. Smith, Mathematical Methods for Scientists and Engineers, Prentice Hall, 
Englewood Cliffs, N.J., 1953. 
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Figura 1.10: O eixo linha O'x'y'z' estão girados de 45º em torno do eixo 2. 


e todos os outros produtos são nulos. Então a matriz transformação é 


coso seno 0 
—senġ coso 0 
0 0 1 


E evidente a partir do exemplo acima que a matriz transformação para a rotação em torno 
de um eixo coordenado diferente, digamos em torno do eixo y de um ângulo 0, será dado 
pela matriz 


cos6 0 —send 


0 1 0 
sen6 O cos6 


Consequentemente a matriz para a combinação de duas rotações, a primeira sendo em torno 
do eixo z (ângulo &) e a segunda sendo em torno do novo eixo y' (ângulo 0) é dada pelo 
produto matricial 


cos6 0 —sen6 cosy seny 0 cosĝ coso cosbsend — sen 
0 1 0 —senġ coso 0 |=| —senQ cos q 0 
send O cos 0 0 1 senĝcosø senĝseno cos 


A multiplicação matricial é, em geral, não comutativa. Então podemos esperar que se 
a ordem das rotações forem trocadas, e dessa maneira a ordem da multiplicação matricial 
na esquerda, o resultado final será diferente. Isto é verdade e o leitor poderá verificar. Isto 
está de acordo com uma observação feita anteriormente, ou seja, de que rotações finitas não 
obedecem a lei de adição vetorial e dessa maneira não são vetores apesar de uma rotação 
única ter uma direção e sentido (o eixo) e um módulo (o ângulo da rotação). Todavia, 
mostraremos mais tarde que rotações infinitesimais obedecem a lei da adição vetorial, e 
podem ser representadas por vetores. 
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1.16 Derivada de um Vetor 


Consideremos um vetor À cujas componentes são funções de uma única variável u. O vetor 
pode representar posição, velocidade, e assim por diante. O parâmetro u é normalmente o 
tempo t, mas pode ser qualquer quantidade que determina as componentes de A: 


A =TA,(u) + 7A, (u) + KA,(u) 


A derivada de A com respeito a u é definida, de uma maneira bastante análoga à derivada 
ordinária de uma função escalar, pelo limite 


dA — AÀ E agr 


Pa Au FJ Au EE 


onde AA, = A,(u + Au) — A,(u), e assim por diante. Logo 


dÀ dA, | (dA, | pç dA: 


o ee e du 


1.2 
du du (4:20) 


A derivada de um vetor, desse modo, é um vetor cujas componentes são derivadas ordinárias. 
Segue-se das equações acima que a derivada de uma soma de dois vetores é igual à soma 
das derivadas, ou seja 
div a dA dB 


Ab a d (1.28) 


Regras para diferenciar produtos vetoriais serão tratadas mais tarde na Seção 1.22. 


1.17 Vetor Posição de uma Partícula 


Em um dado sistema de referência a posição de uma partícula pode ser especificada por um 
único vetor, ou seja, o deslocamento da partícula em relação à origem do sistema de coorde- 
nadas. Este vetor é chamado de vetor posição da partícula. Em coordenadas retangulares, 
Figura 1.11, o vetor posição é simplesmente 


F=TIz+Jy+ks 
As componentes do vetor posição de uma partícula em movimento são funções do tempo, ou 


seja 
c=a(t)  y=y(t) z = z(t) 


1.18 O Vetor Velocidade 


Na Equação (1.27) demos a definição formal da derivada de qualquer vetor em relação a 
algum parâmetro. Em particular, se o vetor for o vetor posição r de uma partícula em 


20 CAPÍTULO 1. CONCEITOS FUNDAMENTAIS 


Figura 1.11: O vetor posição. 


movimento e o parâmetro for o tempo t, a derivada de r relativamente a t será chamada de 
velocidade, que representaremos por v. Logo 
dr 


=z Ii+TA kz (1.29) 


onde os pontos indicam diferenciação em relação a t. (Esta convenção é padrão e será 
usada doravante neste livro). Vamos examinar o significado geométrico do vetor velocidade. 
Suponha que uma partícula esteja em uma certa posição no instante t. Em um instante At 
mais tarde ela terá se movido da posição T(t) para a posição T(t + At). O vetor deslocamento 
durante o intervalo de tempo At é 


AF = (t+ At) — P(t) 


de modo que o quociente Ar/At é um vetor que é paralelo ao deslocamento. A proporção 
que consideramos intervalos de tempo cada vez menores, o quociente Ar/At se aproxima 
de um limite dr/dt que chamamos velocidade. O vetor dr/dt expressa a direção, o sentido 
e o módulo da velocidade. A Figura 1.12 mostra isto graficamente. No intervalo de tempo 
At a partícula move-se ao longo da trajetória de P para P’. Quando At se aproxima de 
zero, o ponto P’ se aproxima de P, e a direção do vetor Ar/At se aproxima da direção da 
tangente à trajetória em P. O vetor velocidade, desse modo, é sempre tangente à trajetória 
do movimento. 

Chamamos o módulo do vetor velocidade simplesmente de velocidade. Em componentes 


retangulares a velocidade é 
v = |V] = yt? + Y? ++ 2? (1.30) 


Se chamarmos a distância escalar ao longo da trajetória de s, então podemos alternativa- 
mente expressar a velocidade como 


ds i As a (Aa)? + (Ay)? + (Az)? 


Ei dE a URL T AO At 


que se reduz à expressão da direita da Equação (1.30). 
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* linha tangente 
ao ponto P 


Figura 1.12: Vetor deslocamento de uma partícula em movimento. 


1.19 Vetor Aceleração 


A derivada da velocidade em relação ao tempo é chamada de aceleração. Representando a 
aceleração por à, temos 


dř dr 
= — = — 1.31 
“O q de Edo 
Em componentes retangulares 
a=TE+Tj+ks (1.32) 


Então aceleração é uma quantidade vetorial cujas componentes, em coordenadas retangula- 
res, são as derivadas segundas das coordenadas de posição de uma partícula em movimento. 
A representação de à em componentes tangencial e normal será discutido na Seção 1.23. 


Exemplos 


1. Vamos examinar o movimento representado pela equação 
F(t) = Tbt + (ct — 2) + KO 


que representa um movimento no plano xy, desde que a componente z é constante e igual a 
zero. À velocidade v é obtida diferenciando-se em relação ao tempo, ou seja, 


A aceleração, do mesmo modo, é dada por 


E dy 
Qu Fê 
dt Jg 
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Então a é no sentido negativo da direção y e tem magnitude constante g. A trajetória 
do movimento é uma parábola, como mostrado na Figura 1.13. (Esta equação realmente 
representa o movimento de um projétil). A velocidade v varia com t de acordo com a equação 


v =b + (c — gt) 


0 x 


Figura 1.13: Vetores posição, velocidade e aceleração de uma partícula movendo-se em uma 
trajetória parabólica. 
2. Suponha que o vetor posição de uma partícula seja dado por 

r= Tbsen wt + Jb coswt + ke 


Vamos analisar o movimento. A distância à origem permanece constante 


IF| = r = Vb? sen? wt + b? cos? wt + @ = V + @ 
Diferenciando r, encontramos 


dr > 
V= Pi Tbw coswt — Jbw sen wt + kO 
Desde que a componente z de v é zero, o vetor velocidade é paralelo ao plano xy. A partícula 
se movimenta em sua trajetória com velocidade constante 


v = |7|= vb2w2 cos? wt + b2w2 sen? wt = bw 


A aceleração 


z= — = —Tbw’ sen wt — Jbw? cos wt 


é perpendicular à velocidade, pois o produto escalar de v e à se anula 


> 


V. E= (bw coswt)(—bw? sen wt) + (—bw sen wt)(—bw? coswt) = 0 


Além disso, a aceleração é perpendicular ao eixo z, como mostrado na Figura 1.14, porque 


a.k = 0. A trajetória real é um círculo de raio b contido no plano z = c. O movimento 
está ilustrado na Figura 1.14. 
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xX 


Figura 1.14: Vetores posição, velocidade e aceleração de uma partícula movendo-se em um 
círculo. 


1.20 Integração Vetorial 


Consideremos o caso em que a derivada em relação ao tempo de um vetor Y seja dada em 
coordenadas retangulares onde cada componente é uma função conhecida do tempo, ou seja, 


T =T filt) +7 f(t) + k fa(t) 


É possível integrarmos em relação ao tempo para obter 
aL T f hdt +T | HOt ER festejar (1.33) 


Trata-se evidentemente do inverso do processo de encontrar o vetor velocidade, quando o 
vetor posição é dado como função do tempo. O mesmo se aplica ao caso em que a aceleração 
é dada como uma função do tempo e uma integração nos leva à velocidade. 


Exemplo 


O vetor velocidade de uma partícula em movimento é dado por 
->C 
V=IÃ+JBt+k F 


onde 4, B e C são constantes. Encontre T. 
Por integração, obtemos 


> f Cdt (e a 
r=1 | Adt+7 f Bidt+k [O =TAL+JBS +KClnt+ão 


O vetor Tọ é a constante de integração. 
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1.21 Velocidade Relativa 


Consideremos duas partículas cujos vetores posição são Yı e r2, respectivamente, como mos- 
trado na Figura 1.15. O deslocamento da segunda partícula com relação à primeira é a 
diferença r3 — rı que chamaremos T12. A velocidade da segunda partícula relativamente à 


primeira é portanto 
os dri2 = d(T2 — Ti) = 
Vig = = = V2 — V1 


dt dt 
que chamamos de velocidade relativa. Transpondo V1, temos 


V2 = Vı + V12 


para a velocidade da partícula 2 em relação a O em termos da velocidade da partícula 1 e 
da velocidade relativa. 


(b) 


Figura 1.15: Vetor posição relativa (a) e vetor velocidade relativa (b) de duas partículas. 


Devemos notar que o módulo da velocidade relativa de duas partículas não é o mesmo 
que a variação temporal da distância entre elas. Esta última quantidade é 


d 


Ce E 
gel = gT A] 


que é diferente, em geral, de |V12|. 


Exemplos 


1. Uma partícula move-se ao longo do eixo x com velocidade v, sendo seu vetor posição 
dado por 74 = T(a + vt) onde a é uma constante. Uma segunda partícula move-se ao longo 
do eixo y com mesma velocidade, e seu vetor posição é r> = J(b+ vt). Então a velocidade 
da segunda partícula em relação à primeira é 


V12 = Və — V1 =v -Tv =v(7 T) 
O módulo da velocidade relativa é então 
Vig = vv2. 


Qual é o valor de d|r12|/dt? 
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2. Um aro de raio b rola no chão com uma velocidade Vo para a frente. Encontre a velocidade 
em relação ao chão de um ponto P qualquer da periferia do aro. Primeiro considere a 
expressão 

rop = Tb cos 0 — Tb sen 0 


onde 
0 = wt 


Isto representa um movimento circular no sentido horário em torno da origem, o centro 
da roda, neste caso. A derivada em relação ao tempo então nos dá a velocidade de P 
relativamente ao centro do aro 


Vre = —Tbw sen 0 — Jbwcos 0 


Mas a velocidade angular é w = vo/b, e desde que a velocidade do centro do aro em relação 
ao solo é Tv, então a velocidade de P em relação ao solo é 


V = Tvo—Ibwsen6 — Tbw cos 0 
= Tvo(1 — sen 0) — J vo cos 0 


A Figura 1.16 mostra os vetores velocidades para vários valores de 6. 


| Sa pa Vel 


~x vwv=0 0 4 
w E y 


Figura 1.16: Vetores velocidade para vários pontos de um aro que rola. 


1.22 Derivadas de Produtos de Vetores 


Frequentemente é necessário tratar com derivadas dos produtos nA, A. Be A x B onde o 
escalar n e os vetores A e B são funções de um único parâmetro u, como na Seção 1.16. Da 
definição geral de derivada, temos 


d(nÃ) ipg n(u + Aw)Ã(u + Au) — n(u)A (u) 
du Au>0 Au 
d(A . B) A(u + Au) . B(u + Au) — A (u) . B(u) 


du Au—>0 Au 
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> > > 


— Atu+ Au) x B(u + Au) — A (u) x B(u) 
— = lim 
du Au>0 Au 


Adicionando-se e subtraindo-se expressões como n(u + Au)A(w) aos numeradores, obtemos 
as seguintes regras: 


d(nÃ) E dn x | dÄ 


i 1.34 
du du du (134) 
dA. B) dA > > dB 
= i | ; 1. 
du du Eres du (1:35) 
dAxB) dA = >» dB 
= B+A 1.36 
du du á s du ( 


Observemos que é necessário preservar a ordem dos termos na derivada do produto ve- 
torial. As etapas são deixadas como exercício para o leitor. 


1.23 Componentes Normal e Tangencial da Aceleração 


Mostrou-se na Seção 1.13 que qualquer vetor pode ser expresso como o produto de seu módulo 
por um vetor unitário de direção idêntica a este. Dessa maneira, o vetor velocidade de uma 
partícula em movimento pode ser escrito como o produto do módulo dessa velocidade v por 
um vetor unitário que dá a direção do movimento da partícula. Logo 


7=uvr (1.37) 


O vetor F é chamado de vetor unitário tangente. A proporção que a partícula se move o 

módulo da velocidade pode mudar e a direção de 7 pode mudar. Vamos usar a regra para 

diferenciação de produto de um escalar por um vetor para obter a aceleração. O resultado é 
dv dor) dr 


ES tt (1.38) 


O vetor unitário 7, sendo de módulo constante, tem uma derivada dT /dt que necessariamen- 
te deve expressar a mudança na direção de 7 com o tempo. Isto está ilustrado na Figura 
117(a). A partícula está inicialmente em um ponto P qualquer de sua trajetória. 

Em um intervalo de tempo At a partícula move-se para outro ponto P’ distante As 
medido ao longo da trajetória. Vamos chamar os vetores unitários tangentes em Pe P’ de 7 
e 7, respectivamente, como mostrado. As direções destes dois vetores diferem de um certo 
ângulo A¢ġ como mostrado na Figura 1.17(b). 

É fácil ver que para pequenos valores de Ag, a diferença AF se aproxima de Aq em 
grandeza. Também, a direção de AF torna-se perpendicular à direção 7 no limite quando 
Ag e As se aproximam de zero. Segue-se que a derivada dT /dy tem módulo unitário e é 
perpendicular a 7. Deveremos então chamá-lo de vetor unitário normal e representá-lo por 


Ema 


n 


dT 


25 = (1.39) 
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AT 


» LAS 


(a) (b) 


Figura 1.17: Vetores unitários tangente e normal. 


A seguir, para encontrar a derivada de F em relação ao tempo dT/dt, usamos a regra da 


cadeia como segue 
dr drdọo  „dọds z2 


dt dọdt “dsdi “p 


onde 
ds 


do 


é o raio de curvatura da trajetória da partícula no ponto P. O valor acima para dT/dt é 
agora levado na Equação (1.38) para nos dar o resultado final 


p 


2 
B (1.40) 


a= ÙT + 
Então a aceleração de uma partícula em movimento tem uma componente de módulo 


Ar UES 


na direção do movimento. Esta é a aceleração tangencial. A outra componente de módulo 
J2 
an = — 
p 


é a componente normal. Esta componente está sempre apontando diretamente para o cen- 
tro de curvatura no lado côncavo da trajetória. Por isso a componente normal é também 
chamada de aceleração centrípeta. 
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(a) (b) 


Figura 1.18: Vetores aceleração para uma partícula se movendo numa trajetória circular. 
(a) velocidade constante; (b) velocidade aumentando. 


Das considerações acima vemos que a derivada relativa ao tempo do módulo da velocidade 
é apenas a componente tangencial da aceleração. O módulo da aceleração total é dado por 


= [2 + (1.41) 


v 
dt 


Por exemplo, se uma partícula move-se num círculo com velocidade constante v, o vetor 
aceleração tem módulo a onde Ro é o raio do círculo. O vetor aceleração aponta sempre 
para o centro do círculo nesse caso. Todavia, se o módulo da velocidade não é constante 
mas aumenta numa taxa v, então a aceleração tem uma componente para a frente desta 
quantidade e é desviada do centro do círculo para o lado do movimento, como ilustrado na 
Figura 1.18. Se a partícula está parando, então o vetor aceleração é desviado na direção 


oposta. 


1.24 Velocidade e Aceleração em Coordenadas Polares 
Planas 
Muitas vezes é conveniente usar coordenadas polares (r, 0) para expressar a posição de uma 


partícula que se move em um plano. Vetorialmente, a posição da partícula pode ser escrita 
como o produto da distância radial r por um vetor unitário radial €,.: 


f=T6, (1.42) 


Quando a partícula se move, ambos r e €, variam, pois ambos são funções do tempo. Dessa 
maneira, se derivarmos em relação ao tempo, teremos 
dr dē, 
F= L = iE, Hra 1.43 
dt dt (E) 
Para calcular a derivada de, /dt vamos considerar o diagrama vetorial na Figura 1.19. Um 
estudo da figura mostra que quando a direção de r varia de uma quantidade Ad, a mudança 
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correspondente Ae, no vetor radial unitário será obtido da seguinte maneira: o módulo 
|Ae,| é aproximadamente igual a A0, e a direção de Ag, é quase perpendicular a €. Vamos 
introduzir outro vetor unitário € cuja direção é perpendicular a €.. Então temos 


Ae, = €&A0 
Dividindo-se por At e tomando-se o limite, obtemos 


dë, d0 


para a derivada em relação ao tempo do vetor radial unitário. Da mesma forma, podemos 
argumentar que a mudança no vetor € é dado pela aproximação 


Aē = —€, A0 


Aqui o sinal negativo é colocado para indicar que a direção da variação Ae, é oposto à 
direção de €, como pode ser visto na figura. Consequentemente, a derivada temporal de eg 
é dada por 


USD quo dt 1.4 
dt S o) 


Figura 1.19: Vetores unitários para coordenadas polar plana. 


Usando a Equação (1.44) para a derivada do vetor radial unitário, podemos finalmente 
escrever a equação para a velocidade como 


Y = 16, + re (1.46) 


Então ř é o valor da componente radial do vetor velocidade, e r0 é o valor da componente 
transversal. 
Para determinar o vetor aceleração, tomamos a derivada em relação ao tempo do vetor 
velocidade. Teremos 
dê r dês 


Fa (70 4 rojat Tô 
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Os valores de de, /dt e deg/dt são dados pelas Equações (1.44) e (1.45) e nos leva à seguinte 
equação para o vetor aceleração em coordenadas polares planas: 


a= (F —rô3e, + (rË + 270)6 (1.47) 


Então o valor da componente radial do vetor aceleração é 


a = —r6? (1.48) 
e a componente transversal é 
à + 200 = 1º (120) (1.49) 
dg =" r0=-—(r i 
f rdt 


O resultado acima mostra, por exemplo, que se uma partícula se move num círculo de raio 
constante b, então 7 = 0, e dessa maneira a componente radial da aceleração tem valor bô? e 
aponta diretamente para o centro da trajetória circular. O valor da componente transversal 
neste caso é b0. Por outro lado, se a partícula se move ao longo de uma linha radial fixa, 
isto é, se 0 é constante, então a componente radial se reduz a 7 e a componente transversal 
se anula. Se r e O ambos variam, então a expressão geral (1.47) dá a aceleração. 


Exemplo 
Uma partícula move-se em uma trajetória espiral dada pelas coordenadas polares 
r = bt? 0=ct 
onde b e c são constantes. Encontre a velocidade e aceleração como função de t. Da Equação 


(1.46), encontramos 


É (ct) = (2bt)6, + (bet?) cg 


(bt?) + ob) 


E 
Da mesma forma, da Equação (1.47), temos 
a = e.(2b — bte) + G9[0 + 2(2bt)c] = b(2 — tc?)6, + 4betêo 


é interessante notar que, neste exemplo, a componente radial da aceleração torna-se negativa 
para t grande, apesar do raio estar sempre crescendo monotonicamente com o tempo. 


1.25 Velocidade e Aceleração em Coordenadas Cilin- 
dricas e Esféricas 


1.25.1 Coordenadas Cilíndricas 


No caso de movimento tridimensional, a posição de uma partícula pode ser descrita em 
termos das coordenadas cilíndricas R, ø, z. O vetor posição pode então ser escrito na forma 


P= Repr + ze, (1.50) 
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onde ep é um vetor radial unitário no plano xy e é, é um vetor unitário na direção z. Um 
terceiro vetor unitário €% é necessário para que os três vetores Cree. constituam uma tríade 
orientada de acordo com a mão direita como ilustrado na Figura 1.20. Notamos que k=8,. 

Os vetores velocidade e aceleração são obtidos por diferenciação, como antes. Isto nova- 
mente envolverá derivadas de vetor unitários. Um argumento semelhante àquele usado no 
caso bidimensional mostra que der /dt = Esh e des/dt = —Er¢. O vetor unitário €, não varia 
sua direção, logo sua derivada relativa ao tempo é zero. 

Tendo em vista estes fatos, vemos facilmente que os vetores velocidade e aceleração são 
dados pelas relações abaixo 


V = REg + Roeg + 3, (1.51) 
ï = (R — Rọ°)Er + (CR$ + R$)E + žē: (1.52) 
Estas equações dão os valores de v e à em termos de suas componentes na tríade girada 
EREJE:z. 
z 
e, 
€p 
ep 
z 


Figura 1.20: Vetores unitários para coordenadas cilíndricas. 


Alternativamente, podemos obter as derivadas dos vetores unitários usando as seguintes 
equações que relacionam as tríades fixa e girada. 


Er = Tcoso+Iseno 
e = -Tseng +J coso (1.53) 
& =k 


Deixamos as etapas como um exercício. O mesmo resultado pode ser também obtido utili- 
zando a matriz rotação como foi feito no Exemplo 2, Seção 1.15. 
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1.25.2 Coordenadas Esféricas 


Quando utilizarmos as coordenadas esféricas r, 0, 4 para especificar a posição de uma partícula, 
escrevemos o vetor posição como o produto da distância radial r pelo vetor radial unitário 
€r da mesma maneira que foi feito quando do uso de coordenadas polares planas. Então 


T= rE, 


A direção €, fica agora especificada por dois ângulos 0 e à. Introduzimos mais dois vetores 
unitários ē e & como mostrado na Figura 1.21. 


Z 


€o 


Figura 1.21: Vetores unitários para coordenadas esféricas. 


A velocidade é dr dë 
T Er 
V= — = ie +r 1.54 
dt dt ER 
O nosso problema seguinte é expressar a derivada de, /dt em termos dos vetores unitários da 
tríade girada. 


Usando a figura, vemos que as seguintes relações entre as duas tríades são válidas. 


& = Tsen 0 coso +7 senl sen q + k cos0 
& = T cos 0 cos o + J cos Ô sen ġ — K sen 0 (1.55) 
€p = -Tseng +] coso 


Estas equações expressam os vetores unitários da tríade girada em termos da tríade fixa. 
Notemos a semelhança entre esta transformação e aquela da segunda parte do Exemplo 2 na 
Seção 1.15. Identificando as rotações corretamente vemos que as duas transformações são, 
de fato, idênticas. Vamos diferenciar a primeira equação em relação ao tempo. O resultado 
é 

de,. 

dt 
A seguir, usando as expressões para € e €p na Equação (1.55), encontramos que a equação 


acima se reduz a Je 
T = dep sen O + 06 (1.56) 


= T(Ê cos 0 cos À — ġ sen 0 sen d) + J(Ô cos 8 sen ġ + P sen 0 cos À) — kôsen6 
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As outras duas derivadas são obtidas por procedimento semelhante. Os resultados são 


t = —ĝëE, + $ cos ës (1.57) 
“e = —QsendE, — ġ cos dE (1.58) 


As etapas são deixadas como um exercício. Voltando agora ao problema de encontrar v, 
inserimos a expressão para de, /dt dada pela Equação (1.56) na Equação (1.54). O resultado 
final é 

7 = 6 + Erosen + Carb (1.59) 


dando o vetor velocidade em termos de suas componentes na tríade girada. 
Para encontrar a aceleração, diferenciamos a expressão acima em relação ao tempo. Ob- 


temos i 
eo dv si E ds d(ró sen 0) des d(r0) | råt 
ma e e dt dad "x 


Usando as fórmulas deduzidas anteriormente para as derivadas dos vetores unitários, encon- 
tramos prontamente que a expressão acima para a aceleração se reduz a 


> 


F Co 


L rósen6 


a = (F— r$? sen?0 — r0?)E, — (rö + 250 — ro? sen 0 cos 0)eg + 
+(ró sen 8 + 27 sen O + 2rĝġ cos 0)E, (1.60) 


dando o vetor aceleração em termos de suas componentes na tríade €., €o, €g- 


1.26 Velocidade Angular 


Uma partícula cujo vetor posição inicial é r sofre um deslocamento caracterizado por uma 
rotação de um ângulo ĝo em torno de um eixo de direção definida por um vetor unitário €, 
Figura 1.22. A partícula se moverá então ao longo de um arco de círculo de raio r sen 6 onde 
0 é o ângulo entre r e é O módulo do deslocamento da partícula é então |ór| = r sen 06%, 
e a direção do deslocamento é perpendicular a ambos r e é Podemos, pois expressar o 
deslocamento vetorialmente como um produto vetorial, ou seja 


ÔT = ÔE XT 
Dessa maneira, a velocidade da partícula é dada por 
r= lim = = ë xT (1.61) 
Podemos agora introduzir o vetor œ definido como o produto 
o = Gê 
chamado de velocidade angular. A velocidade da partícula é então expressa por 


P=0xF (1.62) 
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Vamos agora mostrar que velocidades angulares obedecem a regra de adição vetorial. 
Consideremos o deslocamento causado por uma rotação infinitesimal 94, em torno de um 
eixo €& seguido de uma segunda rotação ô> em torno de um eixo diferente &. A primeira 
rotação muda o vetor posição Y para um novo vetor posição T + d91€j x r. Dessa maneira o 
deslocamento líquido devido às duas rotações é 


Óri2 = 00161 XT + dQ26 X (T+ dQ6 X T) 


Figura 1.22: Deslocamento produzido por uma rotação. O raio da trajetória circular é r sen 0. 
Se as duas rotações angulares são ambas suficientemente pequenas de tal maneira que 
possamos desprezar o produto 901005, então encontramos depois da expansão que 
ÓT19 = (p€ + 0262) xT 


Se trocarmos a ordem das rotações, encontramos, de fato, o mesmo resultado, isto é, r12 = 
ôr21. Em outras palavras, as duas rotações infinitesimais são comutativas. Finalmente, vamos 
dividir por ôt e tomar o limite, como na Equação (1.61) acima. Podemos então escrever 


r= (O +o)xT 


para a velocidade da partícula onde wj = 9164 e üz = 9565. Então podemos descrever o 
movimento da partícula considerando uma única velocidade angular 


W = Wy+ Wo 


Dada pela regra normal de soma vetorial. 
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Exercícios de Treinamento 


1.1 


1.2 


1.3 


1.4 


Dados os dois vetores A = T + Je B= J— k. Encontre 


d) AxBe|AxB 

e) (å + 2B) . (2 — B) 

f) (A + B) x (A — B) 

Dados A =T +7 +K, B =T +2, Č=2-—K. Encontre 
a; A+B-G 

b) A . (B+ Č) e (+B). E 
)JAÃ.(BxOe(AxB). E 

d) A x (B x Č) e (A x B) x E 


Encontre o ângulo entre os vetores A=T +] + keB=T +7 (Dica: Estes dois vetores 
definem uma diagonal do cubo e uma diagonal da face do cubo). 


Dados os vetores que variam com o tempo À = Tcoswt+]senwteB=H+PJ+ BK. 
Encontre 


a) dA/dt e |dĀ/dt| 
b) 2B/dP e |d B/dt?| 
c) d(A . B)/dt e d(A x B)/dt 


Problemas 


1.5 


1.6 


1.7 


1.8 
1.9 


Para quais valores de q os dois vetores A=T+ J+ kq e B = Tq — 27 + 2kq são 
perpendiculares entre si. 


Prove a identidade vetorial A x (B x Č) = (A . OB-(A . B)G. 


Dois vetores A e B representam lados concorrentes de um paralelogramo. Prove que a 
área do paralelogramo é |A x B|. 


Prove a lei trigonométrica dos senos usando métodos vetoriais. 


Três vetores A, Bec representam lados concorrentes de um paralelepípedo. Mostre 
que o volume do paralelepípedo é |A . (B x ©)|. 


>> 


1.10 Expresse o vetor T +7 em termos da tríade T'J’k’ onde os eixos x'z7” são girados em 


torno do eixo y (que coincide com o eixo y’) de um ângulo de 60º. 
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1.11 Mostre que o módulo de um vetor não é alterado por uma rotação. Use a matriz 


cos0 sen 0 
— sen6 cos 0 
0 0 1 


para uma rotação de um ângulo 0 em torno do eixo z. 


1.12 Encontre a matriz transformação para uma rotação de um ângulo O em torno do eixo 
z seguida por uma rotação de um ângulo & em torno do eixo y’. 

1.13 Os dois grupos de vetores à, b, Ceg 7, b', € c' são chamados de recíprocos se à . = 
b.b=e.dt=le qualquer outro bronu escalar misturado como à . b' = 0. Mostre 
que: 


= (E x b)/Q, # = (b x ©)/Q, b’ = (€ x 8)/Q onde Q = 8 . (b x ©). 
1.14 Encontre um conjunto de vetores que sejam recíprocos ao conjunto T, J erT +7 + k. 


1.15 Uma partícula move-se numa trajetória elíptica dada pela equação r = Tb coswt + 
J2bsenwt. Encontre a velocidade escalar em função do tempo. 


1.16 No problema anterior encontre o ângulo entre o vetor velocidade e o vetor aceleração 
no instante t = 7/4w. 


1.17 A posição de uma partícula é dada em coordenadas polares planas por r = be*t, 0 = ct. 
Mostre que o ângulo entre o vetor velocidade e o vetor aceleração permanece constante 
a proporção que a partícula se move para fora numa trajetória em espiral. 


1.18 Uma partícula se move num círculo de raio b. Se a velocidade escalar da partícula 
varia com o tempo de acordo com a equação v = At? para que valor, ou valores de t o 
vetor aceleração forma um ângulo de 45° com o vetor velocidade? 


1.19 Uma partícula se move numa trajetória helicoidal tal que sua posição, em coordenadas 
cilíndricas, é dada por 
R=b 4 = wt z = e 


Encontre a velocidade escalar e o módulo da aceleração como funções do tempo. 


1.20 Mostre que o módulo da componente tangencial da aceleração é dada pela expressão 
Sr 
e que o módulo da componente normal é 


an = /a2 — a? 


1.21 Use o resultado acima para encontrar as componentes normal e tangencial da aceleração 
como função do tempo no Problema 1.19. 
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1.22 Prove que V . à = vi e como consequência que para uma partícula que se move v 
e à são perpendiculares entre sí se a velocidade escalar v for constante. (Sugestão: 
Diferencie ambos os lados da equação V . V = v? com relação a t. Lembre-se que ù não 
é o mesmo que |ãl.) 

1.23 Prove que 


zE- (vxal=r.(vxa) 


1.24 Prove que |Y x a) = v2/p, onde p é o raio de curvatura da trajetória de uma partícula 
em movimento. 


1.25 Usando o fato de que o vetor tangente unitário 7 pode ser expresso como 


F= 


e |< 


encontre uma expressão para o vetor normal unitário em termos de à, a, V, v, e ù. 


1.26 Um aro de raio b é colocado numa montagem sobre mancais e gira como se segue: O 
aro gira com velocidade angular constante wı em torno de seu próprio eixo, o qual, 
por sua vez roda com velocidade angular constante w em torno de um eixo vertical 
de tal maneira que o eixo do aro permanece num plano horizontal e o centro do aro 
fica imóvel. Use coordenadas esféricas para encontrar a aceleração de qualquer ponto 
na periferia do aro. Em particular, encontre a aceleração do ponto mais alto do aro. 
(Sugestão: Use o fato de que coordenadas esféricas podem ser escolhidas tais que r = b, 
0 = wt e d=uwst). 


Respostas de problemas selecionados de números ímpa- 
res 


1.1 Bisa 
(BITER 
(c) 

(d) =1+7 +k, 3 

1.3 arccos 4/2/3 = 85 15 

5 (-1+v17)/4 

1.15 wbv1 + 3cos%wt 

1.19 bw? + 4t, Vbwt + 4e 


HA 


121 pl do 4c2t 
; b2w2+4c2t2 ? vb2w2+4c2t? 


1.25 (av — vo)/|a x 7] 
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Capítulo 2 


Mecânica Newtoniana — Movimento 
Retilíneo de uma Partícula 


Como foi estabelecido na introdução, a dinâmica é a parte da mecânica que lida com as 
leis físicas que governam o movimento real dos corpos. Um dos objetivos fundamentais da 
dinâmica é predizer, dentre todos os movimentos possíveis de um sistema material qual movi- 
mento particular ocorreria numa dada situação. Nosso estudo da dinâmica, neste momento, 
será baseado nas leis do movimento como foram primeiramente formuladas por Newton. 
Num Capítulo posterior veremos maneiras alternativas de expressar as leis do movimento 
numa forma mais avançada utilizando as equações de Lagrange e Hamilton. Entretanto, 
estas equações não constituem teorias diferentes uma vez que podem ser derivadas das leis 
de Newton. 


2.1 As Leis de Newton do Movimento 


O leitor, sem dúvida, já está familiarizado com as leis de Newton do movimento. Elas são 
as seguintes: 


I - Todo corpo permanece em seu estado de repouso ou movimento retilíneo uniforme a 
não ser que seja obrigado, por uma força, a mudar tal estado. 


II - Mudança de movimento é proporcional à força aplicada e ocorre na direção da força. 


III - A cada ação corresponde sempre uma reação em sentido oposto, ou seja, as ações 
mútuas de dois corpos são sempre iguais, em módulo, e com sentidos opostos. 


Examinemos agora estas leis detalhadamente. 


2.2 Primeira Lei de Newton — Referenciais Inerciais 
A primeira lei descreve uma propriedade comum a toda matéria: a inércia. Ela afirma que 


um corpo em movimento descreve uma trajetória retilínea com velocidade constante a menos 
que alguma influência chamada força o impeça de fazê-lo. O fato de um corpo se mover ou 
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não em linha reta com velocidade constante depende não somente de influências externas 
(forças) mas também do sistema de referência utilizado para descrever o movimento. A 
primeira lei realmente nos leva à definição de um tipo particular de sistema de referência 
chamado sistema de referência Newtoniano ou inercial. Tal sistema é aquele no qual a 
primeira lei de Newton é válida. Sistemas de referência acelerados ou em rotação não são 
inerciais. 

A questão que ocorre naturalmente é como podemos determinar se um dado sistema de 
referência é ou não um sistema inercial. A resposta não é simples. A fim de eliminarmos 
todas as forças que poderiam atuar em um corpo é necessário isolá-lo completamente. Isto 
certamente é impossível uma vez que pelo menos forças gravitacionais estariam atuando 
sobre o corpo, a não ser que ele estivesse colocado a uma distância infinita de todos os 
outros corpos. 


Na prática, em casos onde a precisão desejada não seja muito alta, um sistema fixo à 
Terra é aproximadamente inercial. Por exemplo, uma bola de bilhar parece mover-se em 
linha reta com velocidade constante, exceto durante as colisões com outras bolas ou com as 
tabelas. Entretanto, se este movimento fosse medido com bastante precisão observaríamos 
que a trajetória é ligeiramente curvilínea. Isto porque a Terra está girando e um sistema 
fixo a ela não é, na verdade, inercial. Um sistema que usasse o centro da Terra, o centro do 
Sol e uma estrela distante como pontos de referência seria melhor do que o anterior. Mas 
nem mesmo este sistema seria estritamente inercial por causa do movimento orbital da Terra 
em torno do Sol. Prosseguindo, obteríamos uma aproximação ainda melhor se usassemos o 
centro do Sol e duas estrelas distantes como pontos de referência. Aceita-se que o sistema 
de referência por excelência, sob o ponto de vista da mecânica de Newton, seria aquele que 
levasse em conta a média de toda a matéria do universo. 


2.3 Massa e Força — Segunda e Terceira Leis de New- 
ton 


Todos nós estamos familiarizados com o fato de que, além de ser mais difícil levantarmos uma 
pedra grande do que, digamos, um pedaço de madeira, é necessário também um esforço maior 
para colocá-la em movimento (ou para pará-la). Dizemos que a pedra tem maior inércia do 
que a madeira. A medida quantitativa de inércia é chamada massa. Consideremos dois 
corpos A e B. Como medir a inércia de um corpo em relação ao outro? Podemos imaginar 
muitas experiências para responder a essa questão. Se dois corpos interagissem, por exemplo 
através de uma mola, observaríamos, por meio de medidas cuidadosas, que as acelerações 
dos dois corpos teriam sentidos opostos e que a razão entre seus módulos seria constante. 
(Na suposição de que as acelerações fossem medidas em relação a um sistema de referência 
inercial e que apenas a influência da interação dos dois corpos 4 e B estaria presente). 
Podemos expressar este fato importante e fundamental pela equação 


— > MBA (2.1) 
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A constante upa é, de fato, a medida da inércia de B em relação a A. Da Equação (2.1) 
obtemos que uga = l/uag. Podemos expressar uga como a razão 


mB 
HBA = —— 


MA 
e usar a inércia de um corpo padrão como unidade. A razão mp/ma não deve depender 
do corpo escolhido como padrão. Isto seria verdade se, considerando um terceiro corpo C, 


pudéssemos escrever 
HBC 


—— = HBA 
HAC 
Observamos que a razão acima se verifica. Chamamos a quantidade m de massa. 
Estritamente falando, m deve ser chamado de massa inercial, uma vez que sua definição 
é baseada nas propriedades de inércia. Na prática a razão entre as massas é usualmente 
determinada por pesagem. O peso ou força gravitacional é proporcional àquilo que chamamos 
de massa gravitacional de um corpo. Todas as experiências realizadas até hoje indicam que a 
massa inercial é proporcional à massa gravitacional. Por isso, considerando nossos objetivos, 
não precisaremos distinguir os dois tipos de massa. 
O fato fundamental contido na Equação (2.1) pode agora ser escrito na forma 
dv A dv'B 
MA TE = TMB- H 
O produto da massa pela aceleração na equação acima é a mudança de movimento referida 
na segunda lei de Newton e, de acordo com essa lei, é proporcional à força. Em outras 
palavras, podemos escrever a segunda lei como 
F = RA 
dt 


(2.2) 


(2.3) 


onde F é a força e k é uma constante de proporcionalidade. E usual tomar k = 1 e escrever! 


- dy 
F = m— 2.4 
my (2.4) 
A equação acima é equivalente a 
=- d(myv) 
F= 2.5 
H (2.5) 


se a massa fosse constante. De acordo com a teoria da relatividade, a massa de um corpo em 
movimento não é constante mas depende da velocidade do corpo. Então as Equações (2.4) e 
(2.5) não são equivalentes. Entretanto, para velocidades pequenas comparadas à velocidade 
da luz, 3 x 108 m/s, a variação de massa é desprezível. 

O fato fundamental expresso pela Equação (2.2) pode ser interpretado, de acordo com 
a Equação (2.4), como sendo a afirmativa de que dois corpos que interagem diretamente 
exercem forças iguais e opostas um no outro 


h= 


1No sistema MKS a unidade de força, definida pela Equação (2.4) é chamada newton. Então a força de 
1 newton ocasiona uma aceleração de 1 m/s? em um objeto de 1 kg de massa. No sistema CGS a unidade 
de força (1g x 1em/s”) é chamada dina. 
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Isto está incorporado no enunciado da terceira lei de Newton. As forças são chamadas ação 
e reação. 

Existem situações nas quais a terceira lei não é verificada. Quando dois corpos estiverem 
separados por uma grande distância e interagirem um com o outro através de um campo de 
força que se propaga com velocidade finita, tal como a interação entre cargas elétricas em 
movimento, as forças de ação e reação não serão sempre iguais e opostas. Entretanto, em 
tais casos, é possível considerar o campo de força como um terceiro corpo com sua própria 
ação e reação. Então a terceira lei não precisa ser descartada. Veja Seção 4.1 e as referências 
citadas lá. 

Uma grande vantagem do conceito de força é que ele nos permite restringir nossa atenção 
a um único corpo. O significado físico da idéia de força é que, para uma dada situação, usu- 
almente podemos encontrar uma função relativamente simples das coordenadas, chamada 
função força, que quando igualada ao produto da massa pela aceleração descreve correta- 
mente o movimento do corpo. Esta é a essência da mecânica Newtoniana. 


2.4 Momentum Linear 


O produto da massa pela velocidade é chamado momentum linear e é representado pelo 
símbolo p. Então 


PES MA (2.6) 
O enunciado matemático da segunda lei de Newton, Equação (2.5), pode ser escrito 


- dp 
F = — 2 
T (2.7) 


Ou seja, a força é igual à razão temporal da variação do momentum linear. 


A terceira lei, ação e reação, pode ser convenientemente expressa em termos do momen- 
tum linear. Se dois corpos A e B interagem, teremos 


Consequentemente 


Pa + pp = constante 


Então, de acordo com a terceira lei, o momentum linear total de dois corpos que interagem 
permanece sempre constante. 

A conservação do momentum linear de dois corpos que interagem é um caso especial de 
uma regra mais geral, a saber: o momentum linear de qualquer sistema isolado permanece 
constante no tempo. Discutiremos detalhadamente esta regra mais tarde. Ela é conhecida 
como a lei de conservação do momentum linear e é uma das leis mais básicas da física. 
Supõe-se que sua validade se estende àqueles casos nos quais a mecânica Newtoniana falha. 
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2.5 Movimento de uma Partícula 


A equação fundamental do movimento de uma partícula é dada pela segunda lei de Newton, 
Equação (2.4). Quando uma partícula estiver sob a influência de mais de uma força, observar- 
se-á experimentalmente que elas se somam vetorialmente, isto é, 


F=5 F=mo=mã (2.8) 


Conhecendo-se a aceleração da partícula, então a equação do movimento (2.8) dá a força 
que atua nela. Os problemas usuais da dinâmica da partícula, entretanto, são aqueles nos 
quais as forças são funções conhecidas das coordenadas incluindo o tempo, e o objetivo é 
encontrar a posição da partícula em função do tempo. Isto envolve a solução de um conjunto 
de equações diferenciais. Em alguns problemas é impossível obter soluções das equações 
diferenciais do movimento em termos de funções analíticas. Neste caso, devemos utilizar 
métodos de aproximação. Em muitas aplicações práticas, tais como balística, movimento 
de satélites, etc, as equações diferenciais são tão complicadas que é necessário recorrer à 
integração numérica, frequentemente utilizando computadores, para predizer o movimento. 


2.6 Movimento Retilíneo — Aceleração Constante 


Quando uma partícula móvel permanecer em uma única linha reta, o movimento será dito 
retilíneo. Neste caso, sem perda de generalidade poderemos colocar o eixo x na direção do 
movimento. À equação geral do movimento poderá ser escrita 


El rt) = më 


Consideremos alguns casos especiais nos quais a integração da equação do movimento possa 
ser feita por métodos elementares. 

A situação mais simples será aquela na qual a força for constante. Neste caso, a aceleração 
será constante 


dv 
— = — = constante = a 
dt m 
e a solução será facilmente obtida por integração direta em relação ao tempo. 


at + vo (2.9) 


v 


l o 
g = q + vot + xo (2.10) 


onde vo é a velocidade inicial e xo a posição inicial. Eliminando o tempo nas Equações (2.9) 
e (2.10), obtemos 
Za(x — xo) = V — vo (2.11) 


O leitor reconhecerá nas equações acima as equações do movimento uniformemente ace- 
lerado. Existem numerosas aplicações fundamentais. Por exemplo, a aceleração de um corpo 
em queda livre próximo à superfície da Terra, desprezando a resistência do ar, é aproxima- 
damente constante. Representaremos a aceleração de um corpo em queda livre por g. (Em 
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módulo, g = 9,8m/s”). A força da gravidade (peso), cujo sentido é para baixo, é conse- 
quentemente igual a mg. A força gravitacional está sempre presente independentemente do 
movimento do corpo, onde outras forças podem estar atuando. 


Exemplo 


Considere uma partícula que escorrega para baixo num plano liso inclinado de um ângulo 
0 em relação à horizontal como ilustrado na Figura 2.1(a). Colocamos o eixo x na direção 
paralela ao plano inclinado e o orientamos positivamente para baixo, como indicado. 


SE 
u me cosĝð 


(a) (b) 
Figura 2.1: Uma partícula deslizando para baixo em um plano inclinado (a) liso; (b) com 


atrito. 


A componente da força gravitacional na direção x vale mg sen 0. É uma constante e por isso 
o movimento é descrito pelas Equações (2.9), (2.10) e (2.11) onde 


a = — = gsen 0 
m 


Vamos supor que, ao invés de liso, o plano fosse rugoso, isto é, que exercesse uma força de 
atrito f na partícula. Então, como ilustrado na Figura 2.1(b), a força resultante na direção 
x seria igual a mg sen 0 — f. No caso de um corpo deslizando sobre outro observamos que o 
módulo da força de atrito é proporcional ao módulo da força normal N, isto é, 


f=uN 


onde a constante de proporcionalidade u é o coeficiente de atrito cinético. No exemplo em 
questão a força normal N, como pode ser visto na figura, é igual a mg cos 0, portanto 


f = umg cos 6 
Consequentemente, a força resultante na direção «x seria igual a 
mg sen 0 — umg cos 6 


Novamente a força seria constante e as Equações (2.9), (2.10) e (2.11) poderiam ser aplicadas 
com 


F 
a = — = g(sen 0 — u cos 0) 
m 
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A velocidade da partícula crescerá se a expressão entre parênteses for positiva, isto é se 
0 > arctg yu. O ângulo arctg u, usualmente representado por e, é chamado ângulo de atrito. 
Se 0 = e, então a = 0, e a partícula deslizará para baixo com velocidade constante. Se O < e, 
então a será negativo e a partícula eventualmente irá parar. Devemos notar que no caso de 
movimento para cima a força de atrito terá seu sentido invertido em relação ao caso anterior; 
isto é, terá o sentido de x positivo. A aceleração (realmente desaceleração) será então 


a = g(sen 0 + ucos0) 


2.7 O Conceito de Energias Cinética e Potencial 


Ocorre geralmente que a força que atua sobre uma partícula depende de sua posição em 
relação a outros corpos. Isto é verdade no caso de, por exemplo, forças eletrostáticas ou 
gravitacionais e tensão ou compressão elásticas. Se a força independesse da velocidade e do 
tempo, a equação diferencial para movimentos retilíneos seria simplesmente 


F(x) = më 


Usualmente é possível resolver esse tipo de equação diferencial por vários métodos. Uma 
maneira útil e significativa é escrever a aceleração na forma abaixo 
dt dzd dv 
== — VU 
dt  dtdax da 


t= 


e a equação diferencial do movimento pode ser escrita 


dv mdv) dT 
F(x) = M a a dE (2.12) 


A quantidade T = imo’ é chamada energia cinética da partícula. Podemos expressar a 


Equação (2.12) na forma integral 
1 
J Fedr = far = a má? + constante 


A integral f F(x)dx é o trabalho realizado sobre a partícula pela força F(x). Vamos definir 
a função V (x) tal que 
ita F(x) (2.13) 
dx 
A função V(x) é chamada energia potencial. Ela é definida a menos de uma constante 


arbitrária. Em termos de V(x) o trabalho é 


J Eedi = — | ——dx = —V (x) + constante 


consequentemente podemos escrever 


1 
T+V = amu? V(x) = constante = E (2.14) 
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chamamos a constante E de energia total. Dito de outra maneira: Para o movimento uni- 
dimensional, se a força exercida sobre a partícula for uma função da posição apenas, então 
a soma da energia potencial e cinética permanecerá constante durante todo o movimento. 
Neste caso, a força será dita conservativa. (No próximo Capítulo encontraremos uma dis- 
cussão mais completa sobre força conservativa). Forças não conservativas, ou seja, aquelas 
para as quais não existir uma função energia potencial, são usualmente de natureza dissipa- 
tiva, tal como o atrito. 


À V(x) 


região permitida 


pontos de 
retorno 


xy 


Figura 2.2: Gráfico da função energia potencial mostrando a região permitida e os pontos 
de retorno do movimento para um dado valor da energia total E. 


O movimento da partícula pode ser obtido explicitando v na equação de conservação de 
energia (2.14). 

dx 2 

= — = —|E -V 2.15 

= = + SIE - v(a) (2.15) 


que pode ser escrita na forma de uma integral 


J H (2.16) 
JZE- V(a) 
dando t em função de zx. 

Vemos, na Equação (2.15), que a velocidade só será real para valores de x tais que 
V(x) seja menor ou igual à energia total E. Fisicamente, isso significa que a partícula está 
confinada às regiões para as quais V(x) < E. Além disso, a velocidade se anula nos pontos 
onde V(x) = E. Isto é, a partícula deve parar e inverter seu movimento nestes pontos. Eles 
são chamados pontos de retorno do movimento. A Figura 2.2 ilustra essa discussão. 


Exemplo 
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O movimento de um corpo em queda livre discutido anteriormente como exemplo de 
força constante é um caso especial de movimento conservativo. Orientando o eixo x positivo 
para cima, a força gravitacional será —mg e a função energia potencial V = mgz + C. C 
é uma constante arbitrária cujo valor depende meramente da escolha do nível de referência 
para V. Para C = 0, a energia total será 


1 
E = ama? + Mgr 


Vamos supor, por exemplo, que um corpo seja lançado para cima com velocidade inicial vo. 
Fazendo x = 0 no ponto de lançamento 

1 1 

E = mv? = mit? + mge 

I Oe E 
O ponto de retorno é aquele no qual a altura atingida pelo corpo é máxima. Nós o obtemos 
fazendo t = 0 1 
zmo = Mit maz 


ou 3 
vo 


2g 


Integrando a equação de energia obtemos para este movimento 


h = Tmar — 


O leitor deve verificar que esta expressão se reduz à mesma relação entre x e t dada pela 
Equação (2.10) quando se coloca a = —g. 


2.8 Força em Função do Tempo — Conceito de Impul- 
so 


Se for possível expressar a força que atua sobre as partículas em termos de uma função 
explícita do tempo então a equação do movimento será 


dv 
F(t) = m— 
(t) T 
Esta equação pode ser integrada diretamente dando o momentum linear (e portanto a velo- 
cidade) em função do tempo 


/ F(tdt = mo(t) + O (2.17) 


onde C é uma constante de integração. A integral f F(t)dt, chamada impulso, é igual ao 
momentum linear comunicado à partícula pela força F(t). (O uso do conceito de impulso 
será retomado no Capítulo 4). 
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A posição da partícula em função do tempo pode ser obtida por uma segunda integração 


como se segue 
e= fuoa- [|| Ea (2.18) 


Devemos notar que a solução da equação poderá ser escrita como uma simples integral dupla 
apenas nos casos em que a força for uma função explícita do tempo. Para todos os outros 
casos devemos usar os vários métodos para se resolver equação diferencial de segunda ordem 
a fim de obtermos a posição x em termos do tempo t. 


Exemplo 


Um bloco está inicialmente em repouso numa superfície horizontal lisa. No instante t = 0 
aplica-se uma força horizontal uniformemente crescente: F = ct. Encontre a velocidade e o 
deslocamento em função do tempo. 

A equação diferencial do movimento é 


Então 


onde a posição inicial do corpo é na origem a = ui 


2.9 Força Dependente da Velocidade 


Frequentemente ocorre que a força que atua numa partícula é função da velocidade da 
partícula. É o caso de, por exemplo, resistência viscosa exercida em um corpo que se move 
através de um fluido. No caso de resistência exercida por um fluido, observamos que para 
baixas velocidades a força é aproximadamente proporcional à velocidade, enquanto que para 
altas velocidades a proporcionalidade ao quadrado é mais conveniente. Se não houver outras 
forças atuando, a equação do movimento poderá ser escrita 


= M— 
dt 
Integrando uma única vez obtemos t em função de v 


mdu 


t= = t(v) (2.19) 


Podemos omitir a constante de integração porque seu valor depende somente da escolha da 
origem do tempo. Supondo que v possa ser explicitado na equação acima, isto é, 


v = v(t) 
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então uma segunda integração nos dará a posição x em função de t 


J$ 1 v(t)dt = z(t) (2.20) 


Exemplo 


Vamos supor que lançamos um corpo com velocidade inicial vo num plano horizontal 
liso, mas que existe uma resistência do ar proporcional a v; isto é, F(v) = —cv, onde c é 
uma constante de proporcionalidade. (O eixo x está ao longo da direção do movimento). A 
equação diferencial do movimento é 


Integrando, teremos 


v mdv m v 
t= / — = —— ln ( ) 
vo Cu c Vo 


Multiplicando ambos os membros por — < e tomando o exponencial obteremos facilmente v 
em função de t. O resultado é 


ct 


U = Ve m 


A velocidade decresce exponencialmente com o tempo. Uma segunda integração nos dará 


t ct 
LT = | voedt 
0 


Vemos, da equação acima, que o corpo nunca ultrapassará a distância limite mvo/c. 
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2.10 Movimento Vertical num Meio Resistivo 
Velocidade Terminal 


Um objeto caindo verticalmente no ar ou através de qualquer fluido está sujeito à resistência 
viscosa. Se a resistência for proporcional à primeira potência de v (caso linear) poderemos 
expressar esta força como —cv independentemente do sinal de v, uma vez que a força de 
resistência é sempre oposta ao sentido do movimento. A constante de proporcionalidade c 
depende do tamanho e forma do objeto e da viscosidade do fluido. Consideremos o sentido 
para cima positivo. A equação diferencial do movimento será 


Considerando g constante, teremos o caso de força dependente da velocidade, e poderemos 
escrever 


i mdv mdv m mI + cv 
= E e = n 
F(v) Jw =mg — cv c mg + cv 
Facilmente poderemos explicitar v 
m m ét 
paani (= n vo) o (2.21) 
c c 


Depois de um tempo suficientemente grande (t > m/c) o termo exponencial será desprezí- 
vel, e a velocidade se aproximará do valor limite —=2. A velocidade limite de um corpo em 
queda é chamada de velocidade terminal; é aquela velocidade na qual a força de resistência 
seja exatamente igual e oposta ao peso do corpo tal que a força resultante seja nula. O 
módulo da velocidade terminal de uma gota de chuva, por exemplo, é de 3 a 6 metros por 
segundo, dependendo do seu tamanho. 

A Equação (2.21) expressa v em função de t, então uma segunda integração dará x em 


função de t: 
t t 2 ct 
£ — To - | v(t)dt = DE 4 [= dia me) (1-5) (2.22) 
0 


c c2 c 


Vamos designar a velocidade terminal 7# por v, e m/c por T (que chamaremos de tempo 
característico). A Equação (2.21) pode então ser escrita da forma mais significativa 


v = =u; + (vi + voje“? (2.23) 


Então, um objeto que cai a partir do repouso (vo = 0) atingirá a velocidade de (1 — emb) 
vezes a velocidade terminal em um intervalo de tempo 7, (1 — e“2)v, no intervalo 27, e assim 
por diante. Depois de um intervalo de 107 a velocidade será praticamente igual a velocidade 
terminal, a saber 0, 99995v,. 
Se a resistência viscosa for proporcional a v? (caso quadrático), a equação diferencial do 
movimento será, lembrando que consideramos o sentido para cima positivo, 
—mg + cv? = es 


dt 
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O sinal negativo para o termo de resistência se refere à subida (v positivo), e o sinal positivo 
à descida. O duplo sinal é necessário no caso da força que envolve potência par de v. Como 
no caso anterior, a equação diferencial do movimento pode ser integrada dando t em função 
de v: 


d 
t= / = = —rarctg— + to (subida) 
=mg — cu Vt 
d 
r= / 3 d = —rarctgh > + to (descida) 
onde 
m las 
— =7T (tempo característico) 
cg 
e 
RI vi (velocidade terminal) 
c 
Explicitando v, 
to—t 
v = utg ( o ) (subida) (2.24) 
E 
t— to 
v = —utgh (descida) (2.25) 
T 


Se o corpo for abandonado em repouso em t = 0 então tg; = 0. Da definição de tangente 


hiperbólica, 
t et/7T ne: eTHT 
VS SU eoe 


Novamente, depois de um intervalo de tempo de alguns tempos característicos, praticamente 
obtemos a velocidade terminal, por exemplo, para t = 57, a velocidade é 0, 99991v,. A Figura 
2.3 mostra os gráficos da velocidade em função do tempo de queda para leis de resistência dos 
tipos linear e quadrática. É interessante notar que em ambos os casos, o tempo característico 
é igual a v;/g. Por exemplo, se a velocidade terminal de um paraquedas é de 1,2 m/s, o tempo 
característico será de 

1,2m/s 


9,8m/s2 = 0,12s 


Podemos integrar as Equações (2.24) e (2.25) para obtermos expressões explícitas da 
posição em função do tempo. 


2.11 Variação da Gravidade com a Altura 


Na Seção precedente consideramos g constante. Realmente, a atração gravitacional da Terra 
sobre um corpo situado acima de sua superfície varia com o inverso do quadrado da distância. 
(Lei de Newton da gravitação). Então a força gravitacional num corpo de massa m é 


GMm 


r2 


F= 
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V velocidade terminal 
t 
(eb) ; ot 
O “. resistencia linear 
(2) 
g resistencia quadratica 
O 
(ab) 
> 
0 L L 
0 T 2T 3T 
tempo 


Figura 2.3: Gráfico da velocidade em função do tempo de queda para um corpo sujeito a 
resistência do ar dos tipos linear e quadrática. 


onde G é a constante universal da gravitação, M a massa da Terra, e r a distância do centro 
da Terra ao corpo. Pode-se ver por inspeção que este tipo de força é obtida de uma função 
energia potencial que varia com o inverso da distância 


GM 
Ve) = -27 
m 
onde 
= oV 
— Or 
Se desprezarmos a resistência do ar, a equação diferencial para movimento vertical é 
z GMm 
mi (2.26) 
Escrevendo 7 = rÈ podemos integrar em relação a r e obter 
1 GM 
“mes e p (2.27) 
T 


onde E é a constante de integração. Esta é, de fato, a equação da conservação da energia: 
a soma da energia cinética e potencial permanece constante durante o movimento de um 
corpo em queda. 

Vamos aplicar a equação de energia ao caso de um projétil lançado para cima com 
velocidade vo a partir da superfície da Terra. A constante E é dada pelas condições iniciais, 


1 > GMm 
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onde rr é o raio da Terra. A velocidade em função da altura x é obtida combinando-se as 
Equações (2.27) e (2.28). O resultado é 


1 1 

v? = va + 2GM ( — ) (2.29) 
rr+H£ TT 

onde r = rr + «x. Da Equação (2.26), temos que a aceleração da gravidade na superfície da 


Terra 
= GM 


g 2 
tp 


A fórmula para a velocidade pode então ser escrita como 
AAI 
v? = vå — 2gx (1 + =) (2.30) 
TT 


A equação acima se reduz à fórmula familiar para um campo gravitacional uniforme 


v = vê —2gr 
quando x for muito menor do que rr tal que o termo = possa ser desprezado em comparação 
com a unidade. 

O ponto de retorno do movimento do projétil, isto é, a altura máxima atingida, será 
obtido colocando v = 0 e explicitando x. O resultado será 


vê Ny 
Ta = h = L ( — — (2.31) 
2g 2grr 
Novamente obteremos a fórmula usual 
pai 
2g 


se desprezarmos o segundo termo. 

Finalmente, vamos aplicar a fórmula exata (2.31) para obtermos a velocidade de escape, 
isto é, o valor de vo para o qual h seja infinito. É claro que isto ocorre quando a quantidade 
entre parênteses for igual a zero. O resultado é 


Ve = 4\/29rr 


o que dá 
ve % 11 km/s = 11 x 10° m/s 


para a velocidade de escape da superfície da Terra. 

Na atmosfera terrestre, a velocidade média das moléculas de ar (Os e N5) é de aproxima- 
damente 0,5 km/s, que é consideravelmente menor que a velocidade de escape, então a Terra 
retém sua atmosfera. (De acordo com a teoria cinética dos gases, a velocidade média das 
moléculas de um gás é igual a VE onde k é a constante de Boltzmann (1,38x 101º erg/K) e 
T é a temperatura absoluta). A Lua, por outro lado, não tem atmosfera porque a velocidade 
de escape na sua superfície, devido à pequena massa, é consideravelmente menor do que na 
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superfície da Terra; qualquer molécula de oxigênio ou nitrogênio eventualmente escaparia. 
A atmosfera da Terra, entretanto, não contém uma quantidade significativa de hidrogênio, 
embora ele seja o elemento mais abundante no universo como um todo. A atmosfera de 
hidrogênio teria escapado da superfície da Terra a muito tempo porque sua velocidade mo- 
lecular é suficientemente grande (devido à pequena massa da molécula de hidrogênio) tal 
que um número significativo de moléculas teria, em cada instante, velocidade maior do que 
a velocidade de escape. 


2.12 Força Restauradora Linear — Movimento Har- 
mônico 


Um dos casos mais importantes de movimento retilíneo, tanto sob o ponto de vista prático 
quanto teórico, é aquele provocado por uma força linear restauradora. O módulo desta força 
é proporcional ao deslocamento da partícula em relação à sua posição de equilibrio e seu 
sentido se opõe ao deslocamento. Tal força é exercida por uma corda elástica ou por uma 
mola que obedeça à lei de Hooke 


F = —k(X — a) = —kz (2.32) 


onde X é o comprimento total da mola e a o comprimento da mola não deformada. (Sem 
carga). A variável r = X —a é o deslocamento da mola em relação à sua posição de equilíbrio. 
A constante de proporcionalidade k é chamada de constante elástica. Consideremos uma 
partícula de massa m presa a uma mola como mostra a Figura (2.4a); a força que atua na 
partícula é dada pela Equação (2.32). Consideremos a mesma mola dependurada vertical- 
mente e tendo a mesma partícula presa no seu extremo. A força total que atua na partícula 
é 

F = —k(X — a) + mg (2.33) 


onde o sentido positivo é para baixo. Agora, neste último caso, medimos x relativamente 


à nova posição de equilíbrio, isto é, tomamos z = X — a — 8. Assim temos novamente 
F = —kx, e a equação do movimento, em ambos os casos, é 
më + kz = 0 (2.34) 


Esta equação diferencial do movimento é encontrada em uma grande variedade de pro- 
blemas físicos. No exemplo particular que estamos usando aqui, as constantes m e k se 
referem à massa do corpo e à constante elástica da mola, respectivamente, e o deslocamento 
x é um comprimento. A mesma equação é encontrada, como veremos mais tarde, no caso 
de um pêndulo, onde o deslocamento é um ângulo e as constantes envolvem a aceleração da 
gravidade e o comprimento do pêndulo. Novamente, esta equação se aplica a certos tipos de 
circuitos elétricos, onde as constantes representam os parâmetros do circuito e a quantidade 
x representa a corrente elétrica ou a voltagem. 

A Equação (2.34) pode ser resolvida de várias maneiras. É um exemplo de uma impor- 
tante classe de equações diferenciais conhecidas como equação diferencial linear com coefici- 
entes constantes. A equação geral de ordem n deste tipo é 

dr oč, o” 
Cor Pq Ca 
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-o 
S 
X an 3 
k CX 
a x ES 
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ATATATATATATAT o Q 
700000000 | = 
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equilibrio 
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Figura 2.4: Ilustração do oscilador harmônico linear por meio de um bloco de massa m e 
uma mola. (a) Movimento horizontal. (b) Movimento vertical. 


A equação é chamada homogênea se b = 0. Muitas, se não a maioria, das equações diferen- 
ciais da física são equações diferenciais lineares de 22 ordem. A fim de resolver a Equação 
(2.34) utilizaremos o método das tentativas experimentando como solução a função Ae% 
onde q é uma constante a ser determinada. Se x = Ae” for, de fato, uma solução, então 
para todos os valores de t teremos 


d 
m— (Ae?) + k(Ae*) = 0 
Gal Ae) + (Act) 
o que, depois de cancelarmos os fatores comuns, se reduz a equação, (esta equação é chamada 
equação auxiliar) 
my +k=0 


| k 
q = EMP = iwo 
m 


onde i = y —1, e wọ = VE . As soluções das equações diferenciais lineares são aditivas. (Isto 
é, se fi e fo forem soluções, então a soma fı + fo também será uma solução). A solução da 
Equação (2.34) é então 


ou seja 


r= Ae, A tuo (2.35) 
Desde que e"! = cosu + isen u, 
x = a sen wot + b cos wot (2.36) 
ou 
x = Acos(wot + 06) (2.37) 


são formas alternativas da solução. As constantes de integração nas soluções acima são 
determinadas a partir das condições iniciais. Por substituição direta podemos verificar que 
todas as três expressões são soluções da Equação (2.34). O movimento é uma oscilação 
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Figura 2.5: Gráfico do deslocamento em função do tempo para o oscilador harmônico. 


senoidal do deslocamento x. Por esta razão a Equação (2.34) é frequente chamada de equação 
diferencial do oscilador harmônico ou oscilador linear. 

O coeficiente wo é chamado frequência angular. O valor máximo de x é chamado ampli- 
tude da oscilação; é a constante A na Equação (2.37), ou Va? + b? na Equação (2.36). O 
período T da oscilação é o tempo requerido para um ciclo completo, como mostra a Figura 
2.5, isto é, o período é o tempo para o qual o produto wot aumenta de exatamente 27, então 


27 m 
T) = — = 2m, | — 2.38 
e Wo É k 


A frequência linear de oscilação fo é, por definição, o número de ciclos por unidade de tempo, 
portanto 


Wo = 27 fo 
1 1 jk 
Es do, Cera 2.39 
Jo To 2ry/m l ) 


E comum empregar-se a palavra frequência para designar tanto a frequência angular quanto 
a linear; fica claro do contexto a qual delas estamos nos referindo. 


Exemplo 


Uma mola leve se distende de b quando equilibra um bloco de massa m. Encontre o 
movimento resultante em função de t sabendo-se que no instante t = 0 o bloco foi abandonado 
em repouso a uma distância | abaixo de sua posição de equilíbrio. Primeiro, a fim de obtermos 
a constante da mola, notamos que na condição de equilíbrio estático 


F = —kb+ mg = 0 
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tal que 


Portanto a frequência angular de oscilação é 


|k TE 
üg = a o 
m b 
Para obtermos as constantes da equação de movimento 


xz = Acos(wot + 00) 


vamos fazer 
no instante t = 0. Mas 
e então 


logo 


é a expressão desejada. 


2.13 Considerações de Energia no Movimento Harmô- 
nico 
Considere uma partícula movendo-se sob a ação de uma força restauradora linear 


F = —kx. Vamos calcular o trabalho W realizado por uma força externa F, para levar 
a partícula da posição de equilíbrio (x = 0) a uma posição genérica x. Temos 


F=-F=kg 


x k 
W = f Fada Ei (kx)dz = T 
0 
O trabalho W fica armazenado na mola na forma de energia potencial 


k 
V(x)=W = a (2.40) 
A definição de V, Equação (2.13), requer que F = —% = kz. A energia total será então 
dada pela soma das energias cinética e potencial 


1 1 
E= ama? + ska? (2.41) 
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Podemos escrever a velocidade em função do deslocamento 


2E k 


t=/ ->r 
m m 


e integrar a fim de obtermos t em função de x, como segue 


x 
— C 
arccos ( 7) + 


t= [ da a 
= DE ky Vk 


2E 
Aco 
k 


e C é a constante de integração. Depois de integrar obtemos a mesma relação encontrada na 
Seção precedente, exceto que agora temos um valor explícito para a amplitude 4. Poderíamos 
ter obtido a amplitude diretamente a partir da equação da energia (2.41) lembrando que x 
deve situar-se entre 4/2E/k e —,/2E/k a fim de que % seja real. Isto está ilustrado na 
Figura 2.6 que mostra a função energia potencial e os pontos de retorno do movimento para 
diferentes valores da energia total E. 

Da equação da energia vemos que o valor máximo de x, que chamaremos Var, Ocorre 
em x = 0, e então teremos 


onde 


ou 


Figura 2.6: Gráfico da função energia potencial do oscilador harmônico. Mostra-se os pontos 
de retorno definindo a amplitude para dois valores de energia total. 
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2.14 Movimento Harmônico Amortecido 


A análise do oscilador harmônico feita acima é idealizada uma vez que não levamos em conta 
as forças de atrito. Pequenas ou não elas estão sempre presentes nos sistemas mecânicos. 
Analogamente, existe sempre resistência nos circuitos elétricos. Vamos considerar, por exem- 
plo, o movimento de um objeto dependurado em uma mola de constante k. Vamos supor que 
exista uma força viscosa retardadora variando linearmente com a velocidade (como na Seção 
2.8), ou seja, tal como aquela produzida pela resistência do ar. As forças estão indicadas na 
Figura 2.7. 

Sendo x o deslocamento em relação à posição de equilíbrio, então a força restauradora 
exercida pela mola será —kgx, e a força retardadora —cx, onde c é uma constante de propor- 
cionalidade. A equação diferencial do movimento F = më será portanto —ka — ct = më 
ou, rearranjando os termos 

mz+ct+kz=0 (2.42) 


posição de 
LA equilibrio 


Figura 2.7: O oscilador harmônico amortecido. 


Novamente tentaremos uma solução do tipo exponencial Ae”, que deverá satisfazer a 


equação 
d? d 
qt qt AM — 
mygale ) + caplÃe )+ k(Ae?”) = 0 


para qualquer t. Isto ocorrerá se q satisfizer a equação auxiliar mq? + cq + k = 0. As raízes 
são dadas pela fórmula quadrática 


Dees 
PE cve — Amk (2.43) 


2m 


Existem três casos fisicamente distintos 


i. c2> 4mk super-amortecimento 
ii. c2=4mk amortecimento crítico 
iii. œ <4mk sub-amortecido 
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i) No primeiro caso, vamos chamar —y e —% os dois valores reais de q dado pela Equação 
(2.43). A solução geral poderá ser escrita 


E Ape Mt + Age PE (2.44) 


Vemos que o movimento é não oscilatório, o deslocamento x decai a zero de uma maneira 
exponencial a medida que o tempo passa, como mostra a Figura 2.8. 


X amortecimentos 


super-amortecido 


critico 


Figura 2.8: Gráficos do deslocamento em função do tempo para os casos super-amortecido 
e amortecimento crítico do oscilador harmônico. 


ii) No caso de amortecimento crítico as duas raízes são iguais, e a Equação (2.44) não 
representará a solução geral desde que existe apenas uma função e” e uma única constante 
(A, + A2), onde y = c/2m. Para obtermos a solução geral neste caso voltaremos à equação 
diferencial original (2.42). Para raízes iguais esta equação pode ser fatorada 


ER 


Fazendo a substituição u = yx + de teremos 


d 
d 
(5+) e=0 


o que é facilmente integrável dando u = Aje 7. Portanto, da definição de u, yx + dx/dt = 
Ae7% que também pode ser escrito 


dx d 
A = — ] e = — (re 
1 (+ E ASG ) 


Uma segunda integração em relação a t dá Aıt = xe% — As, ou finalmente 


X— e VA + As) (2.45) 
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Esta equação representa um movimento não oscilatório, o deslocamento x decaindo a zero 
assintoticamente com o tempo, Figura 2.8. O amortecimento crítico produz um ótimo retorno 
à posição de equilíbrio para aplicações práticas tais como suspensões de galvanômetros, etc. 
iii) Se a constante de resistência c for suficientemente pequena tal que c? < 4mk, teremos 
o terceiro caso: sub-amortecimento. Neste caso q será complexo. As duas raízes da equação 
auxiliar serão números conjugados complexos e o movimento será dado pela solução geral 


r = Ape rHwts A riu (2.46) 

onde y = 5, € 
w = Pal a we — 7? (2.47) 
x = e7" (asen wit + bcos wt) (2.48) 


onde a = i(A, — A_) e b = (A, + A_). Podemos também escrever a solução como 
x = Ac cos(wit + 0o) (2.49) 


onde 


A=val+b e o= -arctg (5) 

A forma real da solução mostra que o movimento é oscilatório, e que a amplitude Ae 
decai exponencialmente com o tempo. Além disso, notamos que a frequência angular de 
oscilação wı é menor que a do oscilador não amortecido wọ. A frequência wı é chamada 
frequência natural. 

No caso de amortecimento fraco, ou seja se y for muito pequeno comparado a wo, podemos 
usar a relação aproximada 

2 
Wi E wo — Ma (2.50) 
2wo 
que é obtida expandindo-se em série o lado direito da Equação (2.47) pelo uso do teorema 
binomial e retendo somente os dois primeiros termos. 

A Figura 2.9 mostra o gráfico do movimento. Da Equação (2.49) segue que as duas 
curvas x = Ae" e x = — Ae” formam um envoltório da curva do movimento uma vez que 
o cosseno toma valores entre +1 e —1 nos pontos onde a curva toca o envoltório. Os pontos 
de contato estão separados por um intervalo de tempo de meio período, ou 7/w1, mas estes 
pontos não são os máximos e os mínimos do deslocamento x. Fica para o leitor encontrar os 
valores de t nos quais o deslocamento assume valores extremos. 


Considerações de Energia 
A energia do oscilador harmônico amortecido é, em qualquer instante, igual à soma da 


energia cinética imi? e da energia potencial tkz’: 


1 1 
E= zri + ska? 
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Figura 2.9: Gráfico do deslocamento em função do tempo para o oscilador harmônico sub- 
amortecido. 


No oscilador não amortecido E é constante. Vamos diferenciar a equação acima em relação 
a t a fim de obtermos a taxa temporal de variação de E. Temos então 


a = mit + kizx = (më + kr) 


Mas, da equação diferencial do movimento, Equação (2.42) 
mz + kz = —cã 


Consequentemente, 


dE : 
— = —cã 2.51 
E ci (2.51) 


É uma quantidade sempre negativa e representa a razão na qual a energia está sendo dissipada 
na forma de calor pelo atrito. 
Exemplos 


1. Uma partícula de massa m está presa a uma mola de constante elástica k. O amorte- 
cimento é tal que y = wo/4. Encontre a frequência natural. Da Equação (2.47) temos 


TE ao. 
0 i6 V16 YmV16 


2.15. MOVIMENTO HARMÔNICO FORÇADO — RESSONÂNCIA 63 


2. No problema acima, encontre a razão entre as amplitudes de duas oscilações sucessivas. 
Esta razão, da teoria anterior, é dada por 


Age = en 
A 
onde 1 
Tı EE E a ET 
h wi 
Portanto, em nosso problema 
2m [16 27 [16 
Ti == E 
wo V 15 4yV\15 
ou 
m /16 
ea 
E 
Portanto a razão entre as amplitudes de duas oscilações sucessivas é 
ge — (0, 20. 
2.15 Movimento Harmônico Forçado — Ressonância 


Nesta seção estudaremos o movimento do oscilador harmônico amortecido forçado por uma 
força externa harmônica, isto é, uma força que varia senoidalmente com o tempo. Vamos 
supor que esta força tenha frequência angular w e amplitude Fọ e portanto podemos escrever 


Feu = Focos(wt + 0) 
É conveniente, entretanto, usar a forma exponencial 
Fert = FoelwitO) 


ao invés da trigonométrica, embora ambas possam ser usadas. (A forma exponencial é 
equivalente a Fert = Fo cos(wt + 0) + iFosen(wt + 0). A equação diferencial resultante será 
satisfeita se as partes real e imaginária nos dois lados da equação forem iguais). A força total, 
então, será a soma de três forças: a força elástica restauradora — kz, a força de amortecimento 
viscoso —cțt e a força externa Fest. A equação diferencial do movimento será portanto 


—kr — ct + Fert = MË 


ou 
më + ci + kr = Fen = et) (2.52) 


A solução da equação diferencial linear acima é dada pela soma de duas partes, a primeira 
parte sendo a solução da equação diferencial homogênea mz + ct + kr = 0, que já foi 
resolvida na Seção precedente; e a segunda parte sendo qualquer solução particular. Como 
vimos, a solução da equação homogênea representa uma oscilação que eventualmente decai 


64 CAPÍTULO 2. MECÂNICA NEWTONIANA 


a zero — é chamado o termo transiente. Estamos interessados na solução que depende da 
natureza da força aplicada. Uma vez que esta força varia senoidalmente com o tempo e tem 
amplitude constante, é razoável esperar uma solução para a qual o deslocamento x também 
tenha uma dependência temporal senoidal. Portanto, para o regime estacionário, tentaremos 


uma solução da forma 
x = Aet") 


Se esta função tentativa for correta teremos 


d? ; 1 d a / a 1 ; 
i(wt+0') i(wt+O”) i(wt+O) i(wt+0) 
Mmg (Ae ) + n (Ae ) + kAe = He 


para qualquer valor de t. Depois de efetuar as operações indicadas e cancelar os fatores 
comuns teremos: 


—- mw? A + iwcA + kA = Mei) = Fylcos(8 — 6º) + isen(0 — 69] 
Identificando as partes real e imaginária, teremos 
A(k — mw?) = Fo cos 4 (2.53) 


cwA = Fosen ġ (2.54) 


onde chamamos de q a diferença de fase ou ângulo de fase 0 — 0’. Dividindo a segunda 
equação pela primeira e usando a identidade sen &/ cos q = tg À, obtemos 


CW 


tgo = (2.55) 


k — mw? 


Elevando-se ao quadrado ambos os lados das Equações (2.53) e (2.54), somando e lembrando 
a identidade sen?y + cos? = 1, obtemos 


A’ (k -mv P + ew A = FE 
Explicitando a amplitude das oscilações no regime estacionário A, teremos 


Fo 


A= (2.56) 
V(k — mw? + ew? 
e como, wo = 14/k/m e y = c/2m, podemos escrever 
2yw 
ted = —————— 2.57 
go (2 — 03) (2.57) 
A= Rom (2.58) 


MO — w2)2 + 4y2w? 


A equação acima relacionando a amplitude A à frequência da força externa é de im- 
portância fundamental. O gráfico, Figura 2.10, mostra que 4 assume um valor máximo 
numa certa frequência wp, chamada de frequência de ressonância. 
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(1) c=1/2 mo, 
(2)c=1/4 mo, 


amplitude 


(1) 


0 ll l 
20, 


Q) 
frequencia 


Figura 2.10: Gráfico da amplitude em função da frequência da força externa. 


A fim de obtermos a frequência de ressonância calculamos dA/dw da Equação (2.58) 
e igualamos o resultado a zero. Explicitando w da equação resultante encontramos que a 


frequência será dada por 
wW = wr = yw — 27? (2.59) 


No caso de amortecimento fraco, isto é, quando a constante de amortecimento c for muito 
pequena, c < 2//mk, ou, equivalentemente, se y < wo, então a frequência de ressonância 
w, será aproximadamente igual à frequência do oscilador harmônico simples (sem amorte- 
cimento) wo. Se utilizarmos o teorema binomial para expandir o lado direito da Equação 
(2.59) retendo apenas os dois primeiros termos, teremos 
2 
Wr © Wo — 5 (2.60) 
Wo 
As Equações (2.59) e (2.60) devem ser comparadas às Equações (2.47) e (2.50) que dão a 
frequência de oscilação do oscilador harmônico amortecido. Façamos € representar a quan- 
tidade 42 /wo. Então podemos escrever 


1 
w1 MW — 5€ (2.61) 
para o valor aproximado da frequência natural e 
Wr Fº wo — € (2.62) 


para o valor aproximado da frequência de ressonância. 
A amplitude no regime estacionário e na frequência de ressonância, que chamaremos 
Amaz, é obtida das Equações (2.58) e (2.59). O resultado é 


Amaz = Torm = f 


-n fT 
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No caso de amortecimento fraco, desprezamos 7? e escrevemos 


Fi Fi 
Anar z= E = E 
2yMwo  Cwo 


Então a amplitude da oscilação induzida na condição de ressonância se torna muito grande 
se a constante de amortecimento for muito pequena, e vice-versa. Em sistemas mecânicos 
pode ser que seja ou não desejável grandes amplitudes de ressonância. No caso de motores 
elétricos, por exemplo, usa-se borracha ou suporte de molas para minimizar a transmissão 
de vibração. A constante elástica de tal suporte é escolhida de maneira a garantir que a 
frequência de ressonância resultante seja bem diferente da frequência de operação do motor. 

A largura do pico de ressonância é frequentemente de interesse. Consideremos o caso de 
amortecimento fraco y < wo. Então na expressão para a amplitude no regime estacionário, 
Equação (2.58), poderemos fazer as seguintes substituições 


ww? = (wo+w)(wo -— w) 


Q 


2wo(wo — w) 
Vo X Yo 


Estas, juntamente com a expressão para Amaz, Nos permitirá escrever a equação da amplitude 


na forma 
AmarY 


Vu — w)? +7? 


A equação acima mostra que quando |wo — w| = y, ou equivalentemente, se 


A= 


(2.63) 


w=w ty 
então 
A? = Ea 
2 max 


Sabendo-se que a energia é proporcional a 42, vemos que 2% é a diferença de frequência entre 
os pontos para os quais a energia do oscilador é igual a metade da energia na frequência de 
ressonância. Então, y é uma medida da largura da curva de ressonância. Isto está ilustrado 
na Figura 2.10. 

Podemos designar a largura do pico de ressonância de uma outra maneira utilizando-se o 
parâmetro Q, chamado de fator de qualidade, de um sistema ressonante. Q é definido como 
Wr 


g= (2.64) 


ou, para amortecimento fraco 


Me (2.65) 


Então a largura a meia altura da curva energia x frequência é aproximadamente 


Wo 
Aw = 2y x — 
Q 
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ou, como w = 27 f 


dando a largura relativa do pico de ressonância. 


Osciladores de cristal de quartzo são usados para controlar a frequência de estações 
radiodifusoras. O Q do cristal de quartzo em tal aplicação é da ordem de 104. Valores 
tão grandes quanto este garantem que a frequência de oscilação permaneça precisamente na 
frequência de ressonância. 


T/2 = 


(1) c = 1/2 mo, 
(2) c = 1⁄4 mo, 


0 (0) 20 
0 É: 0 
frequencia 


Figura 2.11: Gráfico do ângulo de fase em função da frequência da força externa. 


A diferença de fase à entre a força aplicada e a resposta do sistema é dada pela Equação 
(2.57). Esta equação está representada na Figura 2.11 que mostra à em função de w. Vemos 
que a diferença de fase é pequena para w pequeno, ou seja, a resposta do sistema está em 
fase com a força aplicada a ele. Na frequência de ressonância à aumenta para 7/2, isto é, a 
resposta está defasada de 90º em relação à força aplicada. Finalmente, para valores muito 
grandes de w, o valor de q se aproxima de 7, portanto o movimento do sistema está com 
180º de defasagem em relação à força externa. 
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Análogos Elétrico-Mecânico 


Existe uma analogia precisa entre o movimento de um sistema mecânico de massas e molas 
com forças de atrito do tipo estudado previamente e um circuito constituido de elementos 
indutivo, capacitivo e resistivo, no qual flui uma corrente elétrica. Assim, se a corrente 

= dq/dt (q sendo a carga) flui através de um indutor L, a diferença de potencial nos 
extremos do indutor é Lã e a energia armazenada é ELR’. Portanto indutância e carga são 
os análogos à massa e deslocamento, respectivamente, e diferença de potencial é o análogo à 
força. Similarmente, se um capacitor de capacitância C tiver carga q, a diferença de potencial 
será & e a energia armazenada é lg, Consequentemente, vemos que o recíproco de C é o 
análogo à constante elástica, k, da mola. Finalmente se uma corrente elétrica i flui através 
de um resistor de resistência R, a diferença de potencial entre seus extremos é iR = qR, 
e a razão temporal da energia dissipada é i? R = qºR em analogia com cù? para o sistema 


mecânico, Equação (2.51). A Tabela 2.1 resume a situação 


Tabela 2.1: 
MECÂNICO ELÉTRICO 
x deslocamento q carga 
x velocidade q=1 corrente 
m massa L indutância 
k constante elástica C! recíproco da capacitância 
c constante de amortecimento R resistência 
F força V diferença de potencial 
Exemplos 


1. Determine a frequência de ressonância e o fator de qualidade para o oscilador amor- 
tecido do Exemplo 1 da página 62. 


Temos 
Wr = ywi -— 29 
2 
Wo 


| 
€ 
(<) 

ol N 
II 
sim] 
00] Í 


(2.66) 


para a frequência angular de ressonância. O fator de qualidade é dado por 
Wr a 
Q= === =2W/=-=1,87 
27 Uwo E 4) 


2. Se a frequência da força aplicada no oscilador acima é w9/2, encontre o ângulo de fase 
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à. Da Equação (2.57), temos 


2(wo/4)(wo/2) o 1/4 ed 


EO” GT 


Portanto 


1 
9 = arctg (5) = 184 


2.16 Movimento sob a ação de uma Força Periódica 
não Senoidal 


A fim de determinarmos o movimento de um oscilador harmônico sujeito à uma força ex- 
terna periódica mas não senoidal é necessário utilizar um método mais complicado que o 
utilizado na seção precedente. Neste caso mais geral é conveniente usar o princípio da super- 
posição. Este princípio estabelece que se a força aplicada F(t) que atua no oscilador puder 
ser expandida numa soma 
Fl =F Fut) 
n 


tal que as equações diferenciais 
min + cin + kn = Falt) 
sejam individualmente satisfeitas pelas funções 
tt) 
então a equação diferencial 


më + ci + kz = F(t) = X. lt) 


n 


será satisfeita pela função 


y= X Dali) 


Da linearidade das equações diferenciais do movimento segue imediatamente que o teorema 
acima é valido. 

Em particular, quando a força externa F(t) for periódica de frequência angular w, po- 
deremos desenvolvê-la em série de Fourier. (Consulte um livro texto sobre o método de 
Fourier). De acordo com a teoria de série de Fourier, podemos expressar F(t) numa soma 
de termos em senos e cossenos, ou alternativamentepodemos escreve-la em termos de uma 
soma de exponenciais complexas, isto é, 


F(t) = X Fp" (n = 0, +1, 2,...) 


Os coeficientes são dados por 
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Como na seção precedente, o movimento real será dado pela soma de duas partes, a 
saber, um termo transiente, que desprezaremos, e a solução estacionária. 


x(t) = Ag + Aje + Age +... 


O primeiro termo Ao é uma constante cujo valor depende da forma de F(t). Ao é nulo 
para força simétrica. O segundo termo é a resposta de um oscilador forçado na frequência 
fundamental w. O terceiro é a resposta no segundo harmônico da força aplicada e assim por 
diante. 

Podemos usar a teoria da seção precedente para obtermos as amplitudes A„ em termos 
dos coeficientes Fp. Então, da Equação (2.58) temos 


Falm 


ue 
Vu — n2w?)2 + 4y2n2y? 


(2.67) 


Da análise acima, vemos que o movimento estacionário final é periódico e que ao harmônico 
nw que for mais próximo da frequência de ressonância w, corresponderá a maior de todas as 
amplitudes. Em particular, se a constante de amortecimento y for muito pequena e se algum 
dos harmônicos da força coincidir com a frequência de ressonância de tal maneira que, para 
algum valor de n, tenhamos 

Wr = NW 


então a amplitude A, correspondente a este harmônico será dominante. Consequentemente 
o movimento resultante de um oscilador será aproximadamente senoidal mesmo que a força 
externa aplicada não o seja. 


Exercícios de Treinamento 


2.1 Uma partícula de massa m está inicialmente em repouso. Subitamente, no instante 
t = 0, aplica-se uma força constante Ho. Depois de um intervalo de tempo to a força 
subitamente dobra de valor e permanece constante daí para a frente. Encontre a 
velocidade da partícula e o deslocamento total no instante t = 2to. 


2.2 Encontre a velocidade v e a posição x em função de t de uma partícula de massa m 
que parte do repouso em t = 0 sujeita às seguintes forças: 


a) F = Fo 
b) F= F +bt 
c) F = Fo cos wt 
d) F = kt? 


2.3 Encontre a velocidade v em função do deslocamento x de uma partícula de massa m 
que parte do repouso em x = 0 sujeita às seguintes forças: 


a) F = Fo + kx 
b) F= je te 
c) F = Fo + kv 


PROBLEMAS a 


2.4 A força que atua numa partícula varia com a posição de acordo com a lei de potência 


F(z) = kz” 


a) Encontre a função energia potencial. 
b) Se v = vo em t = 0 e z = 0, encontre v em função de zx. 
c) Determine os pontos de retorno do movimento. 


Problemas 


2.5 Uma partícula de massa m está inicialmente em repouso. Uma força constante Fp atua 
na partícula até o instante t = tọ. A força então aumenta de valor linearmente com 
o tempo tal que depois de um intervalo adicional tọ a força vale 2Fọ. Mostre que a 
distância total percorrida pela partícula no intervalo total de tempo 2to é 


13 Fot 
6 m` 
2.6 Um bloco é lançado para cima em um plano inclinado com velocidade inicial vo. Se a 


inclinação do plano é 0, e o coeficiente de atrito cinético entre o plano e o bloco é u, 
encontre o tempo total necessário para o bloco retornar ao ponto de lançamento. 


2.7 Um bloco desliza numa superfície horizontal lubrificada com óleo pesado tal que o bloco 
sofre uma resistência viscosa que varia com a velocidade de acordo com a equação 


F(v) = —cv” 


Se a velocidade inicial é vo em t = 0, encontre v e o deslocamento x em função do 
tempo t. Encontre também v em função de x. Em particular, mostre que para n = 1/2, 
o bloco não percorrerá mais do que 
3 
Sma 
3c 


2.8 Uma partícula de massa m é abandonada em repouso a uma distância b de um centro 
fixo de força que atrai a partícula de acordo com a lei do inverso do quadrado 


F(x) = —kr? 
Mostre que o tempo necessário para a partícula chegar à origem é: 


ii 
8k 


T 


2.9 Encontre a relação entre a distância de queda e a velocidade de um corpo abandonado 
do repouso em queda livre e sujeito a resistência do ar proporcional: (a) à velocidade 
(b) ao quadrado da velocidade. 
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2.10 Um projétil é lançado verticalmente para cima com velocidade inicial vo. Supondo 
que a resistência do ar seja proporcional ao quadrado da velocidade, mostre que a 
velocidade do projétil ao retornar ao ponto de lançamento é 


Volt 
vg +V? 


no qual 
m 
v = velocidade terminal = ds 
c 
2.11 A velocidade de uma partícula de massa m varia com o deslocamento x de acordo com 
a equação 


yE 
x 


Encontre a força que atua na partícula em função de zx. 


2.12 Sabendo-se que a força que atua em uma partícula é o produto de uma função da 
posição e uma função da velocidade F(x, v) = f(x)g(v), mostre que a equação diferen- 
cial do movimento pode ser resolvida por integração. A equação de movimento pode 
ser resolvida por simples integração se a força for o produto de uma função da posição 
por uma função do tempo? E se for o produto de uma função do tempo por uma 
função da velocidade? 


2.13 A força que atua em uma partícula de massa m é dada por 
F = kvr 


em que k é uma constante. A partícula passa pela origem com velocidade vo em t = 0. 
Encontre x em função de t. 


2.14 Uma partícula executando movimento harmônico simples de amplitude 4 passa pela 
posição de equilíbrio com velocidade vo. Qual é o período de oscilação? 


2.15 Duas partículas de massas mı e ma executam movimento harmônico simples de ampli- 
tude A, e As, respectivamente. Qual é a razão entre seus períodos: T/T, se a energia 
total da partícula 1 for duas vezes a da partícula 2? 


2.16 Uma partícula executando movimento harmônico simples tem velocidade vı quando o 
deslocamento é x, e a velocidade v quando o deslocamento for £2. Encontre o período 
e a amplitude do movimento. 


2.17 Um único objeto de massa m está dependurado em duas molas de constantes kı e ko 
mantidas em posição vertical. Mostre que a frequência angular de oscilação é 


kı + k2 
m 
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se as molas forem ligadas em paralelo e 


kiko 
(kı + ko)jm 


se as molas forem ligadas em série. 


2.18 Uma caixa de massa M contendo um bloco de massa m está dependurada em uma 
mola de constante elástica k. Encontre a força de reação entre o bloco e o fundo da 
caixa em função do tempo sabendo-se que o sistema foi deslocado para baixo de uma 
distância d, contada a partir da posição de equilíbrio, e em seguida abandonado. Para 
que valor de d o bloco perderá o contato com o fundo da caixa na parte mais alta das 
oscilações verticais? Despreze a resistência do ar. 


2.19 Mostre que a razão entre duas elongações máximas sucessivas de um oscilador harmônico 
amortecido é constante. (Nota: Os máximos não ocorrem nos pontos de contato da 
curva de elongação em função do tempo com a curva Ae 7). 


2.20 Sabendo-se que a amplitude de um oscilador harmônico amortecido diminui para 1/e 
do seu valor inicial depois de n ciclos completos, mostre que a razão entre o período 
de oscilação e o período do mesmo oscilador sem amortecimento é dado por 


LN TRE 
To AN2 | Brn? 


2.21 A velocidade terminal de uma bola em queda livre é v;. Quando a bola está dependu- 
rada em uma mola leve a mola se distende de zọ. Mostre que a frequência natural de 
oscilação é dada por 


g P 


wi =- > -2 
zo 4v2 

supondo que a resistência do ar seja proporcional à velocidade. 
2.22 Mostre que a energia do sistema acima cai a 1/e do seu valor inicial em t = v;/9. 


2.23 Mostre que a frequência da força externa, para a qual a amplitude do oscilador 
harmônico forçado é a metade da amplitude de ressonância, é aproximadamente 


wo + yv3. 
Considere o caso de amortecimento fraco. 


2.24 Encontre a frequência para a qual a velocidade do oscilador harmônico forçado é 
máxima. (Sugestão: Maximize a quantidade vmar = wA(w)). 


2.25 Mostre que o fator de qualidade Q de um oscilador harmônico forçado é igual ao fator 
pelo qual deve ser multiplicado a resposta à frequência para se obter a resposta à 
frequência de ressonância. 
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2.26 Resolva a equação diferencial do movimento do oscilador harmônico sujeito a uma 
força amortecida da forma 
Fert = Foge“ cos(wt). 


2.27 Mostre que a série de Fourier para uma onda quadrada periódica é 
4 1 1 
F(t) = —(sen(wt) + 3 sen(3wt) + E sen(Swt) +...) 
T 
onde 


f(t) = +1 para 0< wt<r,2r < wt < 3r, etc 
f(t) = —1 para m< wt < ?2r,3r <wt < 4r, etc 
(2.68) 


2.28 Use o resultado acima para obter o movimento estacionário de um oscilador harmônico 
forçado por uma onda quadrada periódica de amplitude Fo. Em particular encontre 
as amplitudes relativas dos primeiros três termos A, A3 e A; da função resposta z(t) 
no caso em que o terceiro harmônico, 3w, coincide com a frequência de ressonância do 
oscilador. Considere o fator de qualidade Q = 100. 


Respostas de problemas selecionados de números ímpa- 
res 


— 3Foto. — 5 Foto 
21 v= “dt (Qto) o um 


2.3 aos y (2For + kgz?)m-t 
c)r= [kv — Foln (2) | 


Fo 


2.7 v = (A + Bt)” onde A = v”, B = 
x = Cuf — vê) onde C = 2 (2 = n) 


c 


(n—-1),a=(1-n) 
ent£1,2e8=2-n 


ses 
m 
—1 


2 jest O qa) 
2 


bre o duda) 


2.11 F(x) = —-mb?x’ 
2.13 x = atg (bt) onde a = 4/2mvo/k, b = /kvo/2m 
DO e 


Capítulo 3 


Movimento Geral de uma Partícula 
em Três Dimensões 


Passemos agora ao caso geral do movimento de uma partícula no espaço. 


3.1 Momentum Linear 


Vimos que a forma vetorial da equação de movimento de uma partícula é 
> dp 
Far 
dt 
ou, de modo equivalente, 
= — (my) (3.1) 


Essa é a forma abreviada para três equações nas quais as componentes de força podem 
depender das coordenadas, de suas derivadas em relação ao tempo e do tempo. Infeliz- 
mente não existe um método geral para encontrarmos soluções analíticas dessas equações 
em todos os casos possíveis. Existem entretanto muitos tipos de funções força, fisicamente 
importantes, para as quais as equações diferenciais do movimento podem ser resolvidas por 
métodos simples. Algumas dessas funções serão estudadas nas seções que se seguem. 

Nos casos em que F é uma função explícita do tempo, o momentum p pode ser encontrado 
determinando-se o impulso, isto é, por integração relativamente ao tempo, como no caso 
unidimensional: 


/ Fat = B) +C = m?(t) +C (3.2) 
Analogamente, uma segunda integração dará a posição: 
/ vdt = F(t) + C' (3.3) 


Embora o método seja perfeitamente válido, conhecermos de antemão a força como função 
do tempo não é uma situação típica na dinâmica. É claro que no caso especial da força ser 
nula, o momentum e a velocidade são constantes e as equações anteriores são válidas. Mais 
tarde teremos a oportunidade de discutir, em uma forma mais geral, o conceito de momentum 


To 
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constante sob ação de uma força nula, quando estudarmos sistemas de partículas no Capítulo 
4. 


3.2 Momentum Angular 


Consideremos a equação geral do movimento de uma partícula É = dp /dt. Multipliquemos 
ambos os membros pelo operador rx. Obteremos 


= dp 
rxF=rřx — 
dt 
O membro da esquerda da equação acima é, por definição, o momento da força F em relação 
à origem do sistema de coordenadas. O membro da direita é igual à derivada em relação ao 


tempo da quantidade r x p. Para provar isso, diferenciemos 


CEE E É EE 
— (rx p)= 7x x — 
nExBD=TxP+ixo, 
Mas vxp=vxmvy=myxv=bO. Logo, podemos escrever: 
Reco (E x P) (3.4) 
T > —1I . 
Fr 


A quantidade r x p é denominada momentum angular da partícula em relação à origem. 
Nosso resultado, em palavras, é o de que a taxa de variação do momentum angular de uma 
partícula é igual ao momento da força que atua na partícula. 

O conceito importante de momentum angular será particularmente útil no estudo de 
sistemas de partículas e corpos rígidos, Capítulos 4 e 5. 


3.3 O Princípio do Trabalho 


Na equação geral do movimento vamos tomar o produto escalar pela velocidade v em ambos 
os lados 


> dp d(myv 
PE et 
dt dt 
Usando a regra de diferenciação de um produto escalar, encontramos 
d(v . V V 
dv. y) 7. dy 
dt dt 
Considerando a massa constante, vemos que a equação acima é equivalente a 
> d /1 dT 
Fr-S(5 o) 3.5 
“Eua Ge) 


na qual introduzimos a energia cinética T = amv? Além disso, desde que vdt = dr, podemos 


integrar para obter 
[E dr= far (3.6) 


3.4. FORÇAS CONSERVATIVAS E CAMPOS DE FORÇAS TT 


O lado esquerdo da equação acima é uma integral de linha. Ela representa o trabalho 
realizado na partícula pela força F enquanto a partícula se movimenta em sua trajetória. O 
lado direito é a variação da energia cinética da partícula. Portanto a Equação (3.6) estabelece 
que o trabalho realizado sobre a partícula pela resultante da força é igual à variação da energia 
cinética. 


3.4 Forças Conservativas e Campos de Forças 


Geralmente, o valor de uma integral de linha, o trabalho neste caso, depende da trajetória 
de integração (veja Figura 3.1). Em outras palavras, o trabalho feito usualmente depende 
da rota específica que a partícula toma para ir de um ponto a outro. Isto significa que 
se tivermos de calcular o valor da integral de linha precisaremos conhecer de antemão a 
trajetória do movimento da partícula. Todavia, os tipos usuais de problemas que são de 
interesse em dinâmica da partícula são aqueles nos quais a trajetória do movimento não é 
conhecida a priori, ao contrário, a trajetória é uma da incógnitas a ser calculada. O princípio 
do trabalho expresso pela Equação (3.6) pode parecer, então, não ser muito útil aos nossos 
propósitos. Entretanto, na verdade este princípio é muito útil no estudo do movimento de 
uma partícula sob a ação de um tipo particular de força conhecida como força conservativa. 
Por sorte, muitas das forças fisicamente importantes são desse tipo. 


Figura 3.1: O trabalho feito por uma força Féa integral de linha f F. dr. 


Quando a força F for uma função das coordenadas de posição apenas, dizemos que ela 
define um campo de forças estático. Dentre os tipos possíveis de campos, existe uma classe 
importante para a qual a integral f F. dré independente da trajetória de integração. Tais 
campos de força são conservativos. Matematicamente, um campo conservativo é aquele para 
o qual a expressão F . dr é uma diferencial exata. Quando uma partícula move-se num 
campo conservativo, a integral de linha, e dessa maneira a variação da energia cinética, pode 
ser conhecida de antemão. Este conhecimento pode ser usado para predizer o movimento da 
partícula. 


3.5 A Função Energia Potencial para o Movimento Tri- 
dimensional 


Se uma partícula move-se sob a ação de uma força conservativa F, a afirmação de que o 
acréscimo do trabalho F . dr é uma diferencial exata significa que ela pode ser expressa 
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como a diferencial de uma função escalar da posição r, ou seja 
F . d? = —dV (T) (3.7) 


Esta equação é análoga àquela para o caso unidimensional F(x)dz = —dV, Seção 2.7. A 
função V é a energia potencial. O princípio do trabalho, Equação (3.6), é representada então 
de um modo simples por dT = —dV ou 


d(T +V) =0 (3.8) 


Esta relação implica que a quantidade T +V permanece constante durante o movimento da 
partícula. Chamamos esta constante de energia total E e escrevemos 


ma? +VO=-E (3.9) 

No caso de força não conservativa o incremento do trabalho não é uma diferencial exata e 
dessa maneira não pode ser igual a uma quantidade —dV. Um exemplo comum de força não 
conservativa é o atrito. Quando forças não conservativas estão presentes podemos expressar a 
força total como a soma F+F onde F representa forças « conservativas e F’ não conservativas. 
O princípio do trabalho é dado então por dT = F . dr+F' . dr = —dV +F . dT ou d(T +V) = 
F . dr, Vemos que a quantidade T +V não é constante, mas aumenta ou diminui quando a 
partícula se move dependendo do sinal de F . dr. No caso de forças dissipativas a direção 
de F' é oposta daquela de dr, desse modo FE. dré negativo e a energia total T + V diminui 
à proporção que a partícula se move. 


3.6 Gradiente e o Operador Del em Mecânica 


Em coordenadas retangulares, a equação 


é dada por 


Fado + Body + Bida = — E dz — dy — dz 
o que nos leva às seguintes igualdades 
F.=——— Fi = a F,=—— (3.10) 
Em palavras, se o campo de força for conservativo, então as componentes da força serão 


obtidas através das derivadas parciais da função energia potencial com o sinal trocado. 
Podemos escrever agora a força vetorialmente como: 


„oV 
— 


gs] 
| 
= 
| 
= 


(3.11) 
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Esta equação pode ser escrita numa forma abreviada mais conveniente 


F=-VV (3.12) 
Introduzimos aqui o operador vetorial 
> o ðo >09 
= FJ— +k 3.13 
E dx 1 Oy Oz 


V é chamado de operador del. A expressão VV é também chamada de O gradiente de 
V e algumas vezes é escrita grad V. Matematicamente, o gradiente de uma função é um 
vetor que representa sua derivada espacial. Fisicamente, o negativo do gradiente da função 
energia potencial nos dá a direção, o sentido e o módulo da força que atua numa partícula 
localizada num campo criado por outras partículas. O sinal negativo significa que a partícula 
é obrigada a mover-se na direção em que a energia potencial decresce. A Figura 3.2 é uma 
ilustração do gradiente. Nessa figura a função potencial é representada na forma de curvas 
de nível representando energia potencial constante. A força em qualquer ponto é sempre 
normal à superfície equipotencial que passa pelo ponto em questão. 


> 


Figura 3.2: Um campo de força representado por linhas de contorno de energia potencial. 


V = constante 


Rd 


3.7 Condições para a Existência de uma Função Po- 
tencial 


No Capítulo 2 vimos que em se tratando de movimento unidimensional, a força será conser- 
vativa sempre que ela for função apenas da posição. Estamos interessados agora em verificar 
se existe uma afirmação equivalente para o caso geral de movimentos bi e tridimensionais. 
Isto é, se a força que atua em uma partícula for uma função apenas das coordenadas de 
posição, existirá sempre uma função V que satisfaça a Equação (3.10) acima? A resposta 
para esta pergunta é não; somente existirá a função potencial se as componentes da força 
satisfizerem certos critérios. 
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Vamos supor que exista a função potencial, isto é, que as equações (3.10) sejam válidas. 


Temos então 
oF; 8V oF J oV 


dy “OyO dr “Oxdy 


Como a ordem de diferenciação pode ser invertida, as duas expressões são iguais e, 


dr,  0F, ôF, ƏF, oF, OF, (3.14) 
dy ôx dz Or dz Oy 


Estas são condições necessárias para que a função potencial exista; elas expressam a condição 
de que 
F . dr = F,dx + Fdy+ Fodz 


seja uma diferencial exata. É possível mostrar também que elas são condições suficientes, 
isto é, se as Equações (3.14) forem válidas em todos os pontos, então as componentes da 
força serão realmente deriváveis de uma função potencial V (x,y,z), e a soma da energia 
cinética com a energia potencial será constantel. 

Os critérios para que um campo de força seja conservativo serão convenientemente apre- 
sentados em termos do operador del. Nesta aplicação introduzimos o produto vetorial do 
operador del: 


~ no Z (E 2) as 


Oz Ox 


O produto vetorial definido acima é chamado de rotacional de F. De acordo com as Equações 
(3.14), vemos que as componentes do rotacional se anulam se a força F for conservativa. 
Então a condição para uma força ser conservativa pode ser dada numa forma compacta 


VxF=0 (3.16) 


Matematicamente, esta condição exige que a expressão F.d seja uma diferencial exata, 
ou em outras palavras, que a integral f F. d seja independente do caminho de integração. 
Fisicamente, a anulação do rotacional da força F significa que o trabalho feito por F sobre 
a partícula em movimento é independente da trajetória que liga um ponto ao outro. 

Há uma terceira expressão envolvendo o operador del, ou seja o produto escalar V.E. 
Esse é chamado a divergência de F. No caso de campo de força, a divergência dá uma 
medida da densidade de fontes de campo em um certo ponto. A divergência é de particular 
importância na teoria de eletricidade e magnetismo. 

Expressões para o gradiente, rotacional e divergência em coordenadas cilíndricas e esféricas 
são apresentadas abaixo: 


Coordenadas Cilíndricas 


1Veja qualquer livro texto de cálculo avançado, por exemplo, A. E. Taylor, Advanced Calculus, Ginn, 
Boston, 1955. Uma discussão interessante dos critérios de conservação quando este campo contém singula- 
ridades foi feita por Feng na Amer. J. Phys. 37, 616 (1969). 
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Gradiente 
oV 18V 


Cp 4 C, 4 


Wea a ES 


Rotacional 


sm (10F, DRN a a Ore 0FN. ORE Toa 
AR ESS E (x dB Rô)! 


Divergência 


1ô(RFr) | 10F, OF, 
R ôR ` oz 


vV.F= 


Coordenadas Esféricas 


Gradiente əv 19V 1 əv 
PE Aro “rseng dp * 
Rotacional 
E 1 O(Fosend) Oy, 1/1 ôF, Ori. 
F = — — | 
VR r sen 0 o0 e) ” r (= oQ Or di 
1 olr Fo) = oF, 25 
rl ôr a)” 
Divergência 
g F- 1 (r? F.) 1 sendo) 1 09F 
Cr ðr ' rsenð 90 “rsen6 d4 
Exemplos 


1. Encontre o campo de força da função potencial V = z? + «xy + zz. Aplicando o 
operador del. . F 
F =-VV = —T(2r +y +2)—-Jz-krzr 


2. O campo de força F = Ixy +Jxz + kyz é conservativo? O rotacional de F é 
To SE 
VES Ž E l=I(z-2)+]70+k(2— 7x) 
TY TZ YZ 


A expressão final não é nula, logo o campo não é conservativo. 
3. Para que valores das constantes a, b e c a força F = T(ax + by?) +Jcxy é conservativa? 
Tomando o rotacional de F temos 


T po Bi 
VxF= É do à |=k(c-2by 
ax +by? cry 0 
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O que mostra que a força será conservativa se c = 2b. O resultado é independente do valor 
de a. 

4. Mostre que a lei de força inverso do quadrado em três dimensões F = -5ē, é conser- 
vativa calculando o rotacional. Use coordenadas esféricas. 

O rotacional é dado na página anterior. 


1 E Cor Ersen 

E ð ə a 
Rip e r2seng | 9 2 3p 

F, rh, rFọsenð 


Temos F, = -5, Fo = 0, F = 0. O rotacional se reduz então a 
ana e e o) 
F = — — — =0 
IR rsend dp \ r? r oO r? 


que, evidentemente, se anula porque ambas as derivadas parciais são nulas. Então a força 
em questão é conservativa. 


3.8 Forças do Tipo Separável 


Muitas vezes, escolhendo adequadamente o sistema de coordenadas as componentes de um 
campo de força dependem apenas das coordenadas respectivas, isto é 


F =TF,(£) +7F,(y) + K Fi(2) (3.17) 


Forças deste tipo são chamadas de separáveis. Verificamos que o rotacional de tais forças 
é nulo e então, o campo é conservativo independentemente das formas particulares das 
componentes de força desde que cada uma seja função apenas da coordenada envolvida. 
A integração das equações diferenciais do movimento é então muito simples, porque cada 
equação componente é do tipo më = F(x). Nesse caso as equações podem ser resolvidas 
pelos métodos descritos no capítulo anterior de movimento retilíneo. 

Se as componentes da força envolverem o tempo e a derivada das coordenadas com relação 
ao tempo, então a força não é mais necessariamente conservativa. Apesar disso, se a força for 
separável, então as componentes da equação do movimento serão da forma më = F(x, x,t) e 
podem ser resolvidas pelos métodos usados no capítulo anterior. Alguns exemplos de forças 
separáveis, conservativas e não conservativas, serão discutidos nas seções seguintes. 


3.9 Movimento de um Projétil em um Campo Gravi- 
tacional Uniforme 
Um dos mais famosos problemas clássicos de dinâmica é o de movimento de um projétil. Es- 


tudaremos este problema com alguns detalhes porque ele ilustra a maior parte dos princípios 
estudados nas seções anteriores. 


Sem Resistência do Ar 
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Primeiramente, por simplicidade, consideremos o caso de um projétil movendo-se sem 
resistência do ar. Nessa situação idealizada apenas a força da gravidade atua no projétil. 
Se fizermos o eixo z coincidir com a vertical, a equação diferencial do movimento ficará 


m— = —mgk 


Se idealizarmos mais ainda o problema e considerarmos que a aceleração da gravidade g é 
constante, então a função força é visivelmente do tipo separável e é também conservativa 
desde que é um caso especial daquele expresso pela Equação (3.17). Particularizaremos o 
problema considerando que no instante t = 0 o projétil estava na origem com velocidade vo. 
A equação da energia (3.9) fica 


ou, 
v = vê — 292 (3.18) 


dando a velocidade como uma função da altura. Esta é toda a informação que podemos 
obter da equação de energia. 
Para conseguir mais informações, devemos voltar à equação diferencial do movimento. 


Podemos escrevê-la 
E a DR 
ala)? 


que é da forma discutida na Seção (1.20). Pode ser integrada diretamente. Uma primeira 
integração nos dá a velocidade. 


dr > 
— = —gtk + V 
dt g 0 
onde a constante de integração Vo é a velocidade inicial. Outra integração nos leva ao vetor 
posição 
1 e 
r= -39ťk + Vot (3.19) 


A constante de integração rọ, neste caso, é zero porque a posição inicial do projétil era a 
origem. As componentes da equação acima são 


T = Tot 
= Yot 


Lo 
ZE UE Ca 


Aqui to, Wo € Zo são as componentes da velocidade inicial vo. Resolvemos então o problema 
de determinar a posição do projétil como função do tempo. 
Para conhecermos a trajetória do projétil, eliminamos o tempo t das equações x e y e 
obtemos 
dba 
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onde a constante b é 


Isto significa que a trajetória está contida em um plano. Em particular, se y = 0, então a 
trajetória está contida no plano zz. Eliminando o tempo nas equações x e z, encontramos 
que a equação da trajetória é da forma 


z = az — Br? 


onde a = žo/čto e 8 = g/2t8. A trajetória é, deste modo, uma parábola contida no plano 
y = bz. A Figura 3.3 mostra isto. 


Figura 3.3: Trajetória de um projétil movendo-se em três dimensões. 


Resistência do Ar Linear 


Vamos considerar o movimento de um projétil para a situação mais realista onde há uma 
força retardadora devido a resistência do ar. Neste caso o movimento é não conservativo. A 
energia total diminui continuamente como resultado da perda por fricção. 

Por simplicidade, vamos supor que a lei de resistência do ar seja linear de tal maneira que 
a força resistiva varia diretamente com a velocidade v. Será conveniente escrever a constante 
de proporcionalidade como my onde m é a massa do projétil. Temos então duas forças 
atuando no projétil; a resistência do ar myy, e a força da gravidade —mgk, como antes. À 
equação diferencial do movimento é então 


dr n > 
mop RE mgk 
ou, depois de cancelar os m's, 
dr E: = 
o q gh 


dt 
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A integração desta equação é feita convenientemente expressando-a na forma de componentes 
como segue: 


ť = —yt 
= =Yyy 
Z = —y2-g9 


Vemos que as equações estão separadas. Cada uma então pode ser resolvida separadamente 
pelos métodos do capítulo anterior. Usando nossos resultados da Seção 2.9, podemos escrever 
prontamente as soluções 


q = toe 
= ne” (3.20) 
Z 


Ran g —yt 
= ge Eq pr 
(1-0) 
para as componentes da velocidade, e 


g = “a e) 


= Yo = e 
Da (1 ) (3.21) 


a e (2+5) (1- er) = 
É O Y 
para as coordenadas de posição. Aqui, como antes, as componentes da velocidade inicial são 
čo, Yo € Žo e consideramos o projétil inicialmente na origem. 
Como no caso de resistência do ar nula, o movimento está inteiramente contido no plano 
y = br com b = yo/io. Todavia, a trajetória neste plano não é uma parábola, mas uma 
curva que se situa abaixo da parábola correspondente. Veja a ilustração na Figura 3.4. 
Inspecionando as equações de x e y vemos que, para t grande, os valores de x e y se aproximam 
dos valores limites E . 
ss gado 
YJ Y 
Isto significa que a trajetória tem uma assíntota vertical. 
A solução final do movimento de um projétil com resistência do ar linear, veja Equações 
(3.21) pode ser escrita vetorialmente assim 


e a 
r= (2+2) (1.00) -E2 


Por diferenciação vemos que esta equação é realmente solução da equação diferencial do 
movimento. 

É útil considerar o caso em que a resistência do ar for muito pequena, isto é, quando o 
valor da quantidade yt nos fatores exponenciais for muito menor do que a unidade. Para 
isso vamos desenvolver em série o fator exponencial, usando 
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trajetoria sem 


y resistencia do 
ar 
O A 
trajetoria com ' assintota 
resistencia do -© vertical 
ar 


Figura 3.4: Comparação das trajetórias de uma partícula movendo-se com e sem resistência 
do ar. 


onde faremos u = —yt. O resultados, depois de cancelar e colecionar os termos, fica na 
forma 1 
E= Vot — zIËk — AF 


» t ytê v tê ytt 


A quantidade Ar pode ser entendida como uma correção na trajetória do projétil sem re- 
sistência que leva à trajetória verdadeira. 

No movimento real de um projétil na atmosfera, a lei de resistência do ar não é de 
maneira nenhuma linear, mas uma função muito complicada da velocidade. Um cálculo 
preciso da trajetória pode ser feito por meio de métodos de integração numérica utilizando- 
se computadores. 


onde 


ARS 


3.10 O Oscilador Harmônico em duas e três dimensões 


Consideramos nessa seção o movimento de uma partícula que está sujeita a uma força res- 
tauradora linear que está sempre dirigida para um ponto fixo, a origem do nosso sistema de 
coordenadas. Esta força pode ser representada pela expressão 


F = kr 
A equação diferencial do movimento é então expressa por 
d2r 


ma A (3.22) 
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A situação pode ser visualizada, de uma maneira aproximada, por uma partícula presa a 


um conjunto de molas como mostrado na Figura 3.5. Este é o oscilador harmônico tri- 
dimensional. 


04 — 
24 


00000000 
Q 
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= 


Figura 3.5: Modelo de um oscilador harmônico tridimensional. 


O oscilador bi-dimensional 
Para o caso de movimento num plano, a equação diferencial é equivalente às duas equações 
componentes 
mł = —kx 
mij = —ky 


Elas estão separadas, e podemos, imediatamente, escrever as soluções na forma 


x = Acos(wt + a) y = Beos(wt + 8) (3.23) 


onde 


k 
W E= E 
m 


As constantes 4, B, a e 5 ficam determinadas a partir das condições iniciais, em qualquer 
caso. 


Encontramos a equação da trajetória eliminando o tempo t nestas duas equações. Para 
conseguirmos isto, escreveremos a segunda equação na forma 


y = Bcos(wut +a + A) 
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onde À = 8 — a. Desenvolvendo, temos 
y = Blcos(wt + a) cos A — sen(wt + a) sen A] 


Usando a primeira das Equações (3.23), obtemos 


Di aÃ ne 
DT gÂ 1 E sen A (3.24) 


que é uma equação quadrática em x e y. À equação quadrática geral 


ax” + bry + cy + de + ey = f 
pode representar uma elípse, uma parábola, ou uma hipérbole, dependendo do discriminante 
b — 4ac 


ser negativo, zero, ou positivo, respectivamente. No presente caso o discriminante é — (2 sen A/AB)?. 
A Figura 3.6 mostra que a trajetória é elíptica pois o discriminante da equação é negativo. 

No caso particular em que a diferença de fase A seja igual a 7/2 a equação da trajetória 
se reduzirá a 


que é a equação de uma elípse cujos eixos coincidem com os eixos coordenados. Por outro 

lado, se a diferença de fase for 0 ou 7, então a equação da trajetória se reduzirá a uma linha 
b) 2 

reta dada por 


O sinal positivo é válido quando A = 0, e o negativo para À = 7. E possível mostrar que 
no caso geral o eixo da trajetória elíptica formará um ângulo à com o eixo x definido por 


o 2AB cos À 


Rr (3.25) 


tg 29 


Como exercício, deixamos para o leitor a dedução da expressão acima. 


O Oscilador Harmônico Tri-dimensional 


Para o caso de movimento tridimensional a equação do movimento é equivalente às três 
equações 
mk = —ka my = —ky mz = —kz 


que estão separadas. As soluções são da forma (3.23), isto é 
Aı sen wt + B4 cos wt 


A, sen wt + Bo cos wt (3.26) 


As sen wt + Bs cos wt 


x e 
| 
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Figura 3.6: Trajetória elíptica do movimento de um oscilador harmônico bi-dimensional. 


Determinam-se as seis constantes de integração a partir da posição e velocidade iniciais da 
partícula. As Equações (3.26) escritas na forma vetorial ficam 


r= Asenwt + B cos wt 


onde as componentes de A são Ay, A e As e da mesma forma para B. Fica claro que o 
movimento se faz inteiramente em um único plano que é o plano comum aos dois vetores 
constantes À e B, e que a trajetória da partícula nesse plano é uma elípse, como no caso 
bidimensional. A análise relativa à forma da trajetória elíptica feita no caso bidimensional 
também se aplica no caso tridimensional. 


Oscilador não Isotrópico 


Na discussão acima consideramos o oscilador tridimensional isotrópico, onde a força res- 
tauradora era independente da direção de deslocamento. Se a força restauradora depender 
da direção de deslocamento, temos o caso do oscilador não-isotrópico. Escolhendo adequa- 
damente os eixos do sistema de coordenadas, as equações diferenciais do movimento para o 
oscilador não isotrópico são 


mz = —kx 
my = —koy (3.27) 
mz = —ka2 


Aqui temos um caso de três frequências de oscilação diferentes: 
| kı | kə | k3 
Wi = = W2 = === W3 = Ea 
m m m 


q = Acos(wt+a) 


e as soluções são 
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y = Bcos(wst+ 5) (3.28) 
z = Ccos(wst+9) 


Também nesse caso, as seis constantes de integração das equações acima ficarão determi- 
nadas a partir das condições iniciais. A oscilação resultante da partícula fica inteiramente 
contida dentro de uma caixa retangular (cujos lados são 24, 2B e 2C) centrada na origem. 
Se as frequências w1, W2 € w3 forem redutíveis a uma medida comum, isto é, se 

Ef ia (3.29) 

nı Na n3 
onde n1, no e nz são inteiros, a trajetória, chamada de Figura de Lissajous, se fechará, porque 
depois de um tempo 27n1/w, = 2mn5/w, = 2mnz/w3 a partícula retornará à sua posição 
inicial e o movimento se repetirá. (Na Equação (3.29) cancelamos os fatores inteiros comuns). 
Por outro lado, se os w's não forem redutíveis a um fator comum, a trajetória não será 
fechada. Neste caso dizemos que a trajetória pode encher completamente a caixa retangular 
mencionada acima, pelo menos no sentido de que se esperarmos um tempo suficientemente 
longo, a partícula chegará arbitrariamente perto de qualquer ponto dado. 

Em muitos casos a força restauradora líquida exercida em um dado átomo numa substância 
cristalina sólida é aproximadamente linear nos deslocamentos. As frequências de oscilação 
resultantes normalmente ficam na região infravermelho do espectro: 10!? a 10!4 vibrações 
por segundo. 


3.11 Movimento de Partículas Carregadas em Campos 
Elétricos e Magnéticos 


Quando uma partícula eletricamente carregada estiver na vizinhança de outras partículas 
carregadas, ela sofrerá a ação de uma força. Dizemos que esta força F é devida ao campo 
elétrico E que tem origem nestas outras cargas. Escrevemos 


F=qÊ 


onde q é a carga elétrica da partícula em questão. A equação de movimento da partícula é 
então 


m—- = qE 3.30 
ps (3.30) 
ou, em componentes 

m qEzx 

my = qE 

mž = qE; 


As componentes do campo são, em geral, funções das coordenadas de posição x, y e z. No 
caso de campos que variam com o tempo (isto é, se as cargas que produzem E estiverem em 
movimento) as componentes, evidentemente, também envolverão t. 
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Vamos considerar um caso simples, ou seja aquele de um campo elétrico uniforme. Esco- 
lhemos um dos eixos, digamos o eixo z, na direção do campo. Então E, = E, = 0, e E = E,. 
As equações diferenciais do movimento de uma partícula carregada neste campo são 


E 
= y=0 z=  — constante 
m 


Estas equações são identicas às equações para um projétil em um campo gravitacional uni- 
forme. À trajetória, portanto, é uma parábola. 
Os livros textos de Eletromagnetismo mostram que 


VxE=0 


se É for devido a cargas estáticas. Isto significa que o movimento em tal campo é conservativo, 
e que existe uma função potencial q onde E = Vo. A energia boreari] de uma partícula 
de carga q neste campo é q¢, e a energia total é constante e igual a imo’ + qo. 

Na presença de um campo magnético estático (B) a força atuante em uma carga móvel 
é convenientemente expressa pelo produto vetorial 


F = q(¥ x B) (3.31) 


onde v é a velocidade e q é a carga. A equação diferencial do movimento de uma partícula 
movendo-se em um campo puramente magnético é então 


dY un 
m= q(¥ x B) (3.32) 


Desta equação concluimos que a aceleração é sempre perpendicular à direção do movimen- 
to. Isto significa que a componente tangencial da aceleração (ù) é nula, e logo, a partícula 
move-se com o módulo da velocidade constante. Este fato é verdadeiro mesmo quando B é 
uma função dependente da posição r, desde que essa função não varie com o tempo. 


Exemplo 


Vamos examinar o movimento de uma partícula carregada em um campo magnético 
constante uniforme. Vamos escolher o eixo z na direção do campo, isto é, 


B=kB 
A equação diferencial do movimento é 
er o n TTE 
TT (WxB)=9B|& | 2 
00 1 


m(Ti+Jj+k2) =qB(ry-— Jt) 
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Igualando as componentes temos, 


mx = qBy 
my = —qBt (3.33) 
2 =n 


Aqui, pela primeira vez, encontramos um conjunto de equações diferenciais que não são do 
tipo separáveis. A solução é relativamente simples, todavia, porque podemos integrar uma 
única vez no tempo e obter 


mč = gBy+ca 


my = —qBr +c 
2 = constante = žo 
ou 
ït =wy + C1 y = —wa + C2 Z= žo (3.34) 


onde usamos a abreviação w = qB/m. Os C’s são constantes de integração, e Cy = a/m, 
C2 = c&2/m. Usando a expressão de y da segunda das Equações (3.34), na primeira das 
Equações (3.33), obtemos uma equação separada para x, 


+w’ r =w a (3.35) 
onde a = C5/w. A solução é 
xz =a + Acos(wt + 00) (3.36) 


onde A e fo são constantes de integração. Diferenciando esta equação com relação ao tempo, 
temos 
t = — Awsen(wt + bo) (3.37) 


Esta expressão pode ser substituida no lado esquerdo da primeira das Equações (3.34) e a 
equação resultante resolvida para y. O resultado é 


y =b — Asen(wt + 06) (3.38) 


onde b = C4 /w. Para encontrarmos a forma da trajetória do movimento, eliminamos t entre 
as Equações (3.36) e (3.38) para obter 


(z-a? +(y-b? = Aº (3.39) 


A projeção da trajetória do movimento no plano xy é um círculo de raio A centrado 
no ponto (a,b). Como da terceira equação (3.34) a velocidade na direção z é constante, 
concluimos que a trajetória é uma helicóide. O eixo da hélice coincide com a direção do campo 
magnético, como pode ser visto na Figura 3.7. Derivando a Equação (3.38) encontramos 


y = — Aw cos(wt + 00) (3.40) 

Eliminando o tempo das Equações (3.37) e (3.40) obtemos 
PHY = Aw = A? (2) 
m 


(3.41) 
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2 
e xX 


Z 


Figura 3.7: Trajetória helicoidal de uma partícula carregada movendo-se em um campo 
magnético. 


chamando vı = yt? + ġ?, vemos que o raio A da hélice fica 


A Aa (3.42) 
wW qB 
Se a velocidade não tiver componente na direção z, a trajetória será um círculo de raio 
A. Fica evidente que 4 é diretamente proporcional ao módulo da velocidade vı, e que a 
frequência angular w do movimento na trajetória circular é independente da velocidade. w 
é conhecida como a frequência do cíclotron. O cíclotron, inventado por Ernest Laurence, 
depende para sua operação do fato que w é independente da velocidade. 


3.12 Movimento Vinculado de uma Partícula 


Quando uma partícula móvel estiver geometricamente presa no sentido de que deverá per- 
manecer em uma superfície ou curva definida, diremos que o movimento é vinculado. Um 
pedaço de gelo deslizando em uma superfície curva, ou uma bolinha furada deslizando presa 
em um arame, são exemplos de movimento vinculado. O vínculo pode ser completo, como 
na bolinha ou pode ser parcial como no primeiro exemplo. Vínculos podem ser fixos, ou 
podem ser móveis. Neste capítulo estudaremos apenas vínculos fixos. 


A Equação de Energia para Vínculos Lisos 


A força resultante em uma partícula que se move sob a ação de um vínculo pode ser 
expressa como a soma vetorial da força externa F e da força de vínculo R, onde R é a reação 
do agente vinculador sobre a partícula. A equação de movimento pode ser escrita como 


do a 
E era 3.43 
my + (3.43) 
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Multiplicando escalarmente a equação acima por v, temos 


dv 


ma V=P.v+R.? (3.44) 


No caso de vínculo liso — por exemplo, uma superfície sem atrito — a reação R é normal à 
superfície ou curva enquanto a velocidade v é tangente à superfície. Logo, R é perpendicular 
a ve o produto escalar R . v se anula. A Equação (3.44) se reduz portanto a: 


Ai R 
algm s) =F. y 
als 


Consequentemente, se F for conservativa, poderemos integrar como na Seção 3.5, e obteremos 
a mesma relação de energia que a Equação (3.9), isto é 


1 

amu? + V(x,y,2) = constante = E 
Então apesar de permanecer na superfície ou na curva, a partícula move-se de modo a manter 
a energia total constante. E claro que esperávamos que isto acontecesse para vínculos sem 
atrito. 


Exemplo 


Coloca-se uma partícula no topo de uma esfera lisa de raio a. Se a partícula for ligeira- 
mente perturbada, em que ponto ela abandonará a esfera? 

As forças atuantes na partícula são a força da gravidade, para baixo, e a reação R da 
superfície esférica. A equação de movimento é 


> 


dv > 
m— =m} +R 
T 
Vamos escolher eixos coordenados como os mostrados na Figura 3.8. A energia potencial é 


então mgz, e a equação da energia é 
l o 
gz +mgz =E 


Usando as condições iniciais (v = 0 para z = a) encontramos E = mga, logo, quando a 
partícula desliza, sua velocidade é dada por 
Qua 
v“ = 2g(a — z) 
Tomando a componente radial da equação do movimento, podemos escrever a equação de 


força assim 


2 
ZLA —mgcos0 + R = -mg +R 
a a 


Então 5 
zZz mw z 


R=mg- — — = mga — Jya — 2) = DS (32 — 29) 
a a a a a 
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mg 


y 


Figura 3.8: Forças atuantes em uma partícula que desliza em uma esfera lisa. 


Logo R se anula quando z = Za e neste ponto a partícula abandona a esfera. Isto poderia 


também ser conseguido observando que R muda de sinal neste ponto. 


Movimento em uma Curva 


No caso em que uma partícula estiver vinculada a mover-se em uma certa curva, a equação 
de energia juntamente com as equações paramétricas da curva 


v=a(5) y = y(S) z = z(5) 


são suficientes para determinar o movimento. (O parâmetro S é a distância medida ao 
longo da curva a partir de algum ponto de referência arbitrário). O movimento pode ser 
determinado considerando-se que a energia potencial e a energia cinética podem ser escritas 
em termos do parâmetro S apenas; V(S) e mS2, respectivamente. Dessa maneira a equação 
de energia pode ser escrita 
1 so 
mS +V(S) =E 

2 
de onde S (e dessa maneira x, y e z) pode ser obtido por integração. Fazendo o caminho 
contrário, diferenciando-se a equação acima com relação a t e cancelando o fator comum S, 
obtemos a seguinte equação diferencial do movimento para a partícula: 


=- dV 
fi 3.45 
ms + T (3.45) 


Esta equação pode ser colocada na forma equivalente 
ms — F Ss = 0 


onde Fs é a componente da força externa F na direção S. Isto significa que 
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3.13 O Pêndulo Simples 


O pêndulo simples ilustra bem as considerações acima. Uma partícula de massa m presa na 
extremidade de uma haste ou corda leve e inextensível movendo-se em um plano vertical, 
constitui o pêndulo simples. Tal pêndulo é também dinamicamente equivalente a uma bolinha 
furada deslizando presa a um fio com a forma de um aro vertical. 

Na Figura 3.9, 0 é o ângulo entre a vertical e a linha CP onde C é o centro da trajetória 
circular e P é a posição instantânea da partícula. A distância S é medida a partir da posição 
de equilíbrio O. Vemos, da figura, que a componente F's da força da gravidade mg na direção 
de S é igual a —-mgsen6. Se l for o comprimento do pêndulo, então 0 = S/l. A equação 
diferencial de movimento então fica assim 


Figura 3.9: O Pêndulo Simples. 


mõ + mg sen a =0 


ou, em termos de 0, ficaria 


à + send =0 


A energia potencial V pode ser escrita como mgz onde z é a distância vertical da partícula 
medida a partir de O, isto é, 


V = mgz = mgl(1 — cos0) = mgl — mgl cos (5) 


Aqui, então, 
dV (5) 9 
“45 = —mgsen X = mgsen6 = Fs. 


Para conseguirmos uma solução aproximada da equação diferencial de movimento, consi- 
deremos que 0 permaneça pequeno. Neste caso 


sen6 = 0 
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Logo 


0+%9=0 


esta é a equação diferencial do oscilador harmônico. A solução, como já vimos na Seção 2.12 
é 
0 = bo cos(wot + do) 


onde wo = 4/ g/l, ło é a amplitude de oscilação e ġo é um fator de fase. Dentro dos limites 
de validade dessa aproximação o movimento é harmônico simples e o período de oscilação é 


dado por: 
27 l 
To = — = 2m4 — (3.46) 
Wo g 


que é uma fórmula elementar bem conhecida. 


3.14 Solução mais Precisa do Pêndulo Simples e o Os- 
cilador não Linear 


A equação diferencial do movimento do pêndulo simples 


E g 
0+ T sen = 0 
é um caso especial de equação diferencial geral para o movimento sob a ação de uma força 
restauradora não linear, isto é, uma força que varia de maneira diferente da linear. 
A equação do problema unidimensional geral sem amortecimento pode ser escrita 


EFI (3.47) 
onde é é a variável que representa o deslocamento a partir da posição de equilíbrio, tal que 
HO) =0 


Equações diferenciais não lineares normalmente exigem algum método de aproximação para 
suas soluções. Suponha que a função f(£) seja desenvolvida como uma série de potência, 


FCE) = mê + ag? + asg? +... 


e levada na Equação (3.47) que se transformaria em 


2 
E | a1 | a£? | azé’ H... =0 (3.48) 


que é a expansão da equação geral do movimento do oscilador não linear sem amortecimen- 
to. O termo ajé nesta equação é o termo linear. Se este termo for predominante, isto é, 
se ay for muito maior do que os outros coeficientes, então o movimento será aproximada- 
mente harmônico com frequência angular ,/aj. Uma solução mais precisa deve levar em 
consideração os termos restantes não lineares. 
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Como ilustração, vamos voltar ao problema do pêndulo simples. Se usarmos o desenvol- 
vimento em série 
6 ge 
Ca TES 


e se retivermos apenas os dois primeiros termos, obteremos 


A g g 3 
Z — .49 


como uma segunda aproximação para a equação diferencial do movimento. Sabemos que o 
movimento é periódico. Vamos tentar uma solução na forma de uma função senoidal simples 


0 = Acoswt 
Inserindo-a na equação diferencial obteremos 


— Aw? cos wt + ZA cos wt — ZA cos? wt = 0 


ou, depois de usar a identidade trigonométrica 


cos? u = é cos u + : cos gu 
4 4 


teremos, depois de reagrupar os termos 


Excluindo o caso trivial A = 0, vemos que a equação acima não é válida para todos os valores 
de t. Dessa maneira nossa função tentativa A cos wt não pode ser uma solução. Observando 
que o termo em cos 3wt aparece na equação acima somos levados a tentar uma solução na 
forma 

0 = Acoswt + B cos 3wt (3.50) 


que deverá ser uma aproximação melhor do que Acoswt. Verifica-se que isto é verdade. 
Levando esta solução na Equação (3.49), encontramos, depois de um tratamento idêntico ao 
que foi feito anteriormente 


Aê Aê 
(~a + TA — e) cos wt + (o8. + ZB = e) cos 3wt + (termos em 


potências mais altas de B e múltiplos mais altos de wt) = 0 


Novamente a equação não é válida para todos os valores de t, mas nossa aproximação será 
razoavelmente precisa se os coeficientes dos dois primeiros termos em cosseno puderem se 
anular separadamente: 

gA’? gA? 


E E e a e 
pe: gi 0 SDE 2 0 


3.15. SOLUÇÃO EXATA DO MOVIMENTO DO PÊNDULO SIMPLES POR MEIO DE INTEGRAIS E 


Da primeira equação 


A? 
w? = 7 h a Z) (3.51) 
Levando este valor de w? na segunda equação, encontramos 
1 A’ 
BA 


3(64 +274?) 192 
Observando a Equação (3.50) vemos que a amplitude 6o da oscilação do pêndulo é dado por 
A? 
bo = A+ B= A- — 
TET 192 
ou, se 4 for pequeno 
o = A 


O significado da Equação (3.51) pode agora ser entendido. A frequência de oscilação depende 
da amplitude ĝo. Na verdade, podemos escrever 


e, para o período 


T= Fa 27 TE 

j (1— 500) V9 

l 1 
2 14 at.) 
mji ( 160 


1 
To (1 + 1690 +.. ) (3.52) 


Q 


Q 


onde To é o período para amplitude nula. 

Esta análise, apesar de ser grosseira, levanta dois fatos essenciais na oscilação livre sob 
a ação de força restauradora não linear; isto é, o período de oscilação é uma função da 
amplitude de vibração, e a oscilação não é exatamente senoidal mas pode ser considerada 
como uma superposição de uma mistura de harmônicos. Pode ser mostrado que a vibração 
de um sistema não linear forçado por uma força puramente senoidal também será distorci- 
do; isto é, conterá harmônicos. O auto-falante de um rádio ou um sistema de “som”, por 
exemplo, pode introduzir distorção (harmônicos) acima e abaixo daquela produzida pelo 
sistema amplificador eletrônico. 


3.15 Solução Exata do Movimento do Pêndulo Simples 
por Meio de Integrais Elípticas 


Usando a equação de energia potencial do pêndulo simples podemos escrever a equação de 
energia na forma 


z (10) + mgl(1 — cos 0) = E (3.53) 
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Se o pêndulo for deslocado de um ângulo ĝo (a amplitude) e abandonado (fo = 0), então 
E = mgl(1 — cos ĝo). Levando este valor na equação de energia esta se reduzirá a 


: 2 
8? = Feos 0 — cos 60) (3.54) 


Usando a identidade cos 0 = 1 — 2 sen 2 (8 y podemos escrever 


3 = 2g 200 an2? 
0 = 7 [sen z — sen 5) (3.55) 


é conveniente descrever o movimento em termos da variável à definida pela equação 


sen | 5 
(G) 1º 
send = sen (8) = Sen; (3.56) 


Derivando em relação ao tempo, obtemos 


Ed 090 
(cos ġ)ọ = g Cos (5) 5 (3.57) 


Usando as Equações (3.56) e (3.57) podemos prontamente escrever a equação (3.55) em 
termos de à 


É? = (1 — k? sen? 6) (3.58) 


separando as variáveis e integrando a relação entre q e t 


l 
= oh ds; e sen? À g ENO) (3.59) 


A função F(k, 6) = Jé(1— k? sen? 4) ? à do é conhecida como integral elíptica incompleta do 
primeiro tipo. O período do pêndulo é obtido notando-se que O aumenta de O até ĝo em um 


quarto de ciclo. Vemos então que q varia de 0 até 5 no mesmo intervalo de tempo. Dessa 
maneira podemos escrever para o período T 


Ta f =” nji KW (3.60) 


A função K(k) = fè (1 — k2sen2 o) ido = F(k, 7) é conhecida como a integral elíptica 
completa do primeiro tipo. Os valores das integrais elípticas são tabelados. Podemos obter 
uma expressão aproximada desenvolvendo o integrando da Equação (3.60) pelo teorema 


binomial e integrando termo a termo. O resultado é 


T= tp (14 Sono. Dao ami (14804. ) (3.61) 
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Para pequenos valores de 06, temos 


o que nos leva à expressão aproximada 


i 02 


que concorda com o valor de T encontrado na seção anterior. 


Exemplo 


Encontre o período de um pêndulo simples oscilando com uma amplitude de 20°. Use 
tabelas de funções elípticas, e também compare com os valores calculados utilizando-se as 
aproximações descritas anteriormente. 

Para uma amplitude de 20°, k = sen 10° = 0,17365, e 0o/2 = 0,17453 radianos. Os 
resultados ficam 


De tabelas e da Equação (3.60) T = 4/2 K(10°) = Vi (6, 3312) 


Da Equação (3.61) T=2m/i (1+ jsen?10º) = ,/Í (6, 3306) 
> = DO PRA AZ 

Da Equação (3.62) T=2m/i (1+ $) = 1 (6,3310) 

Da fórmula elementar To = 2r} = E (6, 2832) 


3.16 O Problema Isocrônico 


É interessante investigar se existe ou não uma curva de vínculo para a qual uma partícula 
oscilará sob a ação da gravidade isocronicamente, isto é, com um período que é independente 
da amplitude. 

Seja 0 o ângulo entre a horizontal e a tangente a curva de vínculo (Figura 3.10). A com- 
ponente da força gravitacional na direção do movimento é —mg sen 0. A equação diferencial 
do movimento ao longo da trajetória de vínculo (considerado sem atrito) é então 


mõ = —mgsen 0 


Contudo, se a equação acima representar movimento harmônico simples ao longo da curva, 
devemos ter 


mS=-KS 


Desse modo, uma curva de vínculo que satisfizer a equação 
S = Csen 0 


produzirá movimento harmônico simples. 
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mgsenO 
O x 


Figura 3.10: Forças envolvidas no caso isocrônico. 


Podemos encontrar «x e y em termos de 0 a partir das equações acima, 


a z = (cos 0) (C cos 0) 
Logo 
r= focos Odo = Cia + sen 20) (3.63) 
Do mesmo modo W E 
o (sen 0) (C cos 0) 
e então 
y = fe sen 8 cos 0d0 = s cos 20 (3.64) 


Equações (3.63) e (3.64) são as equações paramétricas de uma ciclóide. Então uma curva de 
vínculo na forma de uma ciclóide produzirá um movimento em que S varia harmonicamente 
com o tempo, e o período de oscilação será independente da amplitude. Como corolário, 
vemos que uma partícula que começa a mover-se a partir do repouso numa curva cicloidal 
lisa gasta o mesmo tempo para alcançar o fundo independentemente do ponto de onde 
começa. 

O físico e matemático Holandês Christian Huygens descobriu este fato quando tentava 
melhorar a precisão de pêndulos de relógio. Ele descobriu a teoria das evolutas e mostrou 
que a evoluta de uma ciclóide é também uma ciclóide. Dessa maneira, obrigando o pêndulo a 
descrever uma trajetória cicloidal o período será então, independente da amplitude. Apesar 
de simples, a idéia nunca teve muito uso prático. 


3.17 O Pêndulo Esférico 


Um problema clássico de movimento vinculado é o de uma partícula que é obrigada a mover- 
se em uma superfície esférica lisa, por exemplo, um pequeno corpo deslisando dentro de uma 
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bola esférica lisa. O caso talvez fique melhor ilustrado por um pequeno corpo pesado preso 
na extremidade de uma corda ou haste leve e inextensível que é livre para balançar em 
qualquer direção em torno de um ponto fixo, como na Figura 3.11. Este é o famoso pêndulo 
esférico. 


X 


Figura 3.11: O pêndulo esférico. 


Soluções aproximadas em coordenadas retangulares 


Existem duas forças atuando na partícula, ou seja, a força da gravidade para baixo, e a 
tensão S da corda ou haste vinculadora. A equação diferencial do movimento é 


me =mg +5 


Escolhendo o eixo z na vertical, as componentes retangulares da equação do movimento são 


mü = Sz 
mj = Sy 
mã = S, -mg 


Uma solução aproximada é prontamente conseguida para o caso onde o deslocamento a partir 
da posição de equilíbrio é muito pequeno. O valor da tensão é aproximadamente constante 
e igual a mg, e temos |z| & l, |y| « l, z = 0. Dessa maneira as componentes £ y de § são 
dadas aproximadamente por 
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que podem ser facilmente verificadas a partir da geometria da figura. As equações diferenciais 
do movimento em x e y se reduzem então a 


ë+žs=0 
E g 
iao 
es 


que são similares às equações do oscilador harmônico bidimensional tratado anteriormente 
na Seção 3.10. As soluções, como já vimos são 


q = Acos(wt + a) 


y = B cos(wt + p) 


onde w = A como no pêndulo simples. 
Dentro dos limites de nossa aproximação o movimento é tal que sua projeção no plano 
xy é uma elípse. Existem, claro, casos especiais em que a projeção se reduz a uma linha 


reta, ou a um círculo, dependendo das condições iniciais. 


Solução em coordenadas esféricas 


Usaremos coordenadas esféricas (veja a Figura 3.11) para um tratamento do pêndulo 
esférico mais preciso do que foi feito anteriormente. A tensão S tem apenas componente 
radial, mas o peso mg tem ambas, uma componente radial mgcos0 e uma transversal 
—mg sen 0. Desse modo a equação diferencial do movimento escrita em componentes esféricas 
fica 


ma, = F.=mgcos0-— S 
mas = Fə = -—mgsen0 
mas = F =0 


As três componentes da aceleração a,, ag e ay foram dadas no Capítulo 1, Seção 1.25. Já 
que o vínculo é 
r = | = constante 


podemos ignorar a componente radial da aceleração. As outras duas componentes se reduzem 
a 


ay = lÖ-— lọ? senð cos0 
ad = lg sen 8 + 2lġô cos 8 


Transpondo termos e fazendo os cancelamentos óbvios, as equações diferenciais em 0 e q 
tornam-se a g 
0 — 6º sen 8 cos 0 + send = 0 (3.65) 


1 d 
sen 6 dt 


(ósen20) = 0 (3.66) 
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A segunda equação implica que a quantidade entre parênteses é constante. Vamos chamá-la 
h. Na verdade, ela é o momentum angular (por unidade de massa) em torno do eixo vertical. 
Esta grandeza é constante porque não existe nenhum momento de força em torno deste eixo. 
Podemos, então, escrever 


h 


bj = —— 3.67 
9 sen? 0 l ) 
Levando o valor de é na Equação (3.65), obtemos uma equação na coordenada 0 apenas 


0 
cos 0 o 


ô + Z song — k? (3.68) 


sen 6 


É instrutível considerar neste ponto alguns casos especiais. Primeiro, se o ângulo q for 
constante, então q = 0 e logo h = 0. 
Como consequência, a Equação (3.68) se reduz a 


0 + sen =) 
que é, exatamente a equação diferencial do pêndulo simples. O movimento se faz no plano 
b = do = constante. 
— O segundo caso especial é aquele do pêndulo cônico; O = Oy = constante. Neste caso 
0 = 0 e 8 = 0, logo a Equação (3.68) se reduz a 


g = 2 C0800 + 
j sen 0, — h e 0 
ou 
hº = sen” ĝo sec ĝo (3.69) 


Usando o valor de h da Equação (3.69), encontramos a partir da Equação (3.67) 


4 = 5 sec ĝo (3.70) 
que é a condição de movimento cônico para o pêndulo. 

Esta última equação pode ser obtida também considerando-se as forças atuantes na 
partícula durante seu movimento circular conforme a Figura 3.12. A aceleração é cons- 
tante em módulo, isto é, pó? = (l sen 00) 6%, e aponta diretamente para o centro da trajetória 
circular. Dessa maneira, podemos escrever as componentes horizontal e vertical 


S sen 06 
S cos ĝo 


(ml sen 09) 62 
0 


que se reduz à Equação (3.70) depois de eliminarmos S. 

Vamos agora considerar o caso em que o movimento é quase cônico; isto é; o valor de 6 
permanece próximo do valor 09. Substituindo na Equação (3.68) h? por seu valor dado em 
(3.69), encontramos 


— 3.71 
cos ĝo sen? 6 ( 


.. 4 
js A [seno - sen soort) -0 
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Figura 3.12: O caso cônico do pêndulo esférico. 


E conveniente, a esta altura, introduzir a nova variável £ definida como 
E=0- 


A expressão entre parênteses na Equação (3.71) pode ser desenvolvida como uma série de 
potências em é de acordo com a fórmula conhecida 
2 


FCE = F0) + F'(0)E + PO: ds 


Depois de executar as operações indicadas encontramos que f(0) = 0 e f(0) = 3 cos 6o + 
sec ło. Como estamos interessados no caso de pequenos valores de é, desprezamos potências 
mais altas do que a primeira, e a Equação (3.71) pode ser escrita 


É+ Toe =ü (3.72) 


onde b = 3 cos fo + sec do. O movimento em é ou 0 é dado então por 


E = 0 — Oo = £o cos [o + s) (3.73) 


Então 0 oscila harmonicamente em torno do valor 0, com um período 


[1 l 
To= Sr =n 3.74 
D gb nz cos 0, + sec 0o) Mera 


O valor de é, dado pela Equação (3.67), não difere muito do valor dado pelo movimento 
puramente cônico do, indicando que é aumenta estacionariamente durante a oscilação de 
0 em torno de 06. A Figura 3.13 mostra a trajetória da partícula. Durante uma oscilação 
completa de 0 o valor do ângulo azimutal q aumenta de 


dı E Poli 
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Usando os valores de ġo e Tı dados acima, encontramos 


1 


dy = 27 (3 cos? ĝo + r 


Seja p o raio do círculo para 0 = ĝo, como na Figura 3.12. Então cos? 0, = 1 — A consequen- 


1 
temente ġı = 27 (4 — 392/12) 2. Então à, é ligeiramente maior que 7. Desenvolvendo-se em 
potência de p encontramos que o incremento Ag é dado por 


Pon) 
po me N ; É 
e 


Figura 3.13: Projeção no plano xy da trajetória de movimento do pêndulo esférico. 


No início desta seção provamos que a projeção da trajetória do pêndulo no plano xy era 
aproximadamente uma elípse para ângulo 0 pequeno. Podemos agora interpretar o resultado 
acima dizendo que o eixo maior da elípse não está estacionário, mas precessa na direção de 
& crescente. O eixo da elípse gira de Ad durante cada oscilação completa de 0. A Figura 
3.13 ilustra isso. 


Considerações de energia — Limites do movimento vertical 


Afim de relacionar a amplitude da oscilação vertical do pêndulo esférico com os parâmetros 
do problema, é vantajoso usar a equação de energia. Na nossa notação, a energia potencial 
é 

V = —mgl cos 0 
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Para encontrar a energia cinética usamos o fato de que as componentes da veloci- 
dade em coordenadas esféricas são (7, r0, rósen0). Como r = | = constante, temos 
T = mv” = im(20º + Pg? sen? 0). A equação de energia fica 


mM fia 12 
E = o (ô + q“ sen 0) — mgl cos 0 
Vamos resolver esta equação para 62. Para conseguir isto, usamos a relação, desenvolvida 
anteriormente, 9 = h/ sen? 0. Vamos colocar ainda cos0 = u. O resultado é 


: 2E 2g h? 
8? = | — = i 
e a T a. (3.75) 
0-0 
—0= (95) 


Figura 3.14: Ilustrando o movimento vertical do pêndulo esférico. 


As raízes da equação f(u) = 0 determinam os limites, ou os pontos de retorno da oscilação 
em 6, porque Ê se anula nestas raízes. O limite fica restrito para aqueles valores de 8 que 
tornam f(u) positivo. Dessa maneira a oscilação vertical fica entre dois círculos horizontais, 
veja a Figura 3.14. Se, em particular, as duas raízes reais que ficam entre +1 e —1 forem 
iguais, então o movimento fica confinado a um único círculo horizontal, isto é, temos o caso 
do pêndulo cônico. 


Exercícios de Treinamento 


3.1 Calculando o rotacional encontre quais das seguintes forças são conservativas: 


a) F=TIx+Jy 
b) F =Ty +z 
c) F = Ty -fr 
d) F = Try +yz + kzz 
e) F = yz +7zx + key 
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3.2 


3.3 


3.4 


3.9 


Determine o valor da constante c que tornam as seguintes forças conservativas: 
a) F = Try + Jcr? 

b) F =74 +7c4 +k 

Encontre a força associada a cada uma das seguintes funções energia potencial 
a) V = tyz 

b) V=k(2º+98 +z) 

c) V = kz®yf zI 

d) V 


— kelaz+By+y2) 


Encontre a função energia potencial para aquelas forças do problema (3.1) que forem 
conservativas. 


Uma partícula de massa m move-se no campo de força dado pela função energia 
potencial do exercício 3.3 (a). Se a partícula passa pela origem com velocidade vo, 
qual será sua velocidade se e quando passar pelo ponto (2,2,2)? 


Problemas 


3.6 


Considere as duas funções força 

a) F=Tx+Ty 

b) F=Ty-Jr 

Verifique que (a) é conservativa e (b) é não conservativa mostrando que a integral 
JF . dré independente da trajetória de integração para (a), mas não para (b), tomando 
dois caminhos para os quais a origem (0,0) é o ponto de partida e (1,1) é o ponto de 
chegada. Para uma trajetória tome a linha x = y. Para a outra, tome o eixo y até o 
ponto (0,1) e depois a linha y = 1 até o ponto (1,1). 


3.7 Mostre que a variação da gravidade com a altura pode ser explicada aproximadamente 


3.8 


pela seguinte função energia potencial: 


Z 
V = = 
mg (1-5) 


onde R é o raio da Terra. Encontre a força associada com esta função energia potencial. 
Deste resultado, encontre as componentes da equação diferencial do movimento de um 
projétil submetido a esta força. 


Um projétil é lançado da origem com velocidade inicial vọ com ângulo de elevação 0 
em relação a horizontal. Mostre que se desprezarmos a resistência do ar e se o solo for 
tomado como nível, o projétil atinge o solo a uma distância vá sen 28 da origem. Este é o 
alcance horizontal do projétil. Mostre também que o decréscimo no alcance horizontal 
é aproximadamente 4vgy sen 6 sen 26/39? no caso de resistência do ar proporcional à 
velocidade (caso linear). 
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3.9 Partículas de lama são lançadas da periferia de uma roda que rola sem deslizar. Se 
a roda tiver uma velocidade para a frente vo e se seu raio for b, mostrar que a maior 
altura acima do solo que a lama pode atingir é 


2 2 
v gb 
Rg A 
E 2g 5 208 
Em que ponto da roda girante irá a lama abandoná-la? 


3.10 Um canhão é assentado na parte mais baixa de um morro de declive constante q. 
Mostre que o alcance desta arma medida a partir da encosta do morro é 


2vê cos 0 sen(6 — d) 
g cos? À 


onde 0 é o ângulo de elevação do canhão. Mostre também que o valor máximo do 
alcance medido ao longo da encosta é 


E 
g(1 + sen q) 


3.11 Escreva na forma de componentes a equação diferencial de movimento de um projétil 
para o caso de resistência do ar proporcional ao quadrado da velocidade. Estas equações 
são separáveis? Mostre que a componente x da velocidade é 


& = toe PÉ 
onde S é a distância percorrida pelo projétil ao longo da trajetória. 


3.12 As condições iniciais para um oscilador isotrópico bi-dimensional são: æ(0) = 4, 
(0) = 0, y(0) = B w(0) = wC onde w é a frequência de oscilação. Mostre que o 
movimento se dá inteiramente dentro de um retângulo de dimensões 24 e 2y B2 + C°. 
Encontre a inclinação q do eixo da elipse em termos de 4, Be C. 


3.13 Uma partícula de massa unitária move-se no potencial do oscilador harmônico não 
isotrópico tri-dimensional 
V =r +4 +97 


Se a partícula passa pela origem com velocidade unitária na direção (1,1,1) no instante 
t = 0, determine x, y e z como funções do tempo. 


3.14 Um átomo esta situado em uma rede cristalina cúbica simples. Se a energia potencial 
de interação entre dois átomos quaisquer é da forma cr“ onde c e a são constantes 
e r é a distância entre dois átomos, mostre que a energia total de interação de um 
dado átomo com seus seis vizinhos mais próximos é aproximadamente o potencial do 
oscilador harmônico tri-dimensional 


V=A+B(2 +y +2) 
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onde A e B são constantes. (Nota: Considere que os seis átomos vizinhos estão fixos 
e localizados nos pontos (+d,0,0), (0,+d,0), (0,0,+d), e que o deslocamento (z,y,z) 
da posição de equilíbrio (0,0,0) do átomo é pequeno se comparado com d). Então 


= er" onde rı = (d — £)? + y? + 2? com expressão semelhante para rz, 73, 
EE , re. Use o teorema binomial para obter o resultado desejado. 


3.15 Um eletron move-se em um campo de força composto de um campo elétrico uniforme É 
e um campo magnético uniforme B que faz ângulo reto com E. Seja E= JE e B=kB. 
Considere a posição inicial do eletron na origem com velocidade inicial vo = Tvo na 
direção x. Encontre o movimento resultante da partícula. Mostre que a trajetória do 
movimento é uma ciclóide: 


x = asenwt + bt 
= c(1— coswt) 
2 0=0 


Movimento cicloidal de eletrons é utilizado no magnetron — um tubo eletrônico usado 
para produzir ondas de rádio de alta frequência. 


3.16 Uma partícula foi colocada na superfície externa de uma esfera lisa de raio b a uma 
distância b/2 acima do plano central. Determine em que ponto a partícula abandona 
a superfície, depois que ela começa a deslizar para baixo. 


3.17 Uma bolinha (conta) desliza em um fio liso dobrado na forma de um aro de raio b. Se 
o aro é vertical, e se a bolinha estava inicialmente em repouso, em um ponto nivelado 
com o centro do aro, encontre a velocidade da bolinha no ponto mais baixo do aro. 
Encontre também a reação do aro sobre a bolinha neste ponto. 


3.18 No problema acima determine o tempo que a bolinha gasta para atingir o ponto mais 
baixo do aro. 


3.19 Em uma experiência de laboratório um pêndulo simples é usado para se determinar o 
valor de g. Se a amplitude de oscilação do pêndulo é 30º, encontre o erro cometido ao 


se usar a fórmula elementar 
l 
T= mf 
g 


3.20 Um pêndulo esférico de 1 metro de comprimento está executando pequenas oscilações 
em torno de um ângulo cônico 06. Se 66 = 30°, encontre o período do movimento 
cônico, o período de oscilação de O em torno de 0o e o ângulo de precessão Ag. 


3.21 Prove que as duas raízes reais da equação f(u) = 0, Equação (3.75), que estão entre 
+1 e —1 são iguais no caso de pêndulo cônico. 


3.22 A corda de um pêndulo esférico de comprimento l está inicialmente fazendo um ângulo 
de 90º com a vertical. Dá-se à massa do pêndulo uma velocidade inicial vo perpendicular 
à corda. Se vg = gl/2 encontre o ponto mais baixo atingido pela massa do pêndulo 
durante seu movimento. (Sugestão: A partir das condições iniciais u = 0 é uma raiz 


da equação f(u) = 0). 
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Respostas de problemas selecionados de números ímpa- 
res 


3.1 (a), (b), e (e) são conservativas. 
33 (a) F = —(Tyz +7zz + key) 

(c) F = kr z (Tar! +78y 7! + kyz") 
N T 


3.7 F= -mg + ZZ 
3.11 m = “cas, my = —cys, mz = —mg — ciS 
3.13 x = 62 sen(V2t), y = 2473 sen(vBt), z = 543 sen(V18t) 
3.17 v = v29b, R = 3mg 
AA o 
3.19 T = 4 /} K(15°) 


Capítulo 4 


Dinâmica de um Sistema de Muitas 
Partículas 


Ao estudarmos um sistema ou conjunto de muitas partículas livres, é importante que fique- 
mos interessados nas características gerais do movimento de tal sistema. 


4.1 Centro de Massa e Momentum Linear 


Consideremos um sistema constituído de n partículas de massas M1, ma, ... Mn, cujos vetores 
posição são, respectivamente r4,I>,...r,. Definimos o centro de massa do sistema como o 
ponto cujo vetor posição Tem (Figura 4.1) é dado por: 


E Mırı + mol» Eos Mal = 5 Mil; 


Tem A 


(4.1) 


Mı FM +... Mp m 


onde m = >) nm; é a massa total do sistema. A definição acima é certamente equivalente às 
três equações: 


O miti o Xo miyi X mizi 
E IS Vem + komi E 
m m m 


Tem 


Definimos o momentum linear P do sistema como o vetor soma, dos momenta individuais 
das partículas 
Da Equação (4.1), por diferenciação em relação ao tempo, temos 
P= Y mvi = MV em (4.3) 


o que significa que o momentum linear do sistema de partículas é igual ao produto da 
velocidade do centro de massa pela massa total do sistema. 

Suponha agora que existam forças externas F, Eb, as E; ata E, agindo nas respectivas 
partículas. Em adição suponha, também, que existam forças internas de interação entre (duas 
a duas) quaisquer partículas do sistema. Vamos denotar estas forças por E, significando as 
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y 


Z 


Figura 4.1: Centro de massa de um sistema de partículas. 


forças exercidas na partícula i pela partícula j. A equação do movimento da partícula i é 
então 
F, + 5 Fa z mit; = P; (4.4) 
j=I;jżi 


onde F; é a força externa total agindo em i. O segundo termo da equação acima corresponde 
à soma de todas as forças internas exercidas na partícula à pelas outras partículas do sistema. 
Somando a expressão (4.4) para todas as partículas, temos 


TEAS | 5 Fy=) 6, (4.5) 


Na somatória dupla, acima, para cada força F;; existirá também a força F;; e elas são iguais 
e opostas 


por consequência da terceira lei de Newton. Desta maneira, as forças internas cancelam-se 


duas a duas e a somatória dupla vale zero. Podemos reescrever a Equação (4.5) do seguinte 
modo 


SF; = DE = P = Mãem (4.7) 


em outras palavras a aceleração do centro de massa de um sistema de partículas é a mesma 
que uma única partícula de massa igual à massa de todo o sistema teria se submetida a uma 
única força igual à soma de todas as forças externas. 

Considere, por exemplo, um conjunto de partículas movendo-se em um campo gravita- 
cional uniforme. Uma vez que F, = mi;g, para cada partícula, teremos 


DS FE=D mg=mg 
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lembre-se que g é constante se o campo for uniforme. Dai podemos tirar que 
dan =Ë (4.8) 


Desta maneira, o centro de massa de uma granada que explodiu no ar continua seu caminho 
parabólico que ela possuia antes da explosão. 

Em um caso especial em que não há nenhuma força externa agindo no sistema (ou se 
SE; = 0), temos que äm = 0 e Von = constante. De modo que o momentum linear do 
sistema permanece constante 


5p = P = MV em = constante (4.9) 


Este é o princípio de conservação do momentum linear. Na mecânica Newtoniana a constância 
do momentum linear de um sistema isolado está diretamente ligada, e de fato é uma con- 
sequência da terceira lei de Newton. Mas, mesmo nos casos em que as forças entre as 
partículas não obedecem diretamente à lei da ação e reação, como no caso da força magnética 
entre duas partículas móveis carregadas, o princípio de conservação do momentum linear 
ainda é mantido, quando se leva em conta o momentum linear das partículas e do campo 
eletromagnético!. 


4.2 Momentum Angular do Sistema 


Como foi visto no Capítulo 3 (Equação 3.4), o momentum angular de uma partícula é definido 
pelo produto vetorial rx mv. O momentum angular L de um sistema de partículas é definido 
como a soma vetorial dos momenta angular individuais das partículas 


Calculemos a derivada temporal do momentum angular. 


dL n n 
d a i=1 
aqui a primeira somatória é nula (V; x V; = 0) e na segunda somatória o produto md; 


corresponde à força total atuante na partícula i. Podemos então escrever 


Toile (r 5 fo) = Bux 45) 5 Ren; (4.11) 
i=1 


j=jfi i=1 j=1;jżi 


Na expressão acima é usada a mesma nomenclatura da Seção 4.1, onde E, corresponde à 
resultante externa de forças sobre a partícula à e F;; é a força interna da partícula j sobre a 
partícula à. A somatória dupla da expressão (4.11) é constituida de pares 


(E x Fy) + (E; x Fj) (4.12) 


!Veja por exemplo, W. T. Scott, The Physics of Eletricity and Magnetism, 22 ed. John Wiley & Sons, 
Inc., New York, 1966. 
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Se representarmos o vetor deslocamento da partícula j em relação à partícula à por Tij, 
podemos ver pela Figura 4.2 que 


O 


Figura 4.2: Definição do vetor Ti. 


Como Fy = P podemos então reduzir a expressão (4.12) a 


que é nula, caso as forças internas sejam centrais, isto é, caso elas ajam ao longo da linha 
que une as partículas. Desta forma, a somatória dupla, de (4.11) é nula. De acordo com o 
que foi definido na Seção 1.12, o produto vetorial r; x F; corresponde ao momento da força 
externa F;. A soma Darn X Ë, é portanto, o momento total de todas as forças externas que 
agem sobre o sistema. Se denotarmos o momento externo total por N, podemos escrever a 
Equação (4.11) do seguinte modo 
TN 
dt 
que significa: a taxa de variação do momentum angular do sistema com o tempo é igual ao 
momento total das forças externas agindo sobre o sistema. 
Se o sistema é isolado, então N = 0 e o momentum angular permanece constante tanto 
em módulo quanto em direção e sentido, 


(4.15) 


L= 5r X m;v; = constante (4.16) 


Este fato é estabelecido como Princípio de Conservação do Momentum Angular. 
Do mesmo modo que o momentum linear de um sistema isolado permanece constante, o 
momentum angular de um sistema isolado também permanece constante, mesmo no caso 
de um sistema de cargas móveis onde temos que considerar o momentum angular do campo 
eletromagnético”. 


2Ver nota nº 1. 
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4.3 Energia Cinética do Sistema de Partículas 


A energia cinética total T do sistema de partículas é dada pela soma das energias individuais 


mit; . 7:) (4.17) 


Figura 4.3: Definição de Tr). 


Como mostra a Figura 4.3, podemos expressar cada vetor posição T; na forma 
È = Em + E; (4.18) 


onde T; é a posição da partícula į em relação ao centro de massa. Tomando a derivada em 
relação ao tempo, teremos 


Vi = Vem + Vi (4.19) 


aqui, Vem é a velocidade do centro de massa e V; a velocidade da partícula į em relação ao 
centro de massa. A expressão para T pode então ser escrita 


S imin +7). (Vem + Vi) e a ES Vem . Vi +55 mT 


1 
mas da Equação (4.18), temos 


analogamente, temos 


T 


e a expressão para a energia cinética do sistema fica 


1 1 
T= 5 Nem DD gmi? (4.20) 
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o que significa que a energia cinética total do sistema é dada pela soma da energia cinética 
de translação do centro de massa (1º termo) mais a energia cinética do movimento das 
partículas com relação ao centro de massa (último termo). Esta separação de energias 
cinéticas é conveniente, por exemplo, em Física molecular. Desta maneira, para um sistema 
molecular, a energia cinética total consiste na soma da energia cinética de translação de 
todas as moléculas, mais a energia de vibração e rotação de cada molécula. 


4.4 Movimento de Dois Corpos que Interagem — Massa 
Reduzida 


Consideremos o movimento de um sistema constituido de dois corpos (tomados como partículas) 
que interagem através de força central. 

Vamos tomar um sistema isolado e portanto com centro de massa com velocidade cons- 
tante. Por simplicidade tomaremos o centro de massa como origem. Temos então 


mari + moro = 0 (4.21) 


onde, como mostra a Figura 4.4, os vetores Tj e ro representam a posição das partículas de 
massa my e mo, respectivamente, em relação ao centro de massa. Agora, se R é o vetor 
posição da partícula 1 relativo à partícula 2, então 


Figura 4.4: O problema dos dois corpos. 


R=h-b=5 (1 $ m) (4.22) 


a última passagem sai da Equação (4.21). 
A equação diferencial do movimento da partícula 1 relativa ao centro de massa é 


(4.23) 
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onde f(R) é o módulo da força mútua entre as duas partículas. Usando a Equação (4.22), 
podemos escrever 


PŘ R 
enn ae 4.24 
na = BE (4.24) 
onde né 
1 2 
Ee ta re 4.25 
ES (4.25) 


A quantidade u é conhecida por massa reduzida. A nova equação do movimento (4.24) dá 
o movimento da partícula 1 em relação à partícula 2 (tomada como fixa). Esta equação é 
precisamente a mesma equação para o movimento de uma partícula de massa u movendo-se 
em um campo central de força dada por f(R). Deste modo, o fato de mo» estar em movimento 
relativo ao centro de massa está automaticamente computado quando se substitui mı por pu. 
Se os corpos têm massas iguais, então u = m/2. Por outro lado, se ma é muito maior que 
ma, então u é praticamente igual a m. 
Para dois corpos que se atraem por gravitação nós temos 


Gmm 
F(R) = Se (4.26) 
Neste caso, a equação do movimento é 
= Gmims R 
a E 4.27 


Esta é a mesma equação para o movimento de uma única partícula sujeita a uma força 
central regida por uma força dependente do inverso do quadrado da distância. Uma vez que 
a escolha dos subscritos (de mı ou ms) é arbitrária, devemos concluir que o movimento de 
cada partícula em relação à outra é uma curva cônica com esta fixa no foco. 

Deste modo, tomando o sistema Terra-Lua como um sistema isolado, a Lua descreve uma 
elipse em relação à Terra (com a Terra no foco) ou, se o referencial for a Lua, temos a Terra 
descrevendo uma elipse em torno da Lua (com a Lua no foco). 


4.5  Colisões 


Sempre que dois corpos colidem, a força que cada um exerce no outro durante o contato é 
uma força interna, caso os dois sejam considerados um sistema. O momentum linear total 
fica portanto inalterado. Podemos então escrever 


Bı + B2 = Pi + p2 (4.28) 


ou equivalentemente, 
> > > — 
MV1 + MV2 = Mm + mavs (4.29) 


os subscritos 1 e 2 referem-se às duas partículas e as linhas caracterizam os momenta e 
velocidades das partículas após a colisão. As equações acima são completamente gerais. 
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Elas se aplicam a quaisquer dois corpos, independentemente de seus formatos, rigidez ou 
qualquer outra coisa. 
Com relação ao balanço energético podemos escrever 


2 2 12 12 
Pı | E | Po HQ (4.30) 

2M: 2Moə 2my 2mos 

né 1 1 1 1 
Mai + 5 Mata = zmaw'i | zmz’ -Q (4.31) 


Aqui a quantidade Q é introduzida para indicar a energia total perdida, ou ganha, que ocorre 
resultante da colisão. 

No caso de uma colisão perfeitamente elástica, não existe mudança na energia cinética 
total de modo que Q = 0. Se existir uma perda de energia, então Q é positivo. Esta colisão 
é denominada endoérgica. Pode também ocorrer que haja aumento da energia cinética total. 
Isto pode ocorrer por exemplo se houver alguma explosão no ponto de contato dos dois 
corpos. Neste caso, Q é negativo e a colisão é denominada exoérgica. 

O estudo das colisões é de particular importância em física atômica e nuclear. Neste caso, 
os corpos envolvidos são átomos, núcleos ou partículas elementares como eletrons, protons 
ou outras. 


Colisão Direta 


Consideremos o caso especial de uma colisão frontal de dois corpos ou partículas, no qual 

o movimento ocorre inteiramente em uma linha reta como mostra a Figura 4.5. Neste caso, 

a equação de conservação do momentum linear (4.29) pode ser escrita sem a notação vetorial 
como 

Mv + MV2 = Mv] + mova (4.32) 


o sentido do movimento das partículas é dado pelo sinal de v. 


m] mz 
o~) a 
mı Mo 
AC X 
Z€ > SA Z A 


Figura 4.5: Colisão frontal de duas partículas. 


A fim de calcular as velocidades após a colisão, devemos usar a Equação (4.32), junta- 
mente com a equação do balanço energético (4.31). 
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E conveniente, em problemas deste tipo, introduzir-se o parâmetro e conhecido por co- 
eficiente de restituição. Esta quantidade é definida como a razão entre a velocidade de 
separação V” e a velocidade de aproximação V. Em nossa notação é deve ser escrito como 


on] VM 


= = 4.33 
[va = v| V ( 


O valor numérico de e depende primariamente da composição e constituição física dos dois 
corpos. É fácil verificar que em uma colisão perfeitamente elástica ce = 1. Para tanto, 
devemos fazer Q = 0 em (4.31) e resolver o sistema formado pelas Equações (4.31) e (4.32). 

Em uma colisão perfeitamente inelástica, os dois corpos permanecem unidos após a colisão 
de modo que e = 0. Para muitas situações o valor de € situa-se entre os dois extremos (entre 
0 e 1). Para a colisão entre duas bolas de bilhar e vale aproximadamente 0,95. O valor desse 
coeficiente depende também da velocidade com que os corpos se aproximam. 

Podemos calcular os valores das velocidades finais com o uso da Equação (4.32), junta- 
mente com a definição do coeficiente de restituição, Equação (4.33). O resultado é 


is (mı — emoJvy + (mo + emo)va 
i May T Ma 
e (my + emiw + (ma — em w2 (4.34) 
Ma T Ma 


No caso inelástico (e = 0) verifica-se que vj = v5, isto é, não há recuo relativo. Por outro 
lado, no caso especial em que os dois corpos possuem massas iguais, my = ma, e a colisão é 
perfeitamente elástica, e = 1, obtemos 


os dois corpos trocam suas velocidades como resultado da colisão. 
No caso geral de colisões frontais não elásticas, é facil verificar que a perda de energia Q 
é relacionada com o coeficiente de restituição pela equação 


Q= su =e) 


onde u é a massa reduzida e v a velocidade relativa antes da colisão. A dedução desta 
expressão é deixada como exercício. 


4.6 Colisão Oblíqua e Espalhamento. Comparação das 
Coordenadas do Laboratório e do Centro de Massa 
Vamos voltar nossa atenção para um caso de colisão mais geral do que aquele em que o 


movimento está confinado a uma linha reta. Aqui a notação vetorial para as equações do 
momentum (4.28) e (4.29) tem que ser mantida. 
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Vamos estudar o caso especial de uma partícula de massa mı com velocidade inicial 
vi (partícula incidente) que colide com a partícula de massa ma inicialmente em repouso 
(partícula alvo). Este é o problema típico que encontramos em física nuclear. As equações 
do momentum são 


Bı = Pi + P2 (4.35) 
mui = miyi + mov (4.36) 


O balanço de energia é dado por 


H H 4.37 
2M: ma Ma Q ( ) 
Ei 1 1 1 

zma = zma’ | zma’ | Q (4.38) 


Aqui, como antes, as linhas indicam velocidades e momenta finais. A quantidade Q é de 
fundamental importância em física atômica e nuclear uma vez que ela representa a energia 
emitida ou perdida em uma colisão atômica ou nuclear. Em muitos casos a partícula alvo é 
“quebrada” (fissão) ou se modifica por influência da colisão. Em tais casos as partículas, após 
a colisão, são diferentes das partículas antes da colisão. Em tais casos deve-se, logicamente, 
assinalar de modo diferente as massas após a colisão com ms e m4, por exemplo. 

Mas, qualquer que seja o caso, a lei de conservação do momentum tem que ser verificada. 
De acordo com a teoria da relatividade, no entanto, a massa da partícula varia com a veloci- 
dade de um modo que é descrito por esta teoria, e tal fato deve ser tomado em consideração 
quando tratamos colisões envolvendo velocidades próximas à da luz. No entanto, podemos 
adiantar que a conservação do momentum expressa pela Equação (4.28) é relativisticamente 
correta, considerando que a massa seja uma função da velocidade. 


Coordenadas do Centro de Massa 


Cálculos teóricos em física nuclear são feitos comumente em termos de quantidades rela- 
cionadas a um sistema de coordenadas no qual o centro de massa das partículas que colidem 
esteja parado. Por outro lado, as observações experimentais do espalhamento das partículas 
são feitas em termos das coordenadas de laboratório. Há portanto o interesse de se fazer 
considerações breves à respeito do problema de conversão de um sistema de coordenadas 
para o outro. 

Os vetores velocidade no sistema do laboratório e no sistema do centro de massa estão 
representados na Figura 4.6. Nessa Figura q, é o ângulo de deflexão da partícula incidente, 
depois que ela atingiu o alvo e q» é o ângulo que a direção do movimento do alvo faz com a 
direção original do movimento da partícula incidente. Tanto q, quanto 2 são medidos no 
sistema do laboratório. No sistema do centro de massa, uma vez que o centro de massa deve 
permanecer sobre a linha que une as duas partículas, ambas as partículas aproximam-se do 
centro de massa, colidem e afastam-se do mesmo, em sentidos opostos. O ângulo 0 é o ângulo 
de deflexão da partícula incidente em relação à direção original de movimento, medido no 
sistema do centro de massa. 
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SISTEMA 
CENTRO de MASSA 


Mı 
DN 


E PEDR E EE: PRETA 


SISTEMA 
LABORATORIO 


Figura 4.6: Comparação das coordenadas do Laboratório e do Centro de Massa. 


Da definição de centro de massa, vemos que o momentum linear do conjunto no referencial 
do centro de massa é nulo tanto antes quanto depois da colisão. Desta forma, podemos 
escrever o 

Pı tp = 0 (4.39) 
Pı +P = 0 (4.40) 


As barras são usadas para indicar que a quantidade em questão se refere ao sistema do centro 
de massa. O balanço energético corresponde a 


2) =3 2 72 
a L (4.41) 


2my4 2mo 2m4 2mo 
Podemos eliminar p, e p'o da última equação recorrendo às Equações (4.39) e (4.40). Desta 
forma, a Equação (4.41) pode ser escrita como 


—2 2 
Pi _ Pi 
— = — 4.42 
DR Q (4.42) 
As expressões (4.39) e (4.40) podem ser escritas como 
mA ++ MV2 = 0 (4.43) 
mv ++ mov!» = 0 (4.44) 
A velocidade do centro de massa é 
na (4.45) 
ma + Mg 
então temos E 
5 mMə2V1 
Vi = Vi — Vem = —— 4.46 
V= V -y PE ( ) 


A relação entre os vetores velocidade Vem, Vi e V'ı é mostrada na Figura 4.7. Dela, vemos 
que 

visendi = Usen (4.47) 
U, cos Ô + vem (4.48) 


v| COS Q1 
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Dividindo (4.47) por (4.48) encontramos uma equação que conecta os ângulos de espalha- 
mento nos dois sistemas 


sen 0 
t anoo 4.49 
891 y + cos 0 ( 
onde y é um parâmetro numérico dado por 
Vem Mv 
pe Emei 4.50 
E v v (mı + mə) (ea 


a última passagem vem da Equação (4.45). 


Figura 4.7: Relação entre os vetores velocidade no sistema Laboratório e no sistema Centro 
de Massa. 


Podemos agora calcular o valor de v} em termos da velocidade da partícula incidente a 
partir da Equação (4.42). Esta equação nos dá as informações necessárias para encontrarmos 
y e então determinarmos a relação entre os ângulos de espalhamento. Por exemplo, no caso 
de uma colisão perfeitamente elástica, Q = 0, encontramos — de (4.42) — que D, = Dj ou 
vi = V| e este resultado junto com a Equação (4.46) nos leva ao resultado 


ma 


y= (4.51) 


Ma 


para a colisão elástica. 

Dois casos especiais para as colisões elásticas merecem ser considerados. Primeiro, se a 
massa mə da partícula alvo é muito maior do que a massa m; da partícula incidente. Neste 
caso y é muito pequeno e tg dj = tg0 ou q = 0. Isto significa que o ângulo de espalhamento 
visto no referencial do laboratório ou no referencial do centro de massa são praticamente 
iguais. 

O segundo caso especial é o de massas iguais (mı = ma). Neste caso y = 1 e a relação 


de espalhamento se reduz a 
sen 0 


1+cos6 | 
Qı = 2 
que significa que o ângulo de deflexão medido no referencial do laboratório é exatamente a 
metade daquele medido no referencial do centro de massa. Além disto, uma vez que o ângulo 


tg dy = 
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de deflexão da partícula alvo é 7 — O no referencial do centro de massa, como mostra a Figura 
(4.6), temos que o mesmo ângulo medido no referencial do laboratório é (m — 0)/2. Deste 
modo as duas partículas deixam o ponto de impacto em direções perpendiculares quando 
observadas no sistema do laboratório. 
No caso geral de colisões não elásticas, deixa-se como um problema a demonstração de 
que y pode ser expresso como 
ma 


Ma 


1-8 (1 | a (4.52) 


T Mg 


em que T é a energia cinética da partícula incidente, medida no referencial do laboratório. 
4.7 Impulso em Colisões 
Forças de duração temporal extremamente curta, tais como as exercidas por corpos que 


colidem, são denominadas forças impulsivas. Se concentramos nossa atenção em um dos 
corpos, a equação diferencial do movimento como sabemos é 


d(mv) >= 
=P 4.53 
Er (4.53) 
ou na forma diferencial . 
d(mv) = Fdt (4.54) 


integrando de t = tı até t = to que é o intervalo de tempo durante o qual a força atuou, 
teremos j 
pa 
Amy) = Fdt (4.55) 
ty 
A integral temporal da força é denominada impulso, e é denotada pelo símbolo p. Desta 


forma, a equação acima pode ainda ser escrita como 
A(myv) = P (4.56) 


Em palavras temos: A variação no momentum linear de um corpo sujeito a uma força 
impulsiva é igual ao impulso desta força. 

Podemos imaginar um impulso ideal produzido por uma força que tende para infinito 
com duração que tende para zero de modo que a integral / Fdt seja finita. Tal impulso ideal 
produziria uma mudança instantânea do momentum e da velocidade do corpo sem produzir 
qualquer deslocamento. 


Relação entre Impulso e o Coeficiente de Restituição 


Apliquemos o conceito de impulso ao caso de uma colisão frontal de dois corpos esféricos 
(tratados na Seção 4.5). Iremos dividir o impulso em duas partes denominadas impulsos 
de compressão fÑ. e impulso de restituição P,. Vamos considerar apenas as componentes ao 
longo da direção que une os centros das esferas. Podemos escrever, para fe 


Mv — Mv, = Pe (4.57) 


Mavo — MavVa =p. (4.58) 
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onde vo é a velocidade comum das duas esferas no instante em que a velocidade relativa é 
zero. Similarmente, para a restituição temos 
$ A 
Mvi mv = dr (4.59) 
1 A 
mov, — MW = —Pr (4.60) 


eliminando vo das Equações (4.57) e (4.58) e também das Equações (4.59) e (4.60) obtemos 
o par de equações 


MımMmo(vz = v) = Delma T mə) 


Mım(v’ı = v'a) = Pr(mı T ma) 


Dividindo a segunda equação pela primeira temos 


/ 1 É 
vVo-rvi o Pr 
Vi — Va Pc 


(4.61) 
mas o termo à esquerda é justamente a definição do coeficiente de restituição de modo que 


Dr 

Pc 
O coeficiente de restituição é portanto igual a razão entre o impulso de restituição e o impulso 
de compressão. 


4.8 Movimento de um Corpo com Massa Variável — 
Movimento de um Foguete 


No caso de um corpo cuja massa varia com o tempo, é necessário um certo cuidado ao lidar 
com as equações diferenciais do movimento, o conceito de impulso é bastante útil neste tipo 
de problema. 

Considere o caso geral do movimento de um corpo com massa variável. Seja Font a força 
externa que age no corpo em um dado instante e seja Am a variação de massa ocorrida no 
corpo em um curto intervalo de tempo At. Então, FenAtéo impulso devido à força externa 
e então temos , 

Fen At = (Protal) t+ At = (Protal)t 
para a variação total no momentum de um sistema naquele intervalo At. Se denotarmos 
por v a velocidade do corpo e por u a velocidade de Am relativa ao corpo, podemos então 
escrever 
Eert At = (m + Am) (7 + A?) — [my + Am(7 + T) 


que reduz a 7 

Fen Ãt = MAY + AMAY — Am 
ou, dividindo por At 
AV üAm 
At At 


Fent = (m + Am) 
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e no caso limite, quando At tende para zero, temos 
Feu = my — Um (4.63) 


Aqui a força externa pode representar a força da gravidade ou de resistência do ar ou outra 
qualquer. No caso dos foguetes o termo um representa a propulsão. 

Vamos aplicar a Equação acima a dois casos especiais. Primeiro, suponha que um corpo 
esteja se movendo em meio à neblina de modo que colete massa à medida que se move. Neste 


caso a velocidade inicial da matéria coletada é zero e u = —v logo temos 
= dm? 
E = mi + dh = dr) (4.64) 


para a equação do movimento: ela se aplica somente se a velocidade inicial da massa que está 
sendo recolhida é zero. Em outro caso, a expressão da Equação (4.63) tem que ser usada. 
No segundo caso especial, consideremos o movimento de um foguete. Para tal m é 
negativo por causa da perda de massa na forma de combustível ejetado. Assim um é oposto 
à direção de wu (velocidade relativa do combustível ejetado). Para simplificar, podemos 
resolver a equação do movimento para o caso em que a força externa é zero. Temos assim 


my = mn (4.65) 


separando as variáveis e integrando teremos 


jeja 
f w= f E 


Mo 
v = vo + u ln — 
m 


considerando u constante temos 


aqui mo é a massa inicial do foguete (carga + combustível não queimado), m é a massa do 
mesmo sistema no instante t e u é a velocidade do combustível ejetado relativa ao foguete. 
Devido à função logarítimica, é necessário que se tenha uma considerável razão entre a 
massa do combustível e a carcaça para que se alcance as altas velocidades necessárias para 
a colocação de um satélite em órbita. 


Exercícios de Treinamento 


4.1 Um sistema é constituído de três partículas, cada uma de massa unitária, com posições 
e velocidades instantâneas dadas por 


RetAgaE es 
b=T4E = 
=k 73 =T +7 +k 


Calcule as seguintes grandezas 
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a) A posição instantânea do centro de massa. 
b) À velocidade do centro de massa. 
c) O momentum linear do sistema. 
d) O momentum angular do sistema em relação à origem. 
e) À energia cinética. 

4.2 Um revolver de massa m dispara uma bala de massa ym onde y é uma pequena 
fração. Se vo é a velocidade da bala quando ela abandona o cano do revólver, qual é a 
velocidade de recuo do revólver? 

Problemas 

4.3 Um bloco de madeira de massa m repousa sobre uma superfície horizontal. Uma 
bala de massa ym é disparada horizontalmente com velocidade vo contra o bloco, 
alojando-se no mesmo. Que fração da energia cinética inicial da bala é perdida no 
bloco, imediatamente após o impacto? Se u é o coeficiente de atrito entre o bloco e o 
chão, que distância ele percorrerá até parar? 

4.4 Uma granada é lançada com velocidade inicial vo segundo um ângulo de 45º. Na parte 
mais alta da trajetória ela explode, dividindo-se em dois fragmentos de mesma massa, 
um dos quais move-se verticalmente para baixo com velocidade inicial V2vo. Qual a 
direção e velocidade do outro fragmento após a explosão? 

4.5 Três partículas de mesma massa movem-se sobre uma mesma linha reta. Inicialmente 
elas passam pelas posições —1; 0 e +1 com velocidades 4vo; 2vo e vo. Encontre as 
velocidades finais das partículas, considerando que todos os choques são perfeitamente 
elásticos. 

4.6 Uma bola cai de uma altura h chocando-se com um solo horizontal. Se o coeficiente 
de restituição é e, mostre que a distância vertical total percorrida pela bola até que os 
repiques cessem é 

h(1+ e) 
1— e 
e encontre o tempo total gasto até que a bola pare. 

4.7 Considere a Terra e a Lua como um sistema isolado e mostre que cada uma descreve 
uma elipse em relação ao centro de massa do sistema. 

4.8 Mostre que a energia cinética de um sistema de duas partículas é 


1 
z Vem es pu? 


onde m = my + ma, v é à velocidade relativa e u é a massa reduzida. 
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4.9 Se dois corpos colidem frontalmente, mostre que a perda de energia cinética é igual a 
uvZ(1 — e2)/2 onde u = massa reduzida; v = velocidade relativa antes do impacto e € 
= coeficiente de restituição. 


410 Uma partícula móvel de massa mı colide elasticamente com um alvo de massa ms 
inicialmente em repouso. Se a colisão for frontal, mostre que a partícula incidente perde 
uma fração 4u/m de sua energia cinética inicial (u = massa reduzida; m = my + mo). 


4.11 Mostre que o momentum angular de um sistema de duas partículas é 
Iom X MV em +R x uv 


(m = my +ma; u = massa reduzida; R = vetor posição relativa; v = velocidade relativa 
das duas partículas). 


4.12 Um próton de massa mp, com velocidade inicial v colide com um átomo de hélio, massa 
4mp, inicialmente em repouso. Se o próton deixa o ponto de impacto segundo um 
ângulo de 45º em relação à direção original de seu movimento, encontre as velocidades 
finais de cada uma das partículas. Considere a colisão perfeitamente elástica. 


4.13 Resolva o problema acima no caso de uma colisão inelástica em que Q = 1/4 da energia 
cinética inicial do próton. 


4.14 Relativo ao problema 4.12, encontre o ângulo de espalhamento do próton no referencial 
do centro de massa. 


4.15 Relativo ao problema 4.13, encontre o ângulo de espalhamento do proton no referencial 
do centro de massa. 


4.16 Uma partícula de massa m com momentum inicial pı, colide com uma partícula de 
igual massa em repouso. Se o módulo dos momenta finais das duas partículas são p!, 
e p', respectivamente, mostre que a perda de energia na colisão é dada por 


£ £ 
Q = P2 cosy 
m 


onde 1) é o ângulo entre as trajetórias das duas partículas após a colisão. 


4.17 Encontre a equação do movimento de um foguete disparado verticalmente para cima 
considerando g constante. Encontre a razão entre a massa de combustível e a da carcaça 
para que o foguete alcance uma velocidade final igual a velocidade de escape da Terra 
(ve), considerando-se que a velocidade de exaustão do gás é Kv. (K = constante) e 
a taxa de exaustão é m. Determine o valor numérico da relação combustível-carcaça 
para o caso em que K = 1/4, e m é igual a 1% da massa inicial de combustível por 
segundo. 


4.18 Monte a equação diferencial do movimento de uma gota de chuva que cai em meio à 
neblina; recolhendo mais água enquanto cai. Considere que a gota permanece esférica 
e que a taxa de recolhimento de água é proporcional à área da seção transversal da 
gota multiplicada pela velocidade da mesma. Mostre que, se a gota parte do repouso 
com tamanho inicial infinitesimal então a aceleração é constante e igual a 9/7. 
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4.19 Uma corrente de massa uniformemente distribuída e de comprimento a repousa inici- 
almente sobre uma mesa. Uma parte b de seu comprimento pende de uma das bordas 
da mesa enquanto a outra parte, de comprimento a — b repousa sobre a tábua da mesa. 
Se a corrente é abandonada, mostre que a velocidade da mesma, quando o último elo 
abandona a mesa, é vV2g(a? — b?) /3a?. 


4.20 Um foguete, ao viajar na atmosfera, experimenta uma força de resistência linear do 
ar — KV. Encontre a equação diferencial do movimento quando todas as outras forças 
externas são desprezíveis. Integre a equação e mostre que se o foguete parte do repouso, 
a velocidade final é dada por v = Va[1 — (m/mo)!º] onde V é a velocidade relativa 
de exaustão dos gases, a = |m/K| = constante, mo é a massa inicial do foguete mais 
combustível e m é a massa final do foguete. 


Respostas de Problemas Selecionados 


41 (a) (Mr +7+3k)/3 
(c) 3] +k p 
(d) —37 + 3k 


2442 
vY 


4.3 Fração perdida = (1+79) 1, s = Igu 


4.5 Vo; 200; 4vo 


4.12 Proton: v's = v'y = 0, 66v 
Hélio: v's = 0, 09v, v'y = 0, 16v 


4.14 x 55º 
4.18 ž = g — 32?°/z 


4.20 mù + Vm+ Kv = 0 
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Capítulo 5 


Mecânica dos Corpos Rígidos 
Movimento no Plano 


Um corpo rígido pode ser considerado como um sistema de partículas cujas posições relativas 
são fixas, ou, em outras palavras, em que a distância entre duas partículas quaisquer é cons- 
tante. Esta é a definição de um corpo rígido ideal. Primeiramente, como já foi mencionado 
na definição de uma partícula, não existem verdadeiras partículas na natureza. Em segundo 
lugar, os corpos reais não são estritamente rígidos; eles podem ser deformados (esticados, 
comprimidos ou curvados) quando submetidos a forças externas. No entanto, por enquanto, 
iremos desprezar tais deformações. 


5.1 Centro de Massa de um Corpo Rígido 


Já definimos na Seção 4.1 o centro de massa de um sistema de partículas como sendo o ponto 
(Lem, Yem, Zemi) tal que 


= X timi = X yim; = X zimi 
E Dn E Dnm’ nie Dmi 

Para um corpo rígido extenso, podemos substituir a somatória por uma integral sobre todo 
o volume do corpo de modo que 


qo — herd ae Jo pydv 2. = hpzdv 
T fpdo” E Jopdo” = fpd 


onde p é a densidade e dv o elemento de volume. Se o corpo rígido tem a forma de uma 
casca fina, as equações para o centro de massa passam a ser 


f. pxds I. pyds |, pzds 
Tem => 77, em > TR 1) Zem > 17 
J, pds dem f pds J, pds 
onde ds é o elemento de área e p, a densidade superficial; a integral se estende sobre toda a 
superfície do corpo. 
Analogamente, se o corpo tem a forma de um fio fino, temos 


z = erdi TANE Sı pydi Š St pzdl 
ai Jp pdl ° ™  fepdl’ ™ Irpdi 


(4.1) 


Tem 


(5.1) 


(5.2) 
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Neste caso p é a densidade linear e dl o elemento de comprimento. 

Para corpos homogêneos e uniformes, os fatores densidade (p) são constantes em cada 
caso, e portanto podem ser cancelados em todas as equações acima. 

Se o corpo é composto, isto é se ele é constituido de duas ou mais partes cujos centros de 
massa são conhecidos, então, é claro, da definição do centro de massa que podemos escrever 


AER TIM + Loma +... (5.4) 


Mı +mo +... 


Equações semelhantes fornecem Yem € Zem. Aqui (x1, Y1, 21) são as coordenadas do centro 
de massa do constituinte de massa my, etc. 


Considerações de Simetria 


Se o corpo possui alguma simetria, pode-se tirar vantagem dessa simetria para a localiza- 
ção do centro de massa. Deste modo, se o corpo possui um plano de simetria, isto é, se cada 
partícula m; possui uma imagem especular m/ relativa a um certo plano, então o centro de 
massa encontra-se neste plano. Para provar isto, suponhamos que o plano xy é um plano de 
simetria. Teremos então n 

Zem 5 7 
Emi + mi) 
mas m; = m; e zi = —z;. Então os termos no numerador cancelam-se em pares e Zem = 0; 
ou seja o centro de massa encontra-se no plano xy. 

De modo análogo, se o corpo possui uma linha de simetria, é fácil mostrar que o centro 

de massa é um ponto dessa linha. A prova disto fica como exercício. 


Hemisfério Sólido 


Para encontrar o centro de massa de um hemisfério sólido homogêneo de raio a, nós 
sabemos que, por razões de simetria, o centro de massa encontra-se sobre o raio normal à 
superfície plana do hemisfério. 

Fixamos os eixos coordenados como mostra a Figura 5.1. Temos que o centro de massa 
está sobre o eixo z. Para calcular Zem usaremos elementos de volume circulares de espessura 


dz e raio va? — z? como mostrado. Então 


de modo que 
cedo pTzl@ -— z? )dz 3 


ARS Je orla? — z223)dz 8º = 


Casca Hemisférica 


Para uma casca hemisférica de raio a usamos os mesmos eixos do problema anterior (Fig. 
5.1). Novamente, pela simetria, o centro de massa está localizado sobre o eixo z. Para 
elemento de superfície tomamos uma tira circular de largura ad0. Podemos escrever 


ds=21va2 — z?ad0 
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mas 9 = arcsen (2/a) de modo que 


e portanto ds = 27adz, logo 


= (5.6) 


O 


Figura 5.1: Coordenadas para o cálculo do centro de massa de um hemisfério. 


Aro Semicircular 


Para encontrarmos o centro de massa de um aro semicircular de raio a, usamos os eixos 
mostrados na Figura 5.2. Temos 
dl = ado 


(0) x 
Figura 5.2: Coordenadas para o cálculo do centro de massa de um aro semicircular. 


e z = asen f logo 
— Jo plasend)ado 2a 


Parn JE pad x 


(5.7) 


Lâmina Semicircular 
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No caso de uma lâmina semicircular uniforme, o centro de massa está sobre o eixo z (Fig. 
5.2). Deixamos como problema a prova de que 


4a 


= (5.8) 


Zem 

5.2 Alguns Teoremas sobre o Equilíbrio Estático de um 
Corpo Rígido 

Como já foi visto (Seção 4.1), a aceleração do centro de massa de um sistema é igual à soma 


vetorial das forças externas dividida pela massa. Em particular, se o sistema é um corpo 
rígido e se a soma de todas as forças externas é nula 


Fi +F +... =0 (5.9) 


então o centro de massa, se estiver inicialmente em repouso, permanecerá parado. Desta 
forma a Equação (5.9) expressa a condição de equilíbrio translacional de um corpo rígido. 
Analogamente, o anulamento do torque 


r xFi+r xFo+...=0 (5.10) 


significa que o momentum angular total do corpo não varia (Seção 4.2). Esta é a condição 
de equilíbrio rotacional do corpo rígido que implica no fato de se o corpo estiver inicialmente 
em repouso, não comece a girar. As Equações (5.9) e (5.10) juntas constituem as condições 
necessárias para o equilíbrio completo de um corpo rígido. 


Equilíbrio em um Campo Gravitacional Uniforme 


Consideremos um corpo rígido dentro de um campo gravitacional uniforme, digamos, na 
superfície da Terra. Uma vez que a soma das forças gravitacionais é mg onde m é a massa 
do corpo, podemos escrever a condição para o equilíbrio translacional como 


Fi +F +... +mg= 0 (5.11) 


onde F4, F2, etc, são todas as forças externas além da força gravitacional. Analogamente, a 
condição de equilíbrio rotacional deve ser escrita como 


r x F +i xF +... +A Exmg=o0 (5.12) 
mas g é constante de modo que podemos escrever 


A equação acima estabelece que o momento da força da gravidade sobre todos os pontos 
do corpo é igual ao momento de uma única força mg atuando no centro de massat. A 


10O centro aparente de atuação da força gravitacional é denominado centro de gravidade. Em um campo 
gravitacional uniforme como o que estamos considerando, o centro de massa e o centro de gravidade são 
coincidentes. 
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equação para o equilíbrio rotacional pode ser então escrita como 


r x Fi +r x Fot... +Em X mMmE=0 (5.14) 


Equilíbrio sob Forças Coplanares 


Se as linhas de ação de um conjunto de forças que atuam em um corpo rígido são copla- 
nares, isto é, elas estão todas no mesmo plano, então podemos escrever F =TFa +E, 
etc. 

As equações para o equilíbrio na forma de componentes, Equações (5.9) e (5.10), (que o 
estudante deve lembrar da Física Elementar) são então (equilíbrio translacional) 


Fa + Fr2 +... =0 e Fy + Fyt... =0 (5.15) 
(equilíbrio rotacional) 


Zifya — gif + xao — Yy2Fr2 +... = 0 (5.16) 


5.3 Rotação de um Corpo Rígido em Torno de um Eixo 
Fixo — Momento de Inércia 


O tipo mais simples de movimento de um corpo rígido, além da translação pura, é aquele 
no qual o corpo está obrigado a girar em torno de um eixo fixo. Tomemos o eixo z de um 
sistema de coordenadas apropriado como eixo de rotação. 

A trajetória descrita por uma partícula m; localizada no ponto (x;, Yi, zi) é um círculo 


de raio 4/x? + y? = R; centrado no eixo z. Uma seção transversal paralela ao plano xy é 
mostrada na Figura 5.3. 


Figura 5.3: Seção transversal de um corpo rígido girando em torno do eixo z. 


A velocidade v; da partícula i é dada por 


vi = Rw = y z? + yw (5.17) 
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onde w é a velocidade angular da rotação. Analisando a figura vemos que a velocidade tem 
as seguintes componentes 


ti = —vsend = —wyi (5.18) 
Ji = Vicoso = wri (5.19) 
4 = 0 (5.20) 


onde à é definido como mostra a figura. As equações acima também podem ser obtidas ao 
se tomar as componentes de 


onde w = kw. 
Calculemos a energia cinética de rotação do corpo. Temos 


Res Sma? = ; [> mê) w = Ia? (5.22) 
onde 


A quantidade 7, definida pela equação acima, é de particular importância no estudo de 
corpos rígidos. Ela é denominada momento de inércia. 
Para mostrar como o momento de inércia é importante, calculemos primeiramente o 
momentum angular em relação ao eixo de rotação. 
Uma vez que o momentum angular de uma única partícula é definido como T; x m;v;, a 
componente z é 
mix — yiti) = MmT? + yu = m; R?w (5.24) 


onde fizemos uso das Equações (5.18) e (5.19). A componente z total do momentum angular, 
que vamos representar por L é então dada somando-se sobre todas as partículas 


L=5 miRjw = Iw (5.25) 


Na Seção 4.2 vimos que a taxa de variação do momentum angular de qualquer sistema é 
igual ao momento total das forças externas. 
Para um corpo obrigado a girar em torno de um eixo fixo, temos 
dL d(Iw) 


N Sen 2 
dt dt (mee) 


onde N é o momento total de todas as forças externas relativo ao eixo de rotação (componente 
de N sobre o eixo). Se o corpo é rígido, temos que 1 é constante e podemos escrever 


N=I— (5.27) 


A analogia entre as equações para translação e para rotação em torno de um eixo é 
mostrada a seguir 
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TRANSLAÇÃO ROTAÇÃO 
momentum linear p = mv momentum angular L = Iw 
força F = mò torque N = Iw 
energia cinética T = mw? energia cinética T = ilw? 


Desta forma, o momento de inércia é análogo à massa: ele é uma medida da inércia 
rotacional do corpo em relação a um dado eixo fixo de rotação, do mesmo modo como a 
massa é a medida da inércia translacional do corpo. 


5.4 Cálculo do Momento de Inércia 


Nos cálculos práticos do momento de inércia 5) mRĉ?, para corpos extensos, substituimos 
a somatória por uma integração por todo o corpo, do mesmo modo como fizemos para o 
cálculo do centro de massa. Desta forma podemos escrever 


N= / R?dm (5.28) 


onde dm, elemento de massa, é dado pelo fator densidade multiplicado por uma diferencial 
apropriada (volume, área ou comprimento). É importante lembrar que R é a distância 
perpendicular entre o elemento de massa e o eixo de rotação. 
No caso de um corpo composto, é claro, da definição de momento de inércia, que podemos 
escrever 
fo Td dados (5.29) 


onde J, 15, etc, são os momentos de inércia dos vários i constituintes, relativos ao eixo 
escolhido. 
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Calculemos o momento de inércia para alguns casos particulares importantes. 


Barra Fina 


Para uma barra fina uniforme de comprimento a e massa m, temos, para o eixo perpen- 
dicular à barra passando por uma extremidade (Figura 5.4 (a)) 


Àz AZ 


EEE 


Figura 5.4: Coordenadas para o cálculo do momento de inércia de uma barra (a) em relação 
a uma extremidade (b) em relação ao centro. 


3 


a 1 
I= / vpda = pé = ma? (5.30) 
0 3 3 


a última passagem resulta do fato de que pa = m. 
Se o eixo é perpendicular à barra mas passa pelo seu centro (Fig. 5.4 (b)) temos 
+3 1 1 
I= : Tadaa Emo? 5.31 
‘aT pdz poep (5.31) 


paa 
2 


Aro ou Casca Cilíndrica 


No caso de um aro circular fino ou de uma casca cilíndrica, para o eixo central ou de 
simetria, todas as partículas se encontram à mesma distância do eixo de modo que 


I = ma? (5.32) 


onde a é o raio e m a massa. 


Disco Circular ou Cilindro 


Para calcular o momento de inércia de um disco circular uniforme de raio a e massa m, 
vamos usar coordenadas polares. O elemento de massa, um aro fino de raio r e espessura dr 
é dado por 

dm = p2rrdr 
onde p é a densidade superficial. O momento de inércia em relação ao eixo que passa pelo 
centro do disco, normal ao plano das faces (Figura 5.5) é obtido do seguinte modo 


a 4 
I= / p(r?)(2xrdr) = 2rp— = “ma? (5.33) 
0 4 
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A última passagem envolveu a relação ma2p 


m. 


7 


eixo 


Figura 5.5: Coordenadas para o cálculo do momento de inércia de um disco. 


Obviamente a Equação (5.33) também se aplica para um cilindro de raio a e massa m 
para o eixo central do mesmo. 


Esfera 


Encontremos o momento de inércia de uma esfera sólida uniforme de raio a e massa 
m relativo a um eixo (eixo z) passando por seu centro. Iremos dividir a esfera em discos 
circulares finos, como mostra a Figura 5.6. 

O momento de inércia de um desses discos de raio y é dado pela Equação (5.33), ty’dm 
mas dm = pry?dz de modo que 


Figura 5.6: Coordenadas para o cálculo do momento de inércia de uma esfera em relação ao 
eixo z. 


a] al 8 
LS J- Toy dz = je zT — 2 dz = g (5.34) 
A última passagem deve ser completada pelo leitor. Uma vez que m = ira?p temos 
2 


L= zma? (5.35) 
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Casca Esférica 
O momento de inércia de uma casca esférica fina e uniforme pode ser obtido facilmente 


pela aplicação da Equação (5.34). Se diferenciamos em relação a a, temos 


dl = ` npada 


este resultado é o momento de inércia de uma casca de raio a e espessura da. A massa da 
casca é 4ra? pda de modo que 


2 
I= ama? (5.36) 
para o momento de inércia de uma casca de massa m e raio a. O leitor pode também tentar 
obter este resultado por integração. 
Teorema dos Eixos Perpendiculares 


Considere um corpo rígido na forma de uma chapa plana e fina de formato qualquer. 
Coloquemos esta lâmina no plano xy (Figura 5.7). O momento de inércia relativo ao eixo z 


é dado por 
L= Y mile; +y) = Y mi? + Y my; 


Z 


Figura 5.7: Teorema dos eixos perpendiculares. 


Mas X; mz? é exatamente o momento de inércia 1, desse corpo com relação ao eixo y 
pois z; vale zero para todas as partículas. 
Do mesmo modo, 35; m;y? é o momento de inércia I, desse corpo em relação ao eixo z e 
a equação acima pode ser escrita 
l=; Iy (5.37) 


Este é o teorema dos eixos perpendiculares. Em palavras: o momento de inércia de 
qualquer lâmina plana, relativo a qualquer eixo normal ao seu plano, é igual à soma dos 
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momentos de inércia relativos a dois eixos quaisquer mutuamente perpendiculares, contidos 
no plano da lâmina e que interceptam o primeiro eixo. 

Como exemplo do uso deste teorema, vamos considerar um disco circular fino no plano 
xy (Figura 5.8). Da Equação (5.33) temos 


1 
I= ama” =l; + Iy 


ENN 
Ny 


0, 
\ EM 
N pa 


Figura 5.8: Disco circular fino no plano xy. 


| 


A 
= 


Ate 


N 


e. 


Neste caso, no entanto, sabemos por simetria que 1, = I}. De modo que 


1 
b=; qro (5.38) 


para o momento de inércia de um disco, relativo a um eixo contido no plano do disco e que 
passa pelo seu centro. A Equação (5.38) pode ser também obtida por integração direta. 


Teorema dos Eixos Paralelos 


Tomemos a equação do momento de inércia relativo a um eixo qualquer, por exemplo o 
eixo z 


1=Dmi(o? +18) 


Agora podemos expressar x; e y em termos das coordenadas do centro de massa 
(Zem, Yem, Zem) € das coordenadas relativas do centro de massa (%;, J; Zi) (Figura 5.9) 
da seguinte maneira 

Ti = Lem + Ti e Yi = Yem + J; (5.39) 


temos então, depois de substituições e rearranjos de termos 
LaS mla AT) + e + Yom) + n S ate + 2Yem D MiT; (5.40) 


O primeiro termo à direita é o momento de inércia relativo a um eixo paralelo ao eixo z que 
passa pelo centro de massa. Vamos representá-lo por L.m. A segunda soma à direita é igual 
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Z 


Figura 5.9: Teorema dos eixos paralelos. 


ao produto da massa total do corpo pelo quadrado da distância do centro de massa ao eixo 
z. Representamos esta distância por l. Isto é 2 = £2m + Yn- 
Agora, da definição do centro de massa 


De modo que as duas últimas somas da direita da Equação (5.40) são nulas e o resultado 
final pode ser escrito 
I= Im + M? (5.41) 


Este é o teorema dos eixos paralelos. Ele pode ser aplicado a qualquer corpo rígido 
tanto sólido quanto laminar. Este teorema estabelece que o momento de inércia, de um 
corpo rígido, relativo a um eixo qualquer é igual ao momento de inércia relativo a um eixo 
paralelo ao primeiro, passando pelo centro de massa, mais o produto da massa do corpo pelo 
quadrado da distância entre os eixos. 

Aplicando o teorema acima a um disco circular, temos, das Equações (5.33) e (5.41) 


1 3 
I= ama? + ma” = ama? (5.42) 


para o momento de inércia de um disco circular uniforme, relativo a um eixo perpendicular 
ao seu plano, passando por sua borda. 
Por outro lado, das Equações (5.38) e (5.41) temos 


1 
I= qa + ma? 


para o momento de inércia do mesmo disco, relativo a um eixo paralelo ao seu plano, tangente 
à borda. 


Raio de Giração 
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Com diversos propósitos, é conveniente se expressar o momento de inércia de um corpo 
rígido em termos da distância k denominada raio de giração, onde k é definido pela equação 


I= mk’ ou k= E (5.43) 
m 


Por exemplo, o raio de giração de uma barra fina, relativo a um eixo perpendicular que passa 
por uma extremidade (veja a Equação 7.30) é 
Mia dr 
m v3 
Os momentos de inércia de vários objetos podem ser tabelados simplesmente através dos 
quadrados de seus raios de giração. 


k= 


5.5 O Pêndulo Físico 


Um corpo rígido que está livre para oscilar, sob ação de seu próprio peso, em relação a um 
eixo horizontal fixo é denominado pêndulo físico ou pêndulo composto. Um pêndulo 
físico é mostrado na Figura 5.10 onde O representa a localização do eixo de rotação, e cm é 
o centro de massa. A distância entre O e cm é | como pode ser visto. 

Representando o ângulo entre a linha Ocm e a vertical OA por 0, o momento da força 
gravitacional (atuante no cm) relativo ao eixo de rotação tem módulo 


mgl sen O 
A equação fundamental de movimento 
N = Íw 
pode então tomar a forma , 
—mgl sen 0 = I0 
ou l 
Ö + — sen O = 0 (5.44) 


A equação acima é formalmente idêntica à equação do movimento de um pêndulo simples. 
Para pequenas oscilações, como no caso do pêndulo simples, podemos substituir sen 0 por 0 
ô + To =0 (5.45) 

A solução é 
0 = do cos(27 ft + e) (5.46) 


onde ĝo é a amplitude e € é o ângulo de fase. A frequência da oscilação f é dada por 


Ma (5.47) 


T= zN T 
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Figura 5.10: O pêndulo físico. 


O período T é portanto 
I 
T = 2m4 | — (5.48) 
mgl 
(Para evitar confusão, não usaremos um símbolo específico para representar a frequência 
angular cujo valor é dado pelo produto 2r f). A expressão do período em termos do raio de 


giração k é 
k2 


Desta forma, vemos que o período de um pêndulo composto é igual ao de um pêndulo 
simples de comprimento k?/l. Como exemplo, podemos analisar o caso de uma barra fina e 
uniforme de comprimento a oscilando em torno de um de seus extremos (k? = a2/3). Seu 
período é 


Centro de Oscilação 


Usando o teorema dos eixos paralelos, podemos expressar o raio de giração k em termos 

do raio de giração relativo ao centro de massa kem 

I= Imt mM? 
ou 

mk? = mk2, + ml? 

Cancelando as massas 

k? = kla tl (5.50) 
A Equação (5.49) pode então ser escrita 

k?n +12 


T = 274| 5.51 
[es (5.51) 
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Suponha que o eixo de rotação de um pêndulo físico seja transferido para a posição O' 
à distância |” do centro de massa, como mostra a Figura 5.10. O período da oscilação T” 
relativo ao novo eixo é dado por 


ET 


T' = 2r 
gl! 


(5.52) 


Segue-se que os períodos de oscilação em torno de O e em torno de O” serão iguais se 


ken tP kên +? 
gl ql 


que se reduz a 
W = k?n (5.53) 


O ponto O' é denominado centro de oscilação do ponto O. É claro que O também é centro 
de oscilação do ponto O". Desta maneira, em torno de uma extremidade de uma barra temos 
k?„ = a?/12 e l =a/2. Usando a Equação (5.53) encontramos l’ = a/6 e portanto a barra 
terá o mesmo período quando o eixo for transferido da extremidade para a posição O’ situada 
à distância a/6 do seu centro. 


5.6 Teorema Geral Relativo ao Momentum Angular 


Com o objetivo de se estudar o caso mais geral do movimento de um corpo rígido, no qual 
o eixo de rotação não é fixo, precisamos desenvolver um teorema fundamental relativo ao 
momentum angular. Como vimos na Seção 4.2, a taxa de variação temporal do momentum 
angular de qualquer sistema é igual ao torque aplicado 

dL 


RSN 5.54 
E (5.54) 


ou explicitamente 

=> (E; x mvi) = X(T; x F;) (5.55) 

A 2 

Nas equações acima, todas as quantidades se referem a um sistema inercial de coordenadas. 
Agora vamos introduzir o centro de massa para expressar o vetor posição de cada partícula 

r; em termos da posição do centro de massa Yem e do vetor posição da partícula i em relação 

ao centro de massa T; (como na Seção 4.3) 


A Equação (5.55) fica 


1 1 
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Desenvolvendo esta expressão e usando o fato de que tanto  m;i; quanto X m;v; são nulas, 
podemos reduzir a Equação (5.56) a 


onde aim = Vem. 


Na Seção 4.1 mostramos que a translação do centro de massa de qualquer sistema de 
partículas obedece à equação 


Consequentemente o primeiro termo à esquerda da Equação (5.57) cancela-se com o primeiro 
termo à direita. O resultado final é 


A soma à esquerda na equação acima é exatamente o momentum angular do sistema 
e a soma à direita é o momento total das forças externas relativo ao centro de massa. 


Representando estas quantidades por Le N, respectivamente, temos 


| 
Z| 


(5.60) 


Se Eu 


Este resultado estabelece que a taxa de variação temporal do momentum angular relativo 
ao centro de massa de qualquer sistema é igual ao momento total das forças externas relativo 
ao centro de massa. Isto é verdadeiro mesmo que o centro de massa esteja se acelerando. Se 
tomarmos qualquer outro ponto além do centro de massa como ponto de referência, então 
tal ponto precisa estar em repouso em um sistema de referência inercial (exceto para certos 
casos especiais para os quais não vamos nos reter em discussão). 

Um exemplo do uso deste teorema é dado na Seção 5.8. 


5.7 Movimento Laminar de um Corpo Rígido 


Se o movimento do corpo ocorre de modo que todas as suas partículas se deslocam para- 
lelamente a um determinado plano fixo, então este movimento é denominado laminar. No 
movimento laminar, o eixo de rotação pode mudar de posição mas não muda de direção. A 
rotação em torno de um eixo fixo é um caso especial de movimento laminar. Outro exemplo 
é o de cilindro que rola sobre um determinado plano. 

Se um corpo possui um deslocamento laminar, tal deslocamento pode ser especificado 
como se segue. 

Escolhemos um ponto de referência no corpo, por exemplo, o centro de massa; o ponto 
em questão se desloca de Ar, ao mesmo tempo o corpo gira em relação ao ponto de referência 
de um ângulo Ag. É claro que qualquer deslocamento laminar pode ser especificado desta 
maneira. Consequentemente, o movimento laminar pode ser especificado conhecendo-se a 
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velocidade translacional de um ponto de referência adequado, juntamente com a velocidade 
angular. 
A equação fundamental que governa a translação de um corpo rígido é 


> 


F = mim = Mem = Mem (5.61) 


onde F representa a soma de todas as forças externas que agem sobre o corpo, m é a massa 
e dem à aceleração do centro de massa. 
A aplicação da Equação (5.25) ao caso de um movimento laminar de um corpo rígido nos 
fornece 
L= Lato (5.62) 


para o módulo do momentum angular relativo a um eixo Č que passa através do centro de 
massa, onde w é a velocidade angular da rotação em torno desse eixo. 
A equação fundamental que governa a rotação do corpo, Equação (5.60), fica então 


> 


dL + g 
SE DN 5.63 
H ui (5.63) 


onde N é o momento total das forças aplicadas, relativo ao eixo C. 


5.8 Corpo Rolando em um Plano Inclinado 


Para ilustrar o movimento laminar, vamos estudar o movimento de um objeto que possua 
seção circular (cilindro, bola, etc ....) rolando para baixo em um plano inclinado. Como 
mostra a Figura 5.11, existem três forças agindo no corpo. Elas são: (1) a força da gravidade; 
(2) a força de reação normal do plano sobre o corpo (Fy) e (3) a força de atrito (F), paralela 
ao plano. 


Figura 5.11: Corpo rolando em um plano inclinado. 


Escolhendo os eixos como mostra a Figura, as equações das componentes do movimento 
translacional do centro de massa são 
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Mem = mgsenO — F (5.64) 
Mijem = —mgcos6 + Fy (5.65) 


onde 6 é o ângulo de inclinação do plano em relação à horizontal. Uma vez que o corpo 
permanece em contato com o plano, 


Yem = constante 
e portanto tem = 0, de modo que, da Equação (5.65) temos 
Fy = mgcos6 (5.66) 


A única força que provoca torque em relação ao centro de massa é a força de atrito. O 
módulo desse momento é Fa onde a é o raio do corpo. Deste modo, a equação rotacional 
(5.63) fica 

Int = Fa (5.67) 


Para discutir o movimento mais adiante, precisamos fazer algumas suposições relativas ao 
contato do corpo com o plano. Vamos resolver as equações do movimento para dois casos. 
Movimento sem deslizamento 


Se o contato é muito rugoso tal que não ocorre deslizamento do corpo, temos as seguintes 
relações 


tem = GO 
iem = aġ = aw (5.68) 
Lemn = aĝ = QW 


onde à é o ângulo associado à rotação. A Equação (5.67) pode então ser escrita 


Tem 


Substituindo o valor de F (acima) na Equação (5.64) temos 
.. Tm .. 
Mem = Mg sen O — — Tem 
a 


resolvendo para %.m encontramos 


mgsenô  gsenb 
m + (Ima) 1+ (kên/02) 


Tem EE 


(5.70) 


onde kem é O raio de giração em relação ao centro de massa. O corpo rola, portanto, com 
aceleração linear constante e, por consequência da Equação (5.68), com aceleração angular 
constante. 
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Por exemplo, a aceleração linear de um cilindro uniforme (k2, = a2/2) é 


Considerações sobre energia 


Os resultados acima podem ser também obtidos à partir de considerações energéticas. 
Em um campo gravitacional uniforme, a energia potencial V de um corpo rígido é dada pela 
soma das energias potenciais individuais de todas as partículas 


V = X (migz) = mg2Zem 


onde Zem é a distância vertical do centro de massa a um plano arbitrário de referência. 
Se as outras forças, além da gravitacional, não realizarem trabalho, então o movimento é 
conservativo e podemos escrever 


T +V =T + mMgzem = E = constante 


onde T é a energia cinética. 
No caso de um corpo rolando num plano inclinado, Figura 5.11, a energia cinética de 
translação é mi? /2.e a de rotação é Iemw?/2 de modo que a equação da energia fica 


-2 
3"Lem 3 omw Emil = É 
mas w = Zem/G € Zem = —LTem Sen 8; logo 
1. 1 r? 
ope — grem sen = E 


gem l 9 “em a2 


No caso de rolamento puro a força de atrito não afeta a equação de energia uma vez que é 
perpendicular ao deslocamento e consequentemente não realiza trabalho; logo E é constante. 
Diferenciando em relação ao tempo e agrupando os termos temos 


k2 
Tanta (i + t) — gJžem SEn 0 = 0 
a 


Cancelando o fator comum tem (supondo é claro que tem # 0) e resolvendo para žem encon- 
tramos o mesmo resultado obtido anteriormente com o uso de forças e momentos. 


Ocorrência de Deslizamento 
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Consideremos agora o caso em que o contato com o plano não é perfeitamente rugoso 
mas tem um certo coeficiente de atrito de deslizamento u. Havendo este deslizamento, o 
módulo da força de atrito F é 

F = uFy = umg cos0 (5.71) 


A equação da translação (5.64) fica então: 
Mem = Mg send — umg cos O (5.72) 
e a equação de rotação (5.67) fica 
Tem = umga cos 6 (5.73) 
Da Equação (5.72) vemos novamente que o centro de massa possui aceleração constante 
Žem = g(sen 0 — u cos 0) (5.74) 


e, ao mesmo tempo, a aceleração angular é constante 


o — UMA cos 0 _ uga cos 0 (5.75) 


I cem kon 


Integrando estas duas equações em relação ao tempo e, supondo que o corpo parta do 
repouso (isto é: em t = 0, tem = 0, dem = 0) temos 


tem = g(sen6 — ucos0)t 
(5.76) 


Consequentemente as velocidades linear e angular têm razão constante e podemos escrever 
Tem = yaw 
onde 


m ua? cos 0/k2n a? 


os 2 
senB— cos k?n, (= £ 1) (5.77) 


E 7 
Mas, aw não pode ser maior que tem, de modo que y não pode ser menor que a unidade. 
O caso limite, para o qual não existe deslizamento, é dado por tem = aw ou seja y = 1. 
Resolvendo a Equação (5.77) para u, com y = 1, obtemos o valor crítico de u que é dado 


por 
tan 6 


e a 
Se u for maior que este valor, o corpo rola sem deslizar. Por exemplo, se uma bola é 


colocada em um plano com 0 = 45º, ela irá rolar, sem deslizar, se o coeficiente u for maior 
que tan 45º /(1 + 5/2) ou 2/7. 


(5.78) 
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5.9 Movimento de um Corpo Rígido Sob a Ação de 
uma Força Impulsiva 


No Capítulo anterior introduzimos o conceito de uma força impulsiva agindo em uma partícula. 
Vimos que o efeito de tal força, ou impulso, é o de produzir uma mudança súbita na veloci- 
dade da partícula. Nesta Seção vamos estender o conceito de impulso ao caso do movimento 
laminar de um corpo rígido. 


Movimento Livre 


Suponhamos que um corpo se encontre livre para movimentar-se em um plano e esteja 


submetido a um impulso P. Nestas circunstâncias, de acordo com a teoria geral discutida 
na Seção 5.7, temos que considerar tanto a translação quanto a rotação do corpo. 
Primeiramente, a translação é descrita pela fórmula geral 


F=mvYm 
se F é uma força impulsiva, temos 


f P= smi 


De modo que o resultado do impulso é causar uma variação na velocidade do centro de massa 
igual a 


3 | e 


AV em 


(5.79) 
Em segundo lugar, a rotação do corpo é governada pela equação 

N=L = Ini 
Integrando em relação ao tempo obtemos a relação 


j: Ndt = Im Að (5.80) 


Chamamos a integral acima de impulso rotacional. Usemos o símbolo L para designá-lo. 
O efeito de um impulso rotacional é pois o de causar uma variação na velocidade angular do 
corpo de uma quantidade 


Ep 


Ağ = (5.81) 


Teri 


Se o impulso primário P é aplicado ao corpo de um modo que sua linha de ação está à 
distância b do centro de massa, o momento da força (N = b x F), e consequentemente 


A 
> 


L=bx 


Fo!» 


(5.82) 
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Podemos então expressar a variação na velocidade angular produzida por um impulso 

como A 

b xP 

Jpm 

Em resumo: o efeito de um impulso sobre um corpo rígido, livre para deslocar-se em movi- 

mento laminar é (1) produzir uma mudança súbita na velocidade do centro de massa — efeito 

translacional e (2) produzir uma mudança súbita na velocidade angular do corpo — efeito 
rotacional. 


Ag = (5.83) 


Movimento Forçado 


No caso em que o corpo submetido a uma força impulsiva não está livre, mas é obrigado 
a girar em torno de um eixo fixo, devemos considerar apenas a condição de rotação N = Iw 
de modo que 


fNdt=i =1A 


Na equação acima I é o momento de inércia em torno do eixo fixo de rotação e N é o 
momento em relação a esse eixo. Neste caso, o impulso rotacional L, que é produzido pelo 
impulso primário simples P cuja linha de ação está à distância b do eixo de rotação, é dado 
por 


L = Pb 
de modo que i 
Aw = o (5.84) 


é a variação da velocidade angular em torno do eixo fixo de rotação. 


Efeito de Vários Impulsos Simultâneos 


Se um certo número de impulsos diferentes é aplicado simultaneamente em um corpo 
rígido, as variações resultantes da velocidade do centro de massa e da velocidade angular do 
corpo são obtidas pelas adições dos impulsos e momentos, respectivamente. Desta forma, o 
efeito translacional de vários impulsos simultâneos é obtido pela soma vetorial dos impulsos 
individuais de modo que a Equação (5.79) fica 


A 


EA dera 
Area O (5.85) 


m 
De modo análogo, para o efeito rotacional, a Equação (5.83) se modifica para 


Pbi + Pabat... 


A 
£ a 


(5.86) 


No caso de um corpo que é obrigado a girar em torno de um eixo fixo, existe um impulso 
secundário devido à reação do eixo sobre o corpo quando um impulso externo é aplicado. O 
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movimento é então determinado pela soma de todos os impulsos de acordo com as equações 
acima. 


Exemplos 


1. Impulso Aplicado em uma Barra Livre 


Como ilustração da teoria descrita acima, consideremos uma barra que está livre para 


deslizar sobre uma superfície lisa horizontal. Seja P o impulso aplicado à barra à distância 
b do seu centro de massa em uma direção perpendicular ao comprimento da barra, como 
mostra a Figura 5.12. 


"> 


2a 


Figura 5.12: Impulso aplicado a uma barra livre. 


Se a barra estiver inicialmente em repouso, as equações para a translação e para a rotação 
são respectivamente 


P 
Vem e e 5.87 
em = — (5.87) 
Pb 
= 5.88 
w E (5.88) 
Em particular, se a barra é uniforme de comprimento 2a, então Tm = ma?/3 
a 3b 
= pe. 5.89 
e (5.89) 


e portanto, a velocidade adquirida pelo centro de massa é a mesma qualquer que seja o ponto 
de aplicação do impulso enquanto que a velocidade angular adquirida pela barra depende do 
ponto de aplicação do impulso. Vimos também que a energia cinética final da barra é 


T=5 á Imu? = — Sai 
tn O om 


1 1 Ê? 3/6) 
a 


e é claro que ela depende do ponto de aplicação do impulso. 


2. Impulso Aplicado a uma Barra Obrigada a Girar em Torno de um Eixo Fixo 
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Vamos considerar o caso em que a mesma barra do exemplo 1 esteja obrigada a girar em 
torno de um eixo fixo. Suponhamos que o eixo O esteja localizado em um dos seus extremos 
como mostra a Figura 5.13. Neste caso, temos a seguinte equação rotacional 


L=P(a+b)=Iw (5.90) 
Como T, = (4/3)ma? obtemos 
_ g3(a+b) 


para a velocidade angular adquirida pela barra. Agora, uma vez que a barra gira em torno 
de O, o centro de massa está em movimento. Sua velocidade é 


Vem = AW 
ou 3( b) 
K a + 
É 4ma l ) 


Notamos que é diferente de P/m. A primeira vista este resultado parece contradizer a 
equação geral para translação, Equação (5.79). Na verdade não há contradição porque 
existe um outro impulso que age sobre a barra simultaneamente, como já foi mencionado. 
Este é o impulso de reação exercido sobre a barra pelo eixo no ponto O. Vamos denotá-lo 


por P,. A velocidade adquirida pelo centro de massa é 


P +P, 
E a (5.93) 
m 
A U 
bp ~> T 


a+b 


1 


Figura 5.13: Impulso aplicado a uma barra que é obrigada a girar em torno de uma de suas 


extremidades. O impulso de reação no eixo é Po. 
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Podemos então calcular P, usando o valor de vem dado pela Equação (5.92). Desta forma 
temos 


que nos leva a 


(5.94) 


Para o impulso recebido pela barra devido à reação do eixo. Pela lei da ação e reação, o 
impulso recebido pelo eixo é =P. 

Deve-se notar que o impulso de reação pode ser nulo caso se aplique o impulso primário 
em um ponto escolhido adequadamente. Este ponto é denominado centro de percussão. 
No presente caso, este ponto é tal que b = a/3. 


5.10 Colisões de Corpos Rígidos 


Nos problemas que envolvem colisões de corpos rígidos, as forças que os corpos exercem 
uns sobre os outros durante o contato são sempre iguais e opostas. Consequentemente, 
os princípios de conservação de momentum linear e momentum angular são válidos. Os 
conceitos de impulsos lineares e rotacionais são de grande valia nestes problemas. 


Exemplo 
Colisão entre uma Bola e uma Barra 


Considere, por exemplo, o impacto de uma bola de massa m” com uma barra uniforme 
de comprimento 2a e massa m. Vamos supor que a barra esteja inicialmente em repouso em 
uma superfície lisa horizontal, como nos casos anteriores, e que o ponto de impacto seja à 
distância b do centro de massa como mostra a Figura 5.14. 


> Je 2a 


Figura 5.14: Colisão de uma partícula e uma barra. 


As Equações (5.87) e (5.88) dão o movimento da barra após o impacto em termos do 


impulso P aplicado pela bola sobre a barra. Sabemos também que o impulso recebido pela 
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bola neste impacto é —P. Podemos então escrever as equações para a translação 


Mem (5.95) 


o Tp 


onde Vem é à velocidade do centro de massa da barra após o impacto, vo é a velocidade 
da bola antes do impacto e Vı é a velocidade final da bola. O conjunto das duas equações 
relativas à translação implicam na conservação do momentum linear pois, eliminando P das 
duas temos 

m'vo = mY + MYem (5.97) 


Com o objetivo de determinarmos a rotação da barra após o impacto, podemos usar o 
princípio de conservação do momentum angular. O momentum angular inicial da bola em 
relação ao centro de massa da barra é bm'vo e o final é bm'v,. Para a barra o inicial é zero 
e o final é Iemw. Deste modo 


bm'vo = bmv + Iemw (5.98) 


As equações acima (translacional e rotacional) não são suficientes para calcularmos as três 
velocidades (v1, Vem € w) do movimento final. 

Para calcularmos este movimento completamente, necessitamos ainda de outra equação. 
Esta pode ser a equação do balanço de energia 


mv, + moi, + Lato” + (5.99) 


em que Q é a perda de energia devida à colisão. Alternativamente, podemos usar a equação 
do coeficiente de restituição 


velocidade de afastamento 


ET A x = 
velocidade de aproximação 


No problema em questão temos 
velocidade de aproximação = vo 


Para calcularmos a velocidade de afastamento, precisamos conhecer a velocidade da barra 
no ponto de impacto. Esta é dada pela soma das velocidades, translacional do centro de 
massa e a rotacional daquele ponto relativa ao centro de massa. 

Deste modo, a velocidade do ponto de impacto imediatamente após a colisão é vem + bw. 
Podemos então escrever 


velocidade de afastamento = Vem + bw — vo 


Portanto 
Evo = Vem + bw — vo 
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Temos agora equações em número suficiente para resolver o problema. Fazendo L.m = 
ma? /3 obtemos 


VI = vo = a tem (5.100) 


3b 
W = Vem E 


O leitor deve verificar estes resultados. 


Exercícios de Treinamento 


9.1 


5.2 


5.3 


5.4 


5.9 


Ache o centro de massa de: 
a) Um fio fino e uniforme dobrado na forma de um U em que cada seção de reta tem 
o mesmo comprimento a. 
b) A área limitada pela parábola y = ax? e a reta y = b. 
c) O volume limitado pelo parabolóide de revolução z = a(x? + y?) e o plano z = b. 
d) Um cone reto circular uniforme de altura h. 


Uma esfera sólida uniforme de raio a contém uma cavidade esférica de raio b centrada 
a uma distância c do centro da esfera, em que a > (b+c). Ache o seu centro de massa. 


Ache os momentos de inércia de cada uma das figuras do exercício 5.1, em relação a 
seus eixos de simetria. 


Ache o momento de inércia da esfera do exercício 5.2 em relação ao eixo que passa 
pelos centros da esfera e da cavidade. 


Mostre que o momento de inércia de um octante sólido uniforme de uma esfera de raio 
a é 2ma?/5, em relação a um eixo que passa por uma de suas arestas (Nota: esta é a 
mesma expressão para o momento de inércia relativo ao centro de uma esfera sólida 
de mesmo raio). 


Problemas 


5.6 


5.7 


Um arame de forma semi-circular pende apoiado em um prego. A linha que une seus 
extremos faz um ângulo 0 com a horizontal e o arame está a ponto de escorregar. Qual 
é o coeficiente de atrito entre o arame e o prego? 


Um hemisfério sólido uniforme apoia-se em uma parede vertical e está prestes a es- 
corregar. O lado redondo do hemisfério está em contato com a parede e o chão. Se o 
coeficiente de atrito u é o mesmo para a parede e o chão, encontre o ângulo entre a 
face plana do hemisfério e o chão. 
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5.8 Uma casca hemisférica uniforme apoia-se sobre um plano inclinado de ângulo 0 e está 
prestes a perder o equilíbrio. O lado redondo da casca está em contato com o plano e 
o coeficiente de atrito é u. Encontre a inclinação da casca. 


5.9 Sabe-se que um conjunto de forças E, Eo,... age em certo corpo rígido de modo que 
ele está (a) em equilíbrio translacional e (b) em equilíbrio rotacional em relação a certo 
ponto O. Prove que esse conjunto de forças está também em equilíbrio rotacional em 
relação a qualquer outro ponto O”. 


5.10 Mostre que os momentos de inércia de um sólido paralelepipédico regular uniforme, de 
um cilindro elíptico, e de um elipsóide são respectivamente, (m/3) (a? +62): (m/4) (a? + 
b?) e (m/5) (a? +b?) onde m é a massa e 2a e 2b são os diâmetros principais do sólido, 
perpendiculares ao eixo de rotação sendo que este último passa pelo centro em cada 
caso. 


5.11 Um aro circular de raio a oscila como pêndulo físico em torno de um ponto na sua 
circunferência. Ache o período de oscilação se o eixo de rotação é (a) normal ao plano 
do aro; (b) paralelo ao plano do aro. 


5.12 Mostre que o período de um pêndulo físico é igual a 27(d/9)!2 em que d é a distância 
entre o ponto de suspensão O e o centro de oscilação O. 


5.13 Uma bola sólida uniforme está enrolada por algumas voltas de corda leve. Se a ex- 
tremidade da corda é mantida fixa e a bola cai sob a ação da gravidade, qual é a 
aceleração do centro da bola? 


5.14 Dois homens seguram as extremidades de uma prancha uniforme de comprimento l e 
massa m. Mostre que, se um dos homens solta a sua extremidade repentinamente, o 
peso suportado pelo outro cai de repente de mg/2 para mg/4. Mostre também que a 
aceleração vertical da extremidade livre é 39/2. 


5.15 Uma bola sólida uniforme contém uma cavidade oca esférica no seu centro, cujo raio é 
a metade do da bola. Mostre que a aceleração da bola ao rolar por um plano inclinado 
áspero é 98/101 da aceleração da bola sem cavidade (Nota: isto sugere um método 
para teste não destrutivo). 


5.16 Dois pesos de massas my e ma estão amarrados aos extremos de um fio leve e inex- 
tensível. O fio passa em uma roldana de raio a e momento de inércia J. Encontre a 
aceleração dos pesos (considere mı > ms) desprezando o atrito no eixo da roldana. 


5.17 Um cilindro reto circular uniforme de raio a é equilibrado em cima de um cilindro 
áspero e fixo de raio b (b > a), os eixos dos dois cilindros são paralelos. Se o equilíbrio 
é levemente perturbado, encontre o ponto em que o cilindro móvel perde o contato 
com o fixo. 


5.18 Uma escada está encostada em uma parede vertical lisa. Se o chão também é liso e o 
ângulo inicial entre a escada e ele é 05, mostre que a escada ao cair deslizando perde 
contato com a parede quando o ângulo entre o chão e a escada é arc sen(2/3 sen 6o). 
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5.19 Uma barra comprida e uniforme de comprimento l está parada na vertical, sobre um 
assoalho áspero. A barra é ligeiramente perturbada e cai ao chão. (a) Encontre as 
componentes vertical e horizontal da reação do piso como função do ângulo 6 entre a 
barra e a vertical em um instante qualquer. (b) Encontre também o ângulo para o 
qual a barra começa a escorregar. O coeficiente de atrito entre a barra e o chão é u. 


5.20 Uma bola de bilhar de raio a está inicialmente girando em torno de um eixo horizontal 
com velocidade angular wọ e com velocidade inicial de translação nula. Se o coeficiente 
de atrito entre a bola e a mesa é u, ache a distância percorrida pela bola até desaparecer 
o seu deslizamento. 


5.21 Uma bola é jogada com velocidade vo para subir um plano inclinado áspero com 
inclinação 6 e coeficiente de atrito u. Encontre a posição da bola como função do 
tempo e determine a posição da bola quando passar a ocorrer rolamento puro (sem 
deslizamento). Considere que u é maior que (2/7) tan 0. 


5.22 (a) Um aro circular uniforme repousa sobre uma superfície horizontal livre. Se ele 
é golpeado tangencialmente em um ponto da periferia, em torno de qual ponto ele 
começa a girar? (b) Ache a altura em que uma bola de bilhar deve ser golpeada para 
que role sem deslizar. 


5.23 Mostre que o centro de oscilação de um pêndulo físico é tambem o centro de percussão 
para um impulso aplicado no eixo de rotação. 


9.24 Um pêndulo balístico é feito de uma prancha comprida de comprimento | e massa m. 
Ele se encontra livre para oscilar em torno de um extremo O, e está inicialmente em 
repouso na vertical. Uma bala de massa m’ é disparada horizontalmente atingindo o 
pêndulo à distância | do ponto O vindo a se alojar no pêndulo. Se a amplitude da 
oscilação é ĝo, encontre a velocidade da bala. 


5.25 Duas barras uniformes AB e BC de massas iguais (m) e mesmos comprimentos (l) 
estão levemente ligadas no ponto B. O sistema está inicialmente em repouso sobre 
uma superfície lisa horizontal, os pontos A, Be C estão alinhados. Se um impulso p é 
aplicado em 4, em ângulo reto com as barras, encontre o movimento inicial do sistema. 
(Sugestão: isole as barras). 


5.26 Resolva o problema acima para o caso em que as barras estejam em ângulo reto entre 
si. 


Respostas de Problemas Impares Selecionados 
5.1 (a) a/3 do centro inferior 
b) 3b/5 do vétice 

c) 2b/3 do vértice 

d 


( 
( 
( 
(d) 3h/4 do vértice 
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5.3 (a) Tma?/36 onde m = massa total 
(d) 3ma?/10 onde a é o raio da base. 


5.7 senf = 8u(1 + n)/3(1+ q?) 


5.11 27,/2a/9, 2171/30/29 


5.13 59/7 
5.17 Quando a linha dos centros faz um ângulo cos” !(4/7) com a vertical 


5.19 (a) Horizontal: (3/4)mgsen 0(3 cos 0 — 2) 
Vertical: (1/4)mg(3 cos 6 — 1)2 
(b) O escorregamento começa quando 
|3 sen 9(3 cos 0 — 2)| = u(3 cos6 — 1)? 


5.21 S = vot — (gt2/2) (sen 6 + u cos 0) 
o rolamento puro acontece quando 
S = (208/9)(sen 0 + 6u cos 0) (2 sen O + Tu cos 0)? 


5.25 Vemi = —P/4m, wı = —3p/2ml 
Vem2 = 5P/4M, wz = 9p/ml 
Vp = — p/m 


